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DESCRIPTION 

Method of Simulating Indoor Behavior of Pesticidal 
Compound 
Technical Field 

The present invention relates to a method of 
simulating, when a chemical containing a pesticidal 
compound is used indoors, an indoor behavior of the 
pesticidal compound including an estimation method of 
estimating the indoor behavior of the compound and a safety- 
evaluating method of evaluating its safety in human bodies 
by using the estimation method; and, in particular, to a 
method of simulating an indoor behavior of a pesticidal 
compound in the case where a chemical containing the 
pesticidal compound is residually sprayed, sprayed in an 
indoor space, heated to vaporize, or sprayed over the whole 
floor area. 
Background Art 

Conventionally known is a f ugaciousness (hereinafter 
referred to as Fugacity) model for simulating a behavior 
of a chemical material in global environment. The 
above-mentioned fugacity model utilizes Fugacity whose unit 
is an external force by which the chemical material escapes 
from one medium to another medium, i.e., pressure. For 
example, when the chemical material concentration differs 
between two media A and B, concentrations in the respective 
media are expressed by: 
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N A /V A = f A Z A 
N B /V B = f B Z B 

wherein N is chemical mass, V is volume of medium, f is 
Fugacity, and Z is Fugacity capacity of medium. 
5 Here, while the mass N changes over time according to 

transference and degradation of the chemical material 
between the media A and B; assuming that volume V and Fugacity 
capacity Z are constant, the above-mentioned two 
expressions are represented as: 
10 (df A /dt)V A Z A = dN A /dt 

= - ( Degradation ) A ± (Transference)^ 
(df B /dt)V B Z B = dN B /dt 

= - (Degradation) B ± (Transference ) AB 
When Degradation and Transference in these two 
15 differential equations are given, unknown parameters f A and 

f B can be determined by calculation. When these parameters 
are respectively multiplied by Fugacity capacity Z A and Z B , 
the chemical material concentrations in the respective 
media in a specific period of time can be simulated. 
20 As an apparatus for simulating a behavior of a chemical 

material, Japanese Patent Application Laid-Open No. 64- 
88811 discloses a configuration of closed-space simulator 
which can perform simulation in response to any capacity 
of closed space without actually constructing a closed space 
25 when evaluating temperature change of a specific gas 

component such as carbonic acid gas. 
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As a configuration for evaluating influence of a 
harmful material on human bodies, Japanese Patent 
Application Laid-Open No. 3-89146 discloses a configuration 
of percutaneous absorption evaluation apparatus employing 
a vertical type diffusion cell which is a system closer to 
a clinical state than is a parallel type cell, thereby being 
capable of simultaneously measuring, in real time, the 
process of a chemical being emitted from its base on the 
skin by optoacoustic measurement and the process of the 
chemical infiltrating through the skin by absorptiometry. 

Also, Japanese Patent Application Laid-Open No. 7- 
218496 discloses a configuration of system which uses the 
fact that a dissolution parameter inherently existing in 
a chemical material and a logarithmic value of median lethal 
dose of the chemical material with respect to a mammal are 
in a specific correlation, estimating acute toxicity of the 
chemical material with respect to the mammal. 

Further, simulations of indoor behavior of a 
pesticidal compound when insecticides are sprayed in an 
indoor space, electrically heated to vaporize in a room, 
and sprayed over the whole floor surface are respectively 
disclosed in Y. Matoba et al. , "A SIMULATION OF INSECTICIDES 
IN INDOOR AEROSOL SPACE SPRAYING," Chemosphere f Vol.26, 
No. 6, pp. 1167-1186, 1993; Y. Matoba et al., "INDOOR 
SIMULATION OF INSECTICIDES SUPPLIED WITH AN ELECTRIC 
VAPORIZER BY THE FUGACITY MODEL , " Chemosphere , Vol.28, No . 4 , 
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pp. 767-786, 1994; and Y. Matoba et al . , "INDOOR SIMULATION 
OF INSECTICIDES IN BROADCAST SPRAYING," Chemosphftrfi , 
Vol.30, No. 2, pp. 345-365, 1995. 
Disclosure of Invention 

The above-mentioned simulation models, however, do 
not mention how to solve differential equations, and minute 
time units set when solving the differential equations are 
assumed to be constant. Theoretically, the smaller is the 
minute time unit, the longer becomes the calculation time; 
whereas the solution would not converge when the minute time 
unit is large. Accordingly, in the case where a 
differential equation containing a parameter which changes 
over time is to be solved, when the minute time unit is set 
to a constant value so that the solution does not diverge, 
there is a problem that the processing speed of a computer 
must be enhanced. 

Also, the above-mentioned simulation models fail to 
mention any security with respect to human bodies. 

In order to solve the above-mentioned conventional 
problems, it is an object of the present invention to provide 
a method of simulating an indoor behavior of a pesticidal 
compound, which can process simultaneous differential 
equations accurately in a short time by automatically 
setting a minute time unit. 

In order to achieve the above-mentioned object, the 
method of simulating an indoor behavior of a pesticidal 
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compound in accordance with the present invention comprises 
a step of dividing an indoor environment into predetermined 
media (constituents) and forming a differential equation 
concerning a fugacity of the compound in each of the media; 
5 a step of determining the fugacity of the compound in each 

of the media from the differential equation; a step of 
determining the indoor behavior of the compound from the 
fugacity of the compound in each of the media; and a step 
of changing, in response to a fluctuation in mass balance 

10 of the compound indoors, a minute time unit used when solving 

the differential equation. 

As the indoor environment is divided into 
predetermined media, and exchanges of the chemical compound 
between the media and the like are taken into account, 

15 simulation results close to the actual behavior of the 

compound can be obtained, while the minute time unit can 
be set automatically in response to fluctuation in mass 
balance when solving simultaneous differential equations 
including a parameter which changes over time . Accordingly , 

2 0 when a computer processes the above-mentioned differential 

equation, accurate solutions can be obtained in a short time. 

Preferably, the method of simulating an indoor 
behavior of a pesticidal compound in accordance with the 
present invention further comprises a step of evaluating 

2 5 safety of the compound with respect to a human body according 

to the indoor behavior of the compound. 
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As a consequence of this configuration, the safety of 
the pesticidal compound with respect to the human body can 
be evaluated accurately in a short time. Accordingly, when 
formulating a chemical such as insecticide including the 
5 above-mentioned compound, simulation can be easily repeated 

while changing conditions, thereby making it easier to 
formulate a chemical having a high safety conforming to the 
aimed object. 

Further, the method of simulating an indoor behavior 

10 of a pesticidal compound in accordance with the present 

invention may be such that the above-mentioned compound is 
introduced into an indoor space as a solution containing 
the compound is residually sprayed; whereas the above- 
mentioned media are a spraying site, suspended particles 

15 which are divided into at least one kind according to size, 

indoor air, a floor, a wall, and a ceiling. 

Preferably, in this case, the differential equation 
at the spraying site is a differential equation stating a 
relationship among temporal change of fugacity of the 

2 0 compound at the spraying site, temporal change in volume 

of the spraying site, amount of attachment of the suspended 
particles to the spraying site, amount of transference of 
the compound between the spraying site and another medium, 
and change in amount of degradation of the compound at the 

25 spraying site; the differential equation in the suspended 

particles is a differential equation stating a relationship 
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among temporal change of fugacity of the compound in the 
suspended particles, temporal change in volume of the 
suspended particles, amount of transference of the compound 
between the suspended particles and another medium, and 
5 change in amount of degradation of the compound in the 

suspended particles; the differential equation in the 
indoor air is a differential equation stating a relationship 
among temporal change of fugacity of the compound in the 
indoor air, amount of discharge of the compound outdoors, 

10 amount of transference of the compound between the indoor 

air and another medium, and change in amount of degradation 
of the compound in the indoor air; the differential equation 
at the floor is a differential equation stating a 
relationship among temporal change of fugacity of the 

15 compound at the floor, temporal change in volume of the floor, 

amount of attachment of the suspended particles to the floor, 
amount of transference of the compound between the floor 
and another medium, and change in amount of degradation of 
the compound at the floor; the differential equation at the 

20 wall is a differential equation stating a relationship among 

temporal change of fugacity of the compound at the wall, 
temporal change in volume of the wall, amount of attachment 
of the suspended particles to the wall, amount of 
transference of the compound between the wall and another 

25 medium, and change in amount of degradation of the compound 

at the wall; and the differential equation at the ceiling 
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is a differential equation stating a relationship among 
temporal change of fugacity of the compound at the ceiling, 
temporal change in volume of the ceiling, amount of 
attachment of the suspended particles to the ceiling, amount 
of transference of the compound between the ceiling and 
another medium, and change in amount of degradation of the 
compound at the ceiling. 

The method of simulating an indoor behavior of a 
pesticidal compound in accordance with the present 
invention may be such that the above-mentioned compound is 
introduced into an indoor space as a solution containing 
the compound is spatially sprayed; whereas the above- 
mentioned media are suspended particles which are divided 
into at least one kind according to size, indoor air, a floor, 
a wall, and a ceiling. 

Preferably, in this case, the differential equation 
in the suspended particles is a differential equation 
stating a relationship among temporal change of fugacity 
of the compound in the suspended particles, temporal change 
in volume of the suspended particles, amount of transference 
of the compound between the suspended particles and another 
medium, and change in amount of degradation of the compound 
in the suspended particles; the differential equation in 
the indoor air is a differential equation stating a 
relationship among temporal change of fugacity of the 
compound in the indoor air, amount of discharge of the 
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compound outdoors, amount of transference of the compound 
between the indoor air and another medium, and change in 
amount of degradation of the compound in the indoor air; 
the differential equation at the floor is a differential 
5 equation stating a relationship among temporal change of 

fugacity of the compound at the floor, temporal change in 
volume of the floor, amount of attachment of the suspended 
particles to the floor, amount of transference of the 
Jf! compound between the floor and another medium, and change 

J=i 10 in amount of degradation of the compound at the floor; the 

m differential equation at the wall is a differential equation 

J5J stating a relationship among temporal change of fugacity 

]U of the compound at the wall, temporal change in volume of 

the wall, amount of attachment of the suspended particles 
"1" 15 to the wall, amount of transference of the compound between 

the wall and another medium, and change in amount of 
degradation of the compound at the wall ; and the differential 
equation at the ceiling is a differential equation stating 
a relationship among temporal change of fugacity of the 
20 compound at the ceiling, temporal change in volume of the 

ceiling, amount of attachment of the suspended particles 
to the ceiling, amount of transference of the compound 
between the ceiling and another medium, and change in amount 
of degradation of the compound at the ceiling. 
25 The method of simulating an indoor behavior of a 

pesticidal compound in accordance with the present 
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invention may be such that the above-mentioned compound is 
introduced into an indoor space as a solution containing 
the compound is heated to vaporize; whereas the above- 
mentioned media are condensed particles which are divided 
5 into at least one kind according to generation and extinction, 

high-concentration air, medium-concentration air, low- 
concentration air, a floor, a wall, and a ceiling which is 
divided into at least one kind according to compound 
concentration . 

10 Preferably, in this case, the differential equation 

in the condensed particles is a differential equation 
stating a relationship among temporal change of fugacity 
of the compound in the condensed particles, temporal change 
in volume of the condensed particles, amount of transference 

15 of the compound between the condensed particles and another 

medium, and change in amount of degradation of the compound 
in the condensed particles; the differential equation in 
the high-concentration air is a differential equation 
stating a relationship among temporal change of fugacity 

2 0 of the compound in the high-concentration air, amount of 

discharge of the compound, amount of transference of the 
compound between the high-concentration air and another 
medium, and change in amount of degradation of the compound 
in the high-concentration air; the differential equation 

25 in the medium-concentration air is a differential equation 

stating a relationship among temporal change of fugacity 
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of the compound in the medium-concentration air, amount of 
transference of the compound between the medium- 
concentration air and another medium, and change in amount 
of degradation of the compound in the medium-concentration 
air; the differential equation in the low-concentration air 
is a differential equation stating a relationship among 
temporal change of fugacity of the compound in the low- 
concentration air, amount of discharge of the compound 
outdoors, amount of transference of the compound between 
the low-concentration air and another medium, and change 
in amount of degradation of the compound in the low- 
concentration air; the differential equation at the floor 
is a differential equation stating a relationship among 
temporal change of fugacity of the compound at the floor, 
temporal change in volume of the floor, amount of 
transference of the compound between the floor and another 
medium, and change in amount of degradation of the compound 
at the floor; the differential equation at the wall is a 
differential equation stating a relationship among temporal 
change of fugacity of the compound at the wall, temporal 
change in volume of the wall, amount of transference of the 
compound between the wall and another medium, and change 
in amount of degradation of the compound at the wall; and 
the differential equation at the ceiling is a differential 
equation stating a relationship among temporal change of 
fugacity of the compound at the ceiling, temporal change 
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in volume of the ceiling, amount of transference of the 
compound between the ceiling and another medium, and change 
in amount of degradation of the compound at the ceiling. 

The method of simulating an indoor behavior of a 
pesticidal compound in accordance with the present 
invention may be such that the above-mentioned compound is 
introduced into an indoor space as a solution containing 
the compound is sprayed over the whole floor; whereas the 
above-mentioned media are suspended particles which are 
divided into at least one kind according to size, indoor 
air, a floor, a wall, and a ceiling. 

Preferably, in this case, the differential equation 
in the suspended particles is a differential equation 
stating a relationship among temporal change of fugacity 
of the compound in the suspended particles, temporal change 
in volume of the suspended particles , amount of transference 
of the compound between the suspended particles and another 
medium, and change in amount of degradation of the compound 
in the suspended particles; the differential equation in 
the indoor air is a differential equation stating a 
relationship among temporal change of fugacity of the 
compound in the indoor air, amount of discharge of the 
compound outdoors, amount of transference of the compound 
between the indoor air and another medium, and change in 
amount of degradation of the compound in the indoor air; 
the differential equation at the floor is a differential 
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equation stating a relationship among temporal change of 
fugacity of the compound at the floor, temporal change in 
volume of the floor, amount of attachment of the suspended 
particles to the floor, amount of transference of the 
5 compound between the floor and another medium, and change 

in amount of degradation of the compound at the floor; the 
differential equation at the wall is a differential equation 
stating a relationship among temporal change of fugacity 
of the compound at the wall, temporal change in volume of 

10 the wall, amount of attachment of the suspended particles 

to the wall, amount of transference of the compound between 
the wall and another medium, and change in amount of 
degradation of the compound at the wall ; and the differential 
equation at the ceiling is a differential equation stating 

15 a relationship among temporal change of fugacity of the 

compound at the ceiling, temporal change in volume of the 
ceiling, amount of attachment of the suspended particles 
to the ceiling, amount of transference of the compound 
between the ceiling and another medium, and change in amount 

20 of degradation of the compound at the ceiling. 

The method of simulating an indoor behavior of a 
pesticidal compound in accordance with the present 
invention may be such that the floor is constituted by a 
rug having ears of fiber, whereas a space between the ears 

25 is added to the above-mentioned media. 

Preferably, in this case, the differential equation 
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in the space between the ears is a differential equation 
stating a relationship among temporal change of fugacity 
of the compound in the space between the ears, temporal 
change in volume of the solution containing the compound 
in the space between the ears, amount of attachment of the 
compound into the space portion between the ears by falling, 
amount of transference of the compound between the space 
portion between the ears and another medium, and change in 
amount of degradation of the compound in the space portion 
between the ears . 

Even in the case where a rug having ears of fiber is 
spread on the floor, when the space between the ears is 
further added to the media, the behavior of the compound 
can be simulated accurately, thus allowing various kinds 
of simulations to be performed. 

In order to achieve the above-mentioned object, the 
computer program product of the present invention is a 
computer program product to be used together with an 
information processing apparatus comprising input means for 
receiving a data input from outside, display means, and 
readout means for reading out information from a 
computer-usable storage medium; the computer program 
product comprising a computer-usable storage medium which 
has a program area for storing a program and has a 
computer-readable program materialized in the storage 
medium for causing, according to data input from the input 
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means, the display means to display a result of simulation 
of an indoor behavior of a pesticidal compound; the computer 
program product comprising, in the program area, a program 
for dividing an indoor environment into predetermined media 
5 and forming a differential equation concerning a fugacity 

of the compound, a program for determining the fugacity of 
the compound in each of the media from the differential 
equation; a program for determining the indoor behavior of 
the compound from the fugacity of the compound in each of 

10 the media, and a program for changing, in response to a 

fluctuation in mass balance of the compound indoors , a minute 
time unit used when solving the differential equation. 

As the indoor environment is divided into 
predetermined media, and exchanges of the chemical compound 

15 between the media and the like are taken into account, 

simulation results close to the actual behavior of the 
compound can be obtained, while the minute time unit can 
be set automatically in response to fluctuation in mass 
balance when solving simultaneous differential equations 

2 0 including a parameter which changes over time. Accordingly, 

when a computer processes the above-mentioned differential 
equation, accurate solutions can be obtained in a short time. 

Preferably, the computer program product further 
comprises, in the program area, a program for evaluating 

25 safety of the compound with respect to a human body according 

to the indoor behavior of the compound. 
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As a consequence of this configuration, the safety of 
the pesticidal compound with respect to the human body can 
be evaluated accurately in a short time. Accordingly, when 
formulating a chemical such as insecticide including the 
above-mentioned compound, simulation can be easily repeated 
while changing conditions, thereby making it easier to 
formulate a chemical having a high safety conforming to the 
aimed object. 

The computer program product of the present invention 
may be such that the above-mentioned compound is introduced 
into an indoor space as a solution containing the compound 
is residually sprayed; whereas the above-mentioned media 
are a spraying site, suspended particles which are divided 
into at least one kind according to size, indoor air, a floor, 
a wall, and a ceiling. 

Preferably, in this case, the differential equation 
at the spraying site is a differential equation stating a 
relationship among temporal change of fugacity of the 
compound at the spraying site, temporal change in volume 
of the spraying site, amount of attachment of the suspended 
particles to the spraying site, amount of transference of 
the compound between the spraying site and another medium, 
and change in amount of degradation of the compound at the 
spraying site; the differential equation in the suspended 
particles is a differential equation stating a relationship 
among temporal change of fugacity of the compound in the 
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suspended particles, temporal change in volume of the 
suspended particles, amount of transference of the compound 
between the suspended particles and another medium, and 
change in amount of degradation of the compound in the 
suspended particles; the differential equation in the 
indoor air is a differential equation stating a relationship 
among temporal change of fugacity of the compound in the 
indoor air, amount of discharge of the compound outdoors, 
amount of transference of the compound between the indoor 
air and another medium, and change in amount of degradation 
of the compound in the indoor air; the differential equation 
at the floor is a differential equation stating a 
relationship among temporal change of fugacity of the 
compound at the floor , temporal change in volume of the floor , 
amount of attachment of the suspended particles to the floor, 
amount of transference of the compound between the floor 
and another medium, and change in amount of degradation of 
the compound at the floor; the differential equation at the 
wall is a differential equation stating a relationship among 
temporal change of fugacity of the compound at the wall, 
temporal change in volume of the wall, amount of attachment 
of the suspended particles to the wall, amount of 
transference of the compound between the wall and another 
medium, and change in amount of degradation of the compound 
at the wall; and the differential equation at the ceiling 
is a differential equation stating a relationship among 
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temporal change of fugacity of the compound at the ceiling, 
temporal change in volume of the ceiling, amount of 
attachment of the suspended particles to the ceiling, amount 
of transference of the compound between the ceiling and 
another medium, and change in amount of degradation of the 
compound at the ceiling. 

The computer program product of the present invention 
may be such that the above-mentioned compound is introduced 
into an indoor space as a solution containing the compound 
is spatially sprayed; whereas the above-mentioned media are 
suspended particles which are divided into at least one kind 
according to size, indoor air, a floor, a wall, and a ceiling. 

Preferably, in this case, the differential equation 
in the suspended particles is a differential equation 
stating a relationship among temporal change of fugacity 
of the compound in the suspended particles, temporal change 
in volume of the suspended particles, amount of transference 
of the compound between the suspended particles and another 
medium, and change in amount of degradation of the compound 
in the suspended particles; the differential equation in 
the indoor air is a differential equation stating a 
relationship among temporal change of fugacity of the 
compound in the indoor air, amount of discharge of the 
compound outdoors, amount of transference of the compound 
between the indoor air and another medium, and change in 
amount of degradation of the compound in the indoor air; 
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the differential equation at the floor is a differential 
equation stating a relationship among temporal change of 
fugacity of the compound at the floor, temporal change in 
volume of the floor, amount of attachment of the suspended 
particles to the floor, amount of transference of the 
compound between the floor and another medium, and change 
in amount of degradation of the compound at the floor; the 
differential equation at the wall is a differential equation 
stating a relationship among temporal change of fugacity 
of the compound at the wall, temporal change in volume of 
the wall, amount of attachment of the suspended particles 
to the wall, amount of transference of the compound between 
the wall and another medium, and change in amount of 
degradation of the compound at the wall ; and the differential 
equation at the ceiling is a differential equation stating 
a relationship among temporal change of fugacity of the 
compound at the ceiling, temporal change in volume of the 
ceiling, amount of attachment of the suspended particles 
to the ceiling, amount of transference of the compound 
between the ceiling and another medium, and change in amount 
of degradation of the compound at the ceiling* 

The computer program product of the present invention 
may be such that the above-mentioned compound is introduced 
into an indoor space as a solution containing the compound 
is heated to vaporize; whereas the above-mentioned media 
are condensed particles which are divided into at least one 
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kind according to generation and extinction, high- 
concentration air, medium-concentration air, low- 
concentration air, a floor, a wall, and a ceiling which is 
divided into at least one kind according to compound 
concentration . 

Preferably, in this case, the differential equation 
in the condensed particles is a differential equation 
stating a relationship among temporal change of fugacity 
of the compound in the condensed particles, temporal change 
in volume of the condensed particles, amount of transference 
of the compound between the condensed particles and another 
medium, and change in amount of degradation of the compound 
in the condensed particles; the differential equation in 
the high-concentration air is a differential equation 
stating a relationship among temporal change of fugacity 
of the compound in the high-concentration air, amount of 
discharge of the compound, amount of transference of the 
compound between the high-concentration air and another 
medium, and change in amount of degradation of the compound 
in the high-concentration air; the differential equation 
in the medium-concentration air is a differential equation 
stating a relationship among temporal change of fugacity 
of the compound in the medium-concentration air, amount of 
transference of the compound between the medium- 
concentration air and another medium, and change in amount 
of degradation of the compound in the medium-concentration 
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air; the differential equation in the low-concentration air 
is a differential equation stating a relationship among 
temporal change of fugacity of the compound in the low- 
concentration air, amount of discharge of the compound 
5 outdoors, amount of transference of the compound between 

the low-concentration air and another medium, and change 
in amount of degradation of the compound in the low- 
concentration air; the differential equation at the floor 
is a differential equation stating a relationship among 

10 temporal change of fugacity of the compound at the floor, 

temporal change in volume of the floor, amount of 
transference of the compound between the floor and another 
medium, and change in amount of degradation of the compound 
at the floor; the differential equation at the wall is a 

15 differential equation stating a relationship among temporal 

change of fugacity of the compound at the wall, temporal 
change in volume of the wall, amount of transference of the 
compound between the wall and another medium, and change 
in amount of degradation of the compound at the wall; and 

20 the differential equation at the ceiling is a differential 

equation stating a relationship among temporal change of 
fugacity of the compound at the ceiling, temporal change 
in volume of the ceiling, amount of transference of the 
compound between the ceiling and another medium, and change 

25 in amount of degradation of the compound at the ceiling. 

The computer program product of the present invention 
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may be such that the above-mentioned compound is introduced 
into an indoor space as a solution containing the compound 
is sprayed over the whole floor; while the above-mentioned 
media are suspended particles which are divided into at least 
5 one kind according to size, indoor air, a floor, a wall, 

and a ceiling. 

Preferably, in this case, the differential equation 
in the suspended particles is a differential equation 
bf) stating a relationship among temporal change of fugacity 

jrj 10 of the compound in the suspended particles, temporal change 

l= L { in volume of the suspended particles, amount of transference 

Lip 

m. of the compound between the suspended particles and another 

r-i medium, and change in amount of degradation of the compound 

M in the suspended particles; the differential equation in 

15 the indoor air is a differential equation stating a 

* A \ relationship among temporal change of fugacity of the 

compound in the indoor air, amount of discharge of the 
compound outdoors, amount of transference of the compound 
between the indoor air and another medium, and change in 
2 0 amount of degradation of the compound in the indoor air; 

the differential equation at the floor is a differential 
equation stating a relationship among temporal change of 
fugacity of the compound at the floor, temporal change in 
volume of the floor, amount of attachment of the suspended 
25 particles to the floor, amount of transference of the 

compound between the floor and another medium, and change 
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in amount of degradation of the compound at the floor; the 
differential equation at the wall is a differential equation 
stating a relationship among temporal change of fugacity 
of the compound at the wall, temporal change in volume of 
5 the wall, amount of attachment of the suspended particles 

to the wall, amount of transference of the compound between 
the wall and another medium, and change in amount of 
degradation of the compound at the wall ; and the differential 
equation at the ceiling is a differential equation stating 

10 a relationship among temporal change of fugacity of the 

compound at the ceiling, temporal change in volume of the 
ceiling, amount of attachment of the suspended particles 
to the ceiling, amount of transference of the compound 
between the ceiling and another medium, and change in amount 

15 of degradation of the compound at the ceiling. 

The computer program product of the present invention 
may be such that the floor is constituted by a rug having 
ears of fiber, whereas a space between the ears is added 
to the above-mentioned media. 

20 Preferably, in this case, the differential equation 

in the space between the ears is a differential equation 
stating a relationship among temporal change of fugacity 
of the compound in the space between the ears, temporal 
change in volume of the solution containing the compound 

2 5 in the space between the ears, amount of attachment of the 

compound into the space portion between the ears by falling, 
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amount of transference of the compound between the space 
portion between the ears and another medium, and change in 
amount of degradation of the compound in the space portion 
between the ears . 

Even in the case where a rug having ears of fiber is 
spread on the floor, when the space between the ears is 
further added to the media, the behavior of the compound 
can be simulated accurately, thus allowing various kinds 
of simulations to be performed. 
I 10 Brief Description of Drawings 

| Fig, 1 is a flowchart showing a method of simulating 

an indoor behavior of a pesticidal compound in accordance 
i with an embodiment of the present invention; 

I Fig. 2 is an explanatory view showing a model of indoor 

s 15 environment in the case where an indoor behavior of a 

? 

: 

compound is performed in the first embodiment of the present 
invention; 

Fig. 3A is a graph showing temporal change of 
concentration in indoor air; 
2 0 Fig. 3B is an explanatory view in the case where mean 

concentration in indoor air is determined by use of the graph 
in Fig. 3A; 

Fig. 4A is a graph showing temporal change in residual 
amount ; 

2 5 Fig. 4B is an explanatory view in the case where mean 

residual amount is determined by use of the graph in Fig. 
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Fig. 5 is a flowchart showing a method of setting an 
estimation nick time width; 

Fig. 6A is an explanatory view showing a modeled state 
of a carpet in the case where an indoor behavior of a 
pesticidal compound is estimated in accordance with the 
second embodiment of the present invention; 

Fig. 6B is an enlarged view of a carpet fiber in the 
case where an indoor behavior of a pesticidal compound is 
estimated in accordance with the second embodiment of the 
present invention ; 

Fig. 6C is an enlarged view of a space portion of the 
carpet in the case where an indoor behavior of a pesticidal 
compound is estimated in accordance with the second 
embodiment of the present invention; 

Fig. 7 is a graph showing temporal change of 
concentration in the air in Example 1 in the above-mentioned 
first embodiment; 

Fig. 8A is a graph showing temporal change in 
concentration in the air in the case where ventilation rate 
is changed in Example 2; 

Fig. 8B is a graph showing temporal change in 
concentration in the air in the case where a window is opened 
and closed in Example 3; 

Fig. 9A is a graph showing temporal change in residual 
amount on the floor in Example 1; 
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Fig. 9B is a graph showing temporal change in residual 
amount on the wall and ceiling in Example 1; 

Fig. 10A is a graph comparing actually measured values 
of concentration in the air and their simulation results 
in Example 1; 

Fig. 1 OB is a graph comparing actually measured values 
of residual amount at the spraying site, floor, wall, and 
ceiling and their simulation results in Example 1; 

Fig. 11 is an explanatory view showing a modeled state 
of indoor environment in the case where an indoor behavior 
of a compound is estimated in the third embodiment of the 
present invention ; 

Fig. 12A is a graph comparing actually measured values 
of concentration in the "air and their simulation results 
in Examples 6 to 10; 

Fig. 12B is a graph comparing actually measured values 
of concentration in the air and their simulation results 
in Examples 6 to 10; 

Fig. 12C is a graph comparing actually measured values 
of residual amount on the floor and ceiling and their 
simulation results in Example 6; 

Fig. 13A is an explanatory view showing a modeled state 
of a carpet in the case where an indoor behavior of a 
pesticidal compound is estimated in accordance with the 
fourth embodiment of the present invention; 

Fig. 13B is an enlarged view of a carpet fiber in the 
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case where an indoor behavior of a pesticidal compound is 
estimated in accordance with the fourth embodiment of the 
present invention ; 

Fig. 13C is an enlarged view of a space portion of the 
carpet in the case where an indoor behavior of a pesticidal 
compound is estimated in accordance with the fourth 
embodiment of the present invention; 

Fig. 14 is an explanatory view showing an example of 
size of a typical room, together with an airflow; 

Fig. 15 is an explanatory view showing a modeled state 
of indoor environment in the case where an indoor behavior 
of a compound is estimated in the fifth embodiment of the 
present invention ; 

Fig. 16A is an explanatory view showing a modeled state 
of indoor environment in the case where an indoor behavior 
of the compound is estimated in the fifth embodiment of the 
present invention ; 

Fig. 16B is an enlarged view of an electric heating 
evaporator in the fifth embodiment of the present invention; 

Fig. 17 is an explanatory view showing a modeled state 
of indoor environment in the case where an indoor behavior 
of a compound is estimated in the seventh embodiment of the 
present invention ; 

Fig. 18 is an explanatory view showing a falling state 
of a particle zone; 

Fig. 19 is a conf igurational view of a storage medium 
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20 employed in each of the above-mentioned embodiments of 
the present invention; 

Fig. 20 is a conf igurational view of an information 
processing apparatus for operating a program stored in the 
storage medium 20; and 

Fig. 21 is a conf igurational view of the information 
processing apparatus for operating a program stored in the 
storage medium 20. 

Best Modes for Carrying Out the Invention 
( 1 ) First Embodiment 

With reference to the drawings, the first embodiment 
of the present invention will be explained in the following. 

The method of simulating an indoor behavior of a 
pesticidal compound (hereinafter "simply referred to as 
compound) contained in a chemical (solution) such as 
insecticide in accordance with this embodiment estimates, 
for example, the behavior of the compound when the chemical 
is residually sprayed in order to exterminate indoor vermin 
such as cockroach. Here, residual spraying refers to a 
process in which a chemical is sprayed to a local site such 
as the boundary between a floor and a wall. The chemical 
encompasses all kinds of insecticides and the like, 
including pyrethroid insecticidal compounds, 
organophosphorus compounds, carbamate compounds, and 
insect growth restrainers (IGR). 



First, as shown in Fig. 2, an environment is classified 
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into seven kinds of media consisting of spraying site 1, 
large particles (suspended particles) 2, small particles 
(suspended particles) 3, whole indoor air 4, floor 5, wall 
6, and ceiling 7. 
5 When a chemical is sprayed to the boundary between the 

floor 5 and the wall 6, part of the chemical attaches to 
the spraying site 1, whereas the rest dissipates into the 
air as large and small particles 2 and 3, The chemical 
attached to the spraying site 1 not only infiltrates into 

10 this site but also evaporates. On the other hand, among the 

large and small particles 2 and 3 , the smaller particles 
3, due to their weight, scatter higher than the larger 
particles 2 . When these particles attach onto interior 
materials of the floor 5, wall 6, and ceiling 7, the compound 

15 in the particles moves into the air 4 after infiltrating 

into the interior materials or is degraded by light or 
oxidation. Such large and small particles 2 and 3 and 
evaporated chemical have decreased due to ventilation of 
the room. 

20 Simultaneous differential equations concerning 

Fugacity of the compound in the above-mentioned seven kinds 
of media are formed and are solved by Runge-Kutta-Gill method, 
whereby the indoor behavior of the compound in the chemical 
when the chemical is residually sprayed, i.e., compound 

2 5 distribution in the room, and temporal change in amount of 

attachment to the spraying site 1 , floor 5 , wall 6 , or ceiling 
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7 are estimated and analyzed. 

Here, since the floor 5 mainly comprises three kinds, 
i.e., tatami, flooring, or carpet (rug), the differential 
equation is formed so as to correspond to the kind of floor. 
Explained in this embodiment is the case where the floor 
5 is made of flooring. 

With reference to Fig. 1, a method of simulating an 
indoor behavior of the compound will be explained. The 
j simulation method in accordance with this embodiment can 

10 be mainly divided into a step of dividing an indoor 

environment into predetermined media and forming a 
ijl differential equation concerning a fugacity of the compound 

p in each medium (SI to S2) ; a step of determining the fugacity 

of the compound in each medium from the dif ferential equation 
15 (S3); a step of determining the indoor behavior of the 

111 

" compound from the fugacity of the compound in each medium 

(S4); a step of changing, in response to a fluctuation in 
mass balance of the compound indoors, a minute time unit 
used when solving the differential equation (S21 to S35, 

20 see Fig. 5) ; and a step of evaluating safety of the compound 

with respect to a human body according to the indoor behavior 
of the compound (S5 to S17). 

First, at step (hereinafter abridged as S ) 1, a primary 
condition is inputted. The primary condition comprises 

25 physicochemical properties of the compound (see Table 1), 

indoor environmental behavior properties of the compound 
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(see Table 2 ) , indoor environment (see Table 3 ) , and product 
properties (see Table 4). 
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TABLE 1 



PHYSICOCHEMICAL 
PROPERTY OF COMPOUND 


SYMBOL 


EXAMPLE OF VALUE 


MOLECULAR WEIGHT 

SPECIFIC GRAVITY 

VAPOR PRESSURE 

WATER SOLUBILITY 

OCTANOL/ WATER 
DISTRIBUTION COEFFICIENT 


P 
PS 
CS 
Kqw 


350.46 Cg/mole] 
1.11 Cg/cm3] 
1.60 XI 0" 4 [Pa] 
6.28 XI 0" 3 Cmole/m3] 
7.59 X10 4 



TABLE 2 



INDOOR ENVIRONMENTAL 

BEHAVIOR 
PROPERTY OF COMPOUND 


SYMBOL 


EXAMPLE OF VALUE 


DEGRADATION HALF-LIFE 
IN SPRAYING SITE 
PARTICLE 
AIR 
FLOOR 

WALL / CEILING 


Ti(i=2,3) 
tj 0=6,7) 


13 Cday] 
7.7 Ch] 
7.7 Ch] 

13 Cday] (IN CASE OF FLOORING) 
31 Cday] (IN CASE OF WALLPAPER) 
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Subsequently, using the primary condition, a 
secondary condition is determined by calculation (S2) . The 
secondary condition comprises a temporally unchangeable 
definite factor condition (see Table 5) determined by the 
5 primary condition alone, and a temporal change condition 

accompanying a temporal change (see Table 6), Their 
calculations will be explained later in detail. 
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TABLE 6 



TEMPORAL CHANGE CONDITION 


SYMBOL 


SUSPENDED PARTICLE 

DIAMETER 
FALLING SPEED 
FLOATING NUMBER 
FLOATING HEIGHT 


di 0=2,3) 
vi (i=2,3) 
ni (i=2,3) 
Hi (i=2,3) 


TRANSFERENCE SPEED OF COMPOUND 
IN SPRAYING SITE 

SUSPENDED PARTICLE 
AIR 

FLOOR / WALL / CEILING 


kl 

ki 0=2,3) 
k 4 

kj 0=5,6,7) 


TRANSFERENCE COEFFICIENT OF COMPOUND 

BETWEEN SPRAYING SITE AND AIR 

TRANSFERENCE COEFFICIENT OF COMPOUND 

BETWEEN SUSPENDED PARTICLE AND AIR 

TRANSFERENCE COEFFICIENT OF COMPOUND 

BETWEEN AIR AND FLOOR / WALL / CEILING 


D14 

D i4 0=2,3) 
D 4 j 0=5,6,7) 


VOLUME OF SPRAYING SITE 
VOLUME OF SUSPENDED PARTICLE 
VOLUME OF FLOOR / WALL / CEILING 


Vl 

Vj 0=2,3) 
Vj 0=5,6,7) 
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Using the secondary condition, seven kinds of Fugacity 
are calculated (S3). Namely, seven kinds of differential 
equations concerning the spraying site 1, large and small 
particles 2 and 3, air 4, floor 5, wall 6, and ceiling 7 
5 are simultaneously formed and are solved by Runge- 

Kutta-Gill method, whereby seven kinds of Fugacity are 
computed over time. Here, an estimation nick time width 
(minute time unit) set when solving the differential 
equations is automatically set so as to be varied in response 

10 to a fluctuation in mass balance of the compound. 

At S4, using thus computed Fugacity (f 3/ f 4 ) of the 
small particles 3 and air 4, together with Fugacity ( f 2 ) of 
the large particles 2 when necessary, temporal 
concentration of the compound in the indoor air is computed; 

15 whereas temporal residual amounts of the compound to the 

spraying site 1 and floor 5 are computed by use of Fugacity 
(f lf f 5 ) of the spraying site 1 and floor 5. 

At S5, it is judged whether or not to perform safety 
evaluation in the case where the chemical is inhaled. If 

2 0 the safety evaluation in inhalation is to be performed, then 

an estimated exposure amount in inhalation indicating a 
degree of exposure upon inhalation of the contaminated air 
is computed by use of the above-mentioned temporal 
concentration in the indoor air (S6). Thereafter, an 

25 inhalation safety coefficient is computed according to the 

estimated exposure amount in inhalation (S7). At S8, the 
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inhalation safety coefficient is compared with a reference 
value defined in each country. If the inhalation safety 
coefficient exceeds the reference value, it is judged that 
"there is no problem in safety." By contrast, if the 
inhalation safety coefficient is lower than the reference 
value at S8, it is judged that "there is a problem in safety, " 
and the operation returns to SI, where alteration of the 
primary condition such as alteration of compound, 
alteration of chemical formulation, alternation of using 
condition, or the like is considered. 

If the safety evaluation in inhalation is not selected 
at S5, calculated is an estimated amount of percutaneous 
exposure indicating to what extent the skin is exposed in 
contact with the spraying site 1 or floor 5 to which the 
chemical is attached (S9) . Thereafter, at S10, it is judged 
whether to perform a percutaneous safety evaluation or not. 
If the percutaneous safety evaluation is to be performed, 
a percutaneous safety coefficient is computed according to 
the estimated percutaneous exposure amount (Sll). At S12, 
as with the safety evaluation in inhalation, the 
percutaneous safety coefficient is compared with a 
reference value defined in each country, whereby the safety 
is evaluated. If there is a problem in safety at S12, the 
operation returns to SI, where alteration of the primary 
condition is considered. 

If the percutaneous safety evaluation is not selected 
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at S10, calculated according to the estimated percutaneous 
exposure amount is an estimated oral exposure amount 
indicating the degree of exposure in the case where the 
chemical attached to a hand or the like is taken orally ( S13 ) . 
5 Subsequently, according to the estimated oral exposure 

amount, an oral safety coefficient is computed (S14). This 
0 oral exposure may occur, in particular, when an infant puts 

a chemical-attached hand into the mouth. At S15, as with 

o 

in the safety evaluation in inhalation, the oral safety 

q 10 coefficient is compared with a reference value defined in 

jry each country, whereby the safety is evaluated. If there is 

rfi a problem in safety at SI 5, the operation returns to SI, 

where alteration of the primary condition is considered. 
5=4 Finally, if it is judged to be safe at S8, S12, or S15, 

if* 

% 15 overall safety is evaluated. Here, the sum of respective 

! '~ reciprocals of the previously determined inhalation safety 

coefficient, percutaneous safety coefficient, and oral 
safety coefficient is determined, and the reciprocal of thus 
determined value is defined as an overall safety coefficient 
20 (S16). This overall safety coefficient is compared with a 

reference value, whereby an evaluation similar to the 
previous safety evaluation is effected (S17). 

In the following, the above-mentioned steps of S2 to 
S4, S6, S7, S9, SI 1 , S13, and S14 will be explained in detail . 
25 (i) Secondary Condition Calculation (S2) 

Calculations of the definite factor condition shown 
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in Table 5 and the temporal change condition shown in Table 
6 will be explained. 

According to surface temperature (T d ) of particles, 
properties of particle-constituting solvent (vapor 
pressure P d/ molecular weight Mj, and specific gravity p 
d ), and indoor environment (air diffusion coefficient D air/ 
vapor pressure poo, and room temperature T°o), the 
evaporation constant ( a ) of the large and small particles 
2 and 3 is defined as follows: 



wherein R is a gas constant. 

Here, the upper parts of the above-mentioned Tables 
3 and 4 exemplify the case where the residually sprayed 
aerosol is of oil base. In this case, Too (temperature at 
the site far from particles) and T d are set to room 
temperature (T), whereas poo (vapor pressure of oil at the 
site far from particles) is set to zero. 

In the case where the residually sprayed aerosol is 
of water base, it is necessary for P d , Mj, p d , T d , and P 
oo in expression (1) to be changed to values based on 
properties of water. Specifically, the above-mentioned 
factor values when the room temperature (T) is 298 [K] ( = 
25 °C) and the relative humidity is 60[%RH] are exemplified 
in the lower parts of Tables 3 and 4 . 

Here, T d can be determined, according to the room 
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temperature (T) and relative humidity (H) , from "mass-based 
humidity table" disclosed in "Kagaku Kikai no Riron to Keisan 
(Theory and Calculation of Chemical Machines)" (Second 
Edition) (Saburo Kamei ed., Sangyo Tosho) or the like. On 
the other hand, P d and P°o can be computed from the following 
expressions : 

logio^a = 10.23 - 1750/(T d -38) 
log 10 poo = 8,23 - 1750/(T-38) + log lo 0 
wherein 0 is an indoor relative humidity. 

From the vapor pressure (P 3 )/ water solubility (C s ), 
and octanol/water distribution coefficient (K ow ) of the 
compound, and from the oil component content ( p x ) of the 
floor or wall constituting the spraying site 1, Fugacity 
capacity (Z x ) of the spraying site 1 is expressed as: 



From the vapor pressure (P s ) of the compound and the 
surf ace temperature (T d ) of particles, Fugacity capacity (Z if 
i = 2 , 3 ) of the large and small particles 2 and 3 is expressed 
as : 



Though the vapor pressure of the compound in a liquid state 
(P L 3 ) is computed by use of the melting point T M of the harmful 
material and T d here; in the case where T M cannot be obtained, 
P L 3 may be set identical to the vapor pressure of the compound 




s /p s 



...(2) 
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in a solid state (P s ). 

From the room temperature (T) , Fugacity capacity (Z 4 ) 
of the air 4 is expressed as follows: 

Z.-± ...(4) 

5 From the vapor pressure (P s ), water solubility (C s ), 

and octanol/water distribution coefficient (K ow ) of the 
compound, and from the particle-constituting solvent 
component content ( Pj, j = 5, 6, 7) of the material 
constituting the floor 5, wall 6, and ceiling 7, Fugacity 

10 capacity ( Z j , j = 5 , 6 , 7 ) of the floor 5 , wall 6 , and ceiling 

7 is expressed as: „ 
Zj-PjK^CS/P 3 ...(5) 

Though the ratio ( 7 j) of the large and small particles 
2 and 3 attaching to the floor 5, wall 6, and ceiling 7 may 

15 be assumed to be such that 100% of them attach to the floor 

5 from viewpoint of safety evaluation, it may be set, for 
example, floor : wall : ceiling = 97:2:1 in view of the 
actually measured values in the past. 

The ratio (7J of the large and small particles 2 and 

20 3 falling to the spraying site 1 is determined by the 

above-mentioned ratio (7j)/ and the area (A x ) occupied by 
the spraying site 1 with respect to the size of the room 
<V«). 

The diffusion coefficient (D c ) of the compound in the 
25 floor 5, wall 6, and ceiling 7 can be computed as follows. 



43 



• 



SCC 97-06 



Q 10 

Hie 



m 15 



Namely, D c can be determined when the diffusion coefficient 
(D ca ) of the compound in the air is multiplied by 10~ 9 , whereas 
D ca can be computed from properties of the compound (e.g., 
structural formula, molecular weight, and the like) 
according to Wike and Lee method ( "Handbook of Chemical 
Property Estimation Methods," McGraw-Hill Book Company, 
1982), for example. 

Mainly generated in the seven kinds of media are 
degradation reactions caused by light and oxidation. The 
degradation constant (K x ) of the spraying site 1, 
degradation constant (Ki, i = 2, 3) of the large and small 
particles 2 and 3, degradation constant (K 4 ) of the air 4, 
degradation constant (K 5 ) of the floor 5, and degradation 
constant (K-, j = 6, 7) of the wall 6 and ceiling 7 are defined 
by use of degradation half-life (T 1 ), (TJ, {T*)r ( T $)t 
or ( T j ) as : 



K 5 =0.693 /r 5 . 
Kj = 0.693 /t. 

Here, among the degradation half -lives shown in the 
above-mentioned Table 2, it is difficult to actually measure 
the degradation half-life (T 4 ) of the compound in the air, 
whereby it is not usually easy for the degradation constant 
(K 4 ) to be determined. Therefore, when the actually 
measured value is unavailable, it may be determined as being 



K v =0.693 lx v 
K t - 0.693 /t. 
/C 4 = 0.693 /t 4 



...(6) 
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computed by Atmospheric Oxidation program (Atkinson et al. , 
1984 , Chem. Rev. Vol, 84, pp. 437-470) using the structural 
formula of the compound. Also, the degradation constant in 
the particles and in the interior materials (floor 5, wall 
5 6, and ceiling 7) may be set to zero in view of the safety 

of inhabitants when the actually measured values are 
unavailable. 

Since the volume of particles decreases over time due 
I to evaporation of the particle-constituting solvent 

I 10 component, during the evaporation of particle-constituting 

| solvent component, using the diameter (d oi , i = 2, 3) of 

particles immediately after spraying (t = 0) and the 
i evaporation constant (a) determined from expression (1), 

the diameter (d if i = 2, 3) of the large and small particles 
15 2 and 3 at time t is represented as shown in the following 

expression ( 7 ) : 

d^^dl-lat ...(7) 
4 =#7 d 0i ...(7') 

By contrast, the diameter (d t ) of particles after 
complete evaporation of the particle-constituting solvent 
20 component is represented as shown in the above expression 

(7* ) by use of the ratio of volume of the compound in the 
product (R a ) and the diameter (d oi ) of particles immediately 
after spraying. 

From the compound content (C a ) and specific gravity 
25 ( p ) of the compound, the ratio of volume of the compound 
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in the product (R a ) is expressed as Cjp. 

Since movement of the particles is under the control 
of gravity and air resistance, according to Stokes 1 law using 
the specific gravity ( p d ) of the particle-constituting 
solvent and the diameter (di) of particles determined by 
expression ( 7 ) , the falling speed ( v if i = 2 , 3 ) of the large 
and small particles 2 and 3 before complete evaporation of 
the particle-constituting solvent component is represented 
by the following expression (8) : 

V-^dt-fifa-Tat) .-(8) 

S c =1 + — (6.32 + 2.01exp(- 8.322 xlOV)} ...(8') 

7.6xl0 7 </. 1 FV K J 

Vi-iplpM ...<8") 

wherein g is gravitational acceleration, 7? is the 
coefficient of viscosity of the air 4 , S c is a sliding 
correction coefficient , and/? is a speed coefficient. Here, 
while particles do not conform to the Stokes ' law when they 
become small, the sliding correction coefficient (S c ) , which 
is a coefficient for correcting this phenomenon, is 
represented by the above expression (8*) according to the 
diameter (di) of particles determined by expression (7). 

On the other hand, the falling speed (Vi) of the 
particles after the complete evaporation of particle- 
constituting solvent component is represented by the above 
expression (8") using the specific gravity ( p ) of the 
compound, specific gravity ( p d ) of the particle- 
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constituting solvent, speed coefficient ( j8 ) , and diameter 
(di) of particles determined by expression (7*). 

From transference speed (k 4 ) of the compound in the 
air 4, which will be explained later, floating height of 
5 particles (Hi), size of the room (V 4 = L 4 x w 4 * height), and 

the falling speed (v £ ) of particles determined by expression 
( 8 ) or ( 8 " ) , the number of floating large and small particles 
2 and 3 (n t , i = 2, 3) is expressed as: 

dn i n i n i I /c\\ 

— = v i — ...(9) 

dt H; L 4 KJ 

10 Here, the floating height (Hi) of particles refers to the 

maximum height of the particle distribution from the floor 
5 after the particles are reflected by the spraying site 

1 and then are suspended in the case where the chemical is 
sprayed to the spraying site 1. It is positioned higher as 

15 the particles are finer. 

Using the diffusion coefficient (D c ) of compound in 
the floor 5, wall 6, and ceiling 7, the transference speed 
(k x ) of the compound in the spraying site 1 is expressed as 
(D c /t)°\ 

20 From the transference speed ( k 4 ) of the compound in 

the air 4 indicated in the following, the transference speed 
( ki, i = 2 , 3 ) of the compound in the large and small particles 

2 and 3 is expressed as k 4 /100. 

From the ventilation rate (G) , the transference speed 
2 5 ( k 4 ) of the compound in the air 4 is expressed as GV 4 /A 4 . Here, 
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the above-mentioned A 4 is the cross-sectional area of the 
room with respect to the moving direction of the air 4 
(arrowed direction in Fig. 2). 

Using the diffusion coefficient (D c ) of compound in 
the floor 5, wall 6, and ceiling 7, the transference speed 
(kj, j = 5, 6, 7) of the compound in the floor 5, wall 6, 
and ceiling 7 is expressed as (D c /t) 05 . 

From the transference coefficient (kj) of the compound 
in the spraying site 1, transference speed (k 4 ) of the 
compound in the air 4, area (A t ) of the spraying site 1, 
falling speed (Vi) of particles determined by expression (8) 
or (8"), Fugacity capacity (Z x ) of the spraying site 1 
determined by expression (2), and Fugacity capacity (Z 4 ) of 
the air 4 determined by expression (4), the transference 
coefficient (D 14 ) of the compound between the spraying site 
1 and the air 4 is represented as follows: 

D l4 « ...(10) 

l/(* l AA^i) + l/(MiAZ4) 

From the transference speed (kj of the compound in 

the large and small particles 2 and 3, transference speed 

( k 4 ) of the compound in the air 4, diameter (di) of particles 

determined by expression (7) or (7*), Fugacity capacity (Zi) 

of particles determined by expression (3), and Fugacity 

capacity (Z 4 ) of the air 4 determined by expression (4), the 

transference coefficient (D i4 , i = 2, 3) between the large 

and small particles 2 and 3 and the air 4 is represented 
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as follows : 

Z> i4 7 r ...(11) 

wherein A i (= 7Cd*) is the surface area of particles. 

From the transference speed (k 4 ) of the compound in 
5 the air 4, transference speed (k^) of the compound in the 

floor 5, wall 6, and ceiling 7, surface area (Aj) of the floor 
5, wall 6, and ceiling 7, Fugacity capacity (Z 4 ) of the air 
4 determined by expression (4), and Fugacity capacity (Zj) 
of the floor 5 , wall 6 , and ceiling 7 determined by expression 
10 (5), the transference coefficient (D 4j/ j = 5, 6, 7) of the 

compound between the air 4 and the floor 5, wall 6, and 
ceiling 7 becomes: 

D 4i ...(12) 

4 ' l/(M/ 4 ) + l/(Wi) 

Here, the surface area (A 3 ) of each of the floor 5, wall 6, 

15 and ceiling 7 is determined by the size of the room (V 4 = 

L 4 x W 4 x height) and the area (A t ) of the spraying site 1. 

The volume (V x ) of the spraying site 1 is determined 

by the diffusion coefficient (D c ) in the floor 5, wall 6, 

and ceiling 7. 

20 Using the diameter (d i ) of suspended particles 

determined by expression (7) or (7') and the evaporation 

constant (a) determined by expression (1) , the volume 
change (dVi/dt) of the large and small particles 2 and 3 is 
represented as: 



49 



SCC 97-06 



dV. re „ v 

— l - ad t ...(13) 

dt 2 1 V } 

Namely, as the oil component or moisture in the large and 
small particles 2 and 3 evaporates over time, their volume 
(Vi, i = 2, 3) decreases. 

Assuming that each of the floor 5, wall 6, and ceiling 
7 before residually spraying the chemical is like a thin 
film and that, as the chemical infiltrates into the film, 
its thickness increases so as to enhance the volume thereof; 
from the surface area ( A^ ) of the floor 5, wall 6, and ceiling 
7 and diffusion coefficient (D c ) of the compound in the floor 
5, wall 6, and ceiling 7, the volume (Vj, j = 5, 6, 7 ) of 
floor 5, wall 6, and ceiling 7 is indicated as follows: 
V^JDJLA) ...(14) 

(ii) Fugacity Calculation (S3) 

The behavior of the compound in the spraying site 1 
is expressed in the form of differential equation concerning 
Fugacity (f x ) as: 

^V.Z, = -^DJtA t ZJ v V -change 

dt 

3 

F- change ~*~ Yl ^ " Deposition^) ...(15) 

- D lA ( f v - f 4 ) T ransferenceQ) 

- K i V i Z l f l Degradation 

Here, the terms of V-change, Deposition (i), Transference 
(1), and Degradation respectively indicate volume change 
(increase over time) of the spraying site 1, attachment 
accompanying the falling of the large and small particles 
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2 and 3, amount of transference of the compound between the 
spraying site 1 and the air 4, and change in amount of 
photodegradation of the compound. 

The behavior of the compound in the large and small 
5 particles 2 and 3 is expressed in the form of differential 

equation concerning Fugacity (f if i = 2, 3) as: 

&V X Z X = V- change 

Q F -change ~ A^/l ~ A) Transference^) 

- D i4 (fi~ f 4 ) Transference^) 

- K i V i Z i f t Degradation 

ry Here, the terms of V-change, Transference (1), Transference 

if}' ( 4 )/ and Degradation respectively indicate volume change 

g 10 (decrease over time) of the large and small particles 2 and 

H 3, amount of transference of the compound between the 

spraying site 1 and the air 4, amount of transference of 
the compound between the large and small particles 2 and 

3 and the air 4, and change in amount of photodegradation 
15 of the compound. 

The behavior of the compound in the air 4 is expressed 
in the form of differential equation concerning Fugacity 
(f 4 ) as: 

V 4 Z 4 « -GV A Z 4 f 4 Ventilation 



dt 

F - change - D l4 (f 4 - f ) Transference^) 



3 
7 



n i D M-fi) Transference^) ...(17) 
A,, (A - fj ) Transfer ence( j) 



- K 4 V 4 Z 4 f 4 Degradation 
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Here, the terms of Ventilation, Transference (1), 
Transference (i), Transference (j), and Degradation 
respectively indicate amount of discharge of the compound 
outdoors, amount of transference of the compound between 
5 the air 4 and the spraying site 1, amount of transference 

of the compound between the air 4 and the large and small 
particles 2 and 3, amount of transference between the air 
4 and the floor 5, wall 6, and ceiling 7, and change in amount 
of photodegradation of the compound, 
10 The behavior of the compound in the floor 5, wall 6, 

and ceiling 7 is expressed in the form of differential 
equation concerning Fugacity (f^, j = 5, 6, 7 ) as follows: 

^j-VjZj = -jDJtAJjfj V- change 

3 

F - change + y . . ^ nV-Z. f Deposition (/) . . . (i 8) 

- D 4j ( / ; - f 4 ) Transference^) 

- KjVjZj fj Degradation 

Here, the terms of V-change, Deposition (i), Transference 
15 (4), and Degradation respectively indicate volume change 

(increase over time) of the floor 5, wall 6, and ceiling 
7 , attachment accompanying the falling of the large and small 
particles 2 and 3, amount of transference of the compound 
between the floor 5, wall 6, ceiling 7 and the air 4, and 
20 change in amount of photodegradation of the compound. 

The above-mentioned seven kinds of differential 
equations (15) to (18) are simultaneously formed and are 
solved by Runge-Kutta-Gill method, so as to compute Fugacity 
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(fi to f 7 ). 

When solving these simultaneous differential 
equations, it is necessary to set an estimation nick time 
width (dt) which is a minute time unit. Namely, the 
5 estimation nick time width is used such that solutions of 

the simultaneous differential equations are initially 
determined at a time (t 0 ) / and then solutions of the 
simultaneous differential equations are determined at a 
time (t 0 + dt) to which the estimation nick time width is 

10 added. As solutions are obtained while estimation nick time 

widths are successively added, temporally changing Fugacity 
can be determined. Theoretically, as the set time of the 
estimation nick time width is shorter, more accurate 
solutions can be obtained, though necessitating a very long 

15 calculation time. By contrast , when the set time is too long, 

solutions tend to diverge, thereby generating errors. 

Therefore, in the present invention, the estimation 
nick time width is set shorter when a very large change occurs 
in a chemical, whereas it is set longer when there is no 

20 large change. 

Specifically, mass balance is always confirmed such 
that the amount of input of the chemical and the resulting 
solution coincide with each other, and the estimation nick 
time width is set longer when the mass balance does not 

25 fluctuate greatly, whereas it is set shorter when the mass 

balance starts fluctuating. For example, when the 
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fluctuation of mass balance is set to an accuracy of ±5%, 
the estimation nick time width is always set so as to 
constantly satisfy the relationship of: 

compound input amount/ ( existing amount + degrading 
amount + discharging amount) = 0.95 to 1.05 
namely, such that the fluctuation of mass balance lies within 
the range of ±5%. 

Here, the above-mentioned compound input amount is the 
compound emission amount (A a ). Since the temporal amounts 
of compound in the seven kinds of media are determined by 
the simultaneous differential equations, they are summed 
up so as to compute the existing amount as shown in the 
following. Also, the degrading amount and discharging 
amount are as follows : 



3 7 

existing amount = V Y Z l f l + ^ w./.V)Z. + ^ fk^k^k 
degrading amount - + K.nJV.Z^^ KJ k V k Z k 

3 

discharging amount =G V 4 Z 4 f 4 + ^ G nV i Z i f i 

i-2 



.(19) 



With reference to the flowchart of Fig. 5, the method 
of setting the estimation nick time width (dt) will be 
explained . 

First, an initial value of estimation nick time width 
(dt ) is inputted (S21 ) . Then, an upper limit set value (e.g. , 
0.1[%]) which is the upper limit of difference in mass 
balance, and a lower limit set value (e.g., 10" 6 [%]) which 
is the lower limit of difference in mass balance are inputted 
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(S22) . Thereafter, Fugacity and mass balance at t = t 0 are 
calculated (S23 and S24), and Fugacity and mass balance at 
t = t + dt (or t 0 + dt) are calculated (S25 and S26). 

It is judged whether the fluctuation in mass balance 
is within the range of ±5% or not (S27). If the fluctuation 
in mass balance is within the range of ±5%, it is judged 
whether the difference between the mass balance at t = t 

(or t 0 ) and the mass balance at t = t + dt (or t 0 + dt) is 
at least the upper limit set value or not (S28). If it is 
judged to be at least the upper limit set value at S28, then 
solutions become more accurate when the estimation nick time 
width (dt) is made shorter. In this case, the estimation 
nick time width (dt) is multiplied by 1/2 so as to change 
its setting (S29). When the difference is judged to be 
smaller than the upper limit set value at S28, it is judged 
at S30 whether the difference between the mass balance at 
t = t (or t 0 ) and the mass balance at t = t + dt (or t 0 + 
dt) is at most the lower limit set value or not. 

When the mass balance difference is not greater than 
the lower limit set value at S30, since solutions are not 
influenced by longer estimation nick time width (dt), the 
estimation nick time width (dt) is doubled so as to change 
its setting (S31). Subsequently, it is judged whether the 
estimation nick time width (dt) changed at S31 is at most 
a maximum value (e.g., 0.1 [hour]) of estimation nick time 
width ( dt ) or not ( S32 ) . When the estimation nick time width 
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(dt) is not greater than the maximum value at S32, since 
solutions do not diverge, the estimation nick time width 
(dt) set at S31 is used. When the estimation nick time width 
(dt) is greater than the maximum value at S32 , since 
solutions may diverge, the estimation nick time width (dt) 
is reset to the maximum value (S33) . When the mass balance 
difference is greater than the lower limit set value at the 
above-mentioned S3 0, namely, when it lies between the lower 
limit set value and upper limit set value, the calculations 
are continued without changing the estimation nick time 
width (dt). 

After the step of S29, S30, S32, or S33, the operation 
returns to S25 so as to effect calculation again, and this 
process is repeated till the aimed time is attained. 

When the fluctuation in mass balance exceeds the range 
of ±5% at S27 , on the other hand, the mass balance fluctuation 
is so much that calculation is preferably effected with an 
estimation nick time width (dt) shorter than that in the 
case where the fluctuation is within the range of ±5%. 
Accordingly, the calculation is stopped once (S34), the 
lower limit set value is reset to a lower level (S35), and 
then the operation returns to the step of S23. 

Thus, when the estimation nick time width is variably 
set without being held constant, while monitoring the mass 
balance fluctuation, Fugacity can be computed accurately 
and efficiently. 
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(iii) Computation of Temporal Concentration in Indoor Air 
and Residual Amount (S4) 

The temporal concentration of the compound in the 
indoor air is computed when Fugacity (f 4 ) of the air, 
determined by the above-mentioned item (ii), multiplied by 
Fugacity capacity (Z 4 ), and Fugacity (f 3 ) of small particles 
multiplied by Fugacity capacity (Z 3 ) are summed up. Here, 
the large particles may be inhaled by a human body depending 
on the kind of chemical. In such a case, calculation is 
effected by use of Fugacity (f 2 ) and (f 3 ) of large and small 
particles . 

The residual amount of the compound is computed when 
Fugacity (f x ) of the spraying site 1 multiplied by Fugacity 
capacity (ZJ, and Fugacity (f 5 ) of the floor 5 multiplied 
by Fugacity capacity (Z 5 ) are summed up. Here, in the case 
where there is substantially no possibility of the spraying 
site 1 coming into contact with the skin, Fugacity (f 5 ) of 
the floor 5 alone may be used for calculation. 

(iv) Calculation of Estimated Exposure Amount in Inhalation 
and Inhalation Safety Coefficient (S6 and S7) 

The above-mentioned temporal concentration in the 
indoor air forms a curve shown in Fig. 3A, for example. This 
concentration curve is integrated, an accumulated 
concentration of the compound during a specific period (t x 
to t 2 ) is determined (see Fig. 3B) , and the mean concentration 
in the indoor air is computed from thus determined value. 
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While an arbitrary period is set as the specific period 
depending on the object, an appropriate period is usually 
set in view of the method of use of the product and the test 
period of toxicity data. 

Then, according to the above-mentioned mean 
concentration in the indoor air, amount of respiration, and 
exposure time, the estimated exposure amount in inhalation 
is determined. Namely, calculation of: 

estimated exposure amount in inhalation [mg/kg/day] 
= mean concentration in indoor air [mg/m 3 ] x amount of 
respiration [mVkg/min] x exposure time [min/day] 
is effected. Here, as the above-mentioned amount of 
respiration, a published value or actually measured value 
may be used. Also, when amounts of respiration are 
respectively set for adult and child, more appropriate 
estimated exposure amounts in inhalation can be obtained. 
In the case where the inhaled harmful material is not totally 
absorbed into the body but is partially discharged by 
respiration, a more appropriate estimated exposure amount 
in inhalation can be obtained when the inhalation ratio is 
taken into account. 

The inhalation safety coefficient is computed from a 
non-influential amount concerning inhalation toxicity 
examined by an animal experiment beforehand and the 
estimated exposure amount in inhalation determined above. 
Namely, it is expressed as: 
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inhalation safety coefficient = 

inhalation non-influential amount [mg/kg/day ] / 
estimated exposure amount in inhalation [mg/kg/day] 
(v) Calculation of Estimated Percutaneous Exposure Amount 
and Percutaneous Safety Coefficient (S9 and Sll) 

The above-mentioned residual amount forms a curve such 
as that shown in Fig- 4A, for example. This residual amount 
curve is integrated, the accumulated residual amount of the 
compound during a specific period (t x to t 2 ) is determined 
(see Fig. 4B) , and the mean residual amount is computed from 
thus determined value. While an arbitrary period is set as 
the specific period depending on the object, an appropriate 
period is usually set in view of the method of use of the 
product and the test period of toxicity data. 

Then, according to the mean residual amount, skin 
attachment ratio, contact area, and body weight, the 
estimated percutaneous exposure amount is determined. 
Namely, calculation of: 

estimated percutaneous exposure amount [mg/kg/day] = 
(mean residual amount [mg/m 2 ] x skin attachment ratio [%] 
x contact area [m 2 /day ] ) /body weight [kg] 

is effected. Here, as the contact area, a published value 
(e.g., 4 [mVday]) may be used. The skin attachment ratio 
is a ratio of the compound attaching to the skin when the 
latter is in contact with the floor 5 where the compound 
exists. As this value, a published value or a value 
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experimentally obtained from a model may be used. 

A model experiment method for the skin attachment 
ratio is as follows. A weight ( 8 cm x 8 cm * 8 cm; 4.2 kg) 
is placed on a denim cloth (8 cm * 10 cm) with a pressure 
similar to that of an infant in contact with a floor, and 
the denim cloth is pulled on the floor at a speed (120 cm/15 
sec) similar to the moving speed of the infant. The denim 
and floor are analyzed so as to compute the compound 
contained in the denim and floor. From the ratio 
therebetween, the skin attachment ratio is obtained. It has 
been confirmed that the skin attachment ratio obtained by 
this method is identical to or slightly higher than that 
determined from analyzed values of a hand and a floor when 
the hand is actually pressed against the floor, thereby 
proving this model experiment method to be useful for 
evaluating exposure of inhabitants. 

The percutaneous safety coefficient is computed from 
the non-influential amount concerning percutaneous 
toxicity examined by an animal experiment beforehand and 
the estimated percutaneous exposure amount determined above. 
Namely, it is expressed as: 

percutaneous safety coefficient = 

percutaneous non-influential amount [mg/kg/day]/ 
estimated percutaneous exposure amount [mg/kg/day] 

Nevertheless, in general, percutaneous non- 
influential amount has not often been determined, and there 
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are not many published values . Accordingly, a more accurate 
value can be determined from the estimated percutaneous 
exposure amount, and oral non-influential amount and 
percutaneous absorption ratio for which many published 
values exist, according to the following expression: 

percutaneous safety coefficient = oral non- 
influential amount [mg/kg/day ] / ( estimated percutaneous 
exposure amount [mg/kg/day] * percutaneous absorption ratio 
[%]) 

Here, when the percutaneous absorption ratio is unknown, 
employed is a national guideline (e.g., 10%) which usually 
exists • 

(vi) Calculation of Estimated Oral Exposure Amount and Oral 
Safety Coefficient (S13 and S14) 

From the estimated percutaneous exposure amount 
obtained in the above-mentioned item (v) , hand surface area 
ratio, and oral transference ratio, the estimated oral 
exposure amount from hand to mouth is determined. Namely, 
calculation of: 

estimated oral exposure amount [mg/kg/day] = 
estimated percutaneous exposure amount [mg/kg/day] x hand 
surface area ratio [%] * oral transference ratio [%] 
is effected. Here, the hand surface area ratio is expressed 
by (hand surface area/body surface area), for which a 
published value (e.g., 5[%]) may be used. The oral 
transference ratio is a hypothetical value, which is set 
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to 100%, for example. 

In the case where oral exposure might occur via 
tableware or food contaminated with the residually sprayed 
compound, it is required that the estimated oral exposure 
5 amount from tableware or food to mouth be added to the 

estimated oral exposure amount from hand to mouth to yield 
the total estimated oral exposure amount. For example, the 
estimated oral exposure amount from tableware is obtained 
• S when, according to the tableware residual amount indicating 

H 10 the amount of the harmful material remaining in tableware, 

jjssSs 

543 tableware use area which is the sum of tableware surface 

t y 

^ areas, and oral transference ratio from tableware, 

^ calculation of: 

,~ estimated oral exposure amount [mg/kg/day] = 

Uz 15 tableware residual amount [mg/m 2 ] x tableware use area 

[mVday] x oral transference ratio [%]/body weight [kg] 
is effected. Here, the tableware residual amount is 
expressed by (mean floor residual amount x tableware 
contamination ratio) . As the tableware contamination ratio, 
2 0 an actually measured value (e.g. , 9% ) or a hypothetical value 

may be used. 

The oral safety coefficient is computed from the 
non-influential amount concerning oral toxicity examined 
by an animal experiment beforehand and the estimated oral 
2 5 exposure amount determined above. Namely, it is expressed 

as : 
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oral safety coefficient = 

oral non-influential amount [mg/kg/day ] / estimated 
oral exposure amount [mg/kg/day] 

As mentioned in the foregoing, in the method of 
5 estimating an indoor behavior of a pesticidal compound in 

this embodiment, the environment is initially divided into 
the spraying site 1, large and small particles 2 and 3, whole 
indoor air 4 , floor 5 , wall 6 , and ceiling 7 ; and differential 
equations concerning Fugacity of the compound in them are 
10 simultaneously formed and are solved, so as to estimate the 

indoor behavior of the compound when a chemical containing 
S the compound is residually sprayed. 

=i Since the environment is thus considered to be seven 

;»j kinds of media, and exchanges between the individual media 

15 and the like are taken into account, simulation results close 

to the actual behavior of the compound can be obtained. 

The estimation nick time width set when solving the 
simultaneous differential equations is variably set, while 
constantly confirming mass balance of the compound indoors 
20 after the residual spraying, so that the amount of input 

of the chemical indoors and the resulting solution coincide 
with each other. 

Accordingly, since mass balance of the compound after 
the residual spraying is always monitored such that the 
25 amount of input of the compound indoors and the resulting 

solution coincide with each other, thereby varying the 
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estimation nick time width; the estimation nick time width 
is set longer when the mass balance fluctuates a little, 
whereas it is set shorter when the mass balance starts 
fluctuating greatly. Namely, when solving simultaneous 
differential equations including a parameter accompanying 
temporal change, the estimation nick time width is 
automatically set in response to the fluctuation in mass 
balance. Consequently, when processed by a computer, an 
accurate solution can be obtained in a short time. 

The method of evaluating safety of a pesticidal 
compound in accordance with this embodiment uses the 
estimated result mentioned above to evaluate the safety of 
the compound with respect to the human body when the chemical 
is residually sprayed. 

Accordingly, the safety of the compound with respect 
to the human body can be evaluated accurately in a short 
time. As a consequence, when formulating a chemical such 
as insecticide containing the compound, simulation can be 
easily repeated while changing conditions, thereby making 
it easier to formulate a chemical having a high safety 
conforming to the aimed object. 

Though the kind of the floor 5 is assumed to be flooring 
in this embodiment, differential equations similar to those 
mentioned above may also be formed in the case of tatami. 
In this case, however, when determining fugacity capacity 
(Z 5 ) of the floor 5 in equation (5), it is necessary for the 
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particle-constituting solvent component content (see Table 
3 ) to be changed from that of the flooring to that of tatami. 
Also, when determining the degradation constant (K 5 ) of the 
floor 5 in equation (6), the decomposition half -life is 
required to be changed from that of the flooring to that 
of tatami. 

Though Fugacity is determined by use of Runge- 
Kutta-Gill method in this embodiment, other methods may be 
used for solving differential equations, Runge-Kutta-Gill 
method, however, is preferably used since a program for the 
above-mentioned differential equations can be easily made 
by Basic, Also in the case where differential equations are 
solved by a method other than Runge-Kutta-Gill method, 
similar effects can be obtained when the estimation nick 
time width is set as mentioned above. 

In the following, examples (Examples 1 to 5 ) of the 
above-mentioned method of estimating an indoor behavior of 
a pesticidal compound and the safety evaluation method using 
the same will be explained with reference to the accompanying 
drawings . 

The evaluation of the method of estimating an indoor 
behavior of a pesticidal compound in these examples is 
effected as the results obtained by estimation and the 
actually measured values obtained by measurement of the 
actual indoor behavior of the compound are compared with 
each other. 
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First, the measurement of the actual indoor behavior 
of the compound will be explained. 

The employed aerosol can (300 [ml]) contains 0.9 [g] 
of d-phenothrin (trademark: Smithrin; C 23 H 26 0 3 : (3- 
phenoxyphenyl )methyl ( 1R) -cis-trans-2 , 2-dimethyl-3- ( 2- 
methyl-l-prophenyl ) cyclopropane carboxylate) and 1.1 [g] 
of d-tetramethrin (trademark: Neo-Pyramin Forte; C 19 H 25 0 4 : 
(1,3,4,5,6, 7-hexahydro-l , 3-dioxo-2H-isoindol-2- 
yl )methyl ( 1R) -cis-trans-2 , 2-dimethyl-3- ( 2-methyl-l- 
prophenyl ) cyclopropane carboxylate) . 

The indoor environment is assumed to be a six-mat room 
(9.72 [m 2 ]) of a typical apartment in Japan. In conformity 
to the description of the above-mentioned aerosol can, 
spraying was effected at a rate of 60 [sec/m 2 ] from a distance 
of about 20 [cm] from the floor and wall. The spraying site 
had a width of 10 [cm] or 15 [cm] from the boundary between 
the floor and the wall toward the floor, and a width of 10 
[cm] or 15 [cm] from the boundary toward the wall [see Fig. 
2 ] . The total spraying period was 2 . 5 minutes for 2.52 [m 2 ] 
in the region having a width of 10 [cm], and 3.7 minutes 
for 3.78 [m 2 ] in the region having a width of 15 [cm]. 

Five kinds of measurement were effected while the room 
conditions were changed as shown in Table 7. Namely, 
measurement in the case with continuous use (Example 1), 
case with changing ventilation rate (Example 2), case with 
an open window (Example 3), and case with air-conditioner 
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on (Example 4), and measurement for skin attachment ratio 
(Example 5) were performed. 




CQ 




EXAMPLE 5 


FLOOR, HAND, 
DENIM CLOTH 


CN 


OFF 


tO CN 
CN VO 


CLOSED 


w 
Z 


1.54 


EXAMPLE 4 


AIR 


<o 

CN 


ON 


to o 

CN 


CLOSED 


w 


1.56 


EXAMPLE 3 


AIR 


to 

<N 


OFF 


o 

CN VO 


5min - OPEN 
OR 2h- OPEN 


BASICALLY NONE 


1.86, 10.6 


EXAMPLE 2 


AIR 




OFF 


to o 

CN 


CLOSED 


w 
Z 


0.50 OR 1.58 
OR 4. 14 


EXAMPLE 1 


AIR, FLOOR, WALL, 
CEILING 


X 
J* 

<L> 

B 

«o 

CN 


OFF 


<ri o 

CN 


CLOSED 


SYNTHETIC LIGHT 
AND NATURAL LIGHT 
1320 


OO 

»o 




MEASUREMENT SITE 


SPRAYING PERIOD [min/day] 


AIR - CONDITIONER STATE 


ROOM TEMP [°C] 
HUMIDITY [%RH] 


WINDOW STATE 


ILLUMINATION CONDITION 
ILLUMINANCE (lux) 


VENTILATION RATE [H] 



68 



SCC 97-06 

In Example 1, spraying was effected four times at a 
rate of once in two weeks. The air-conditioner was not 
actuated, and the windows of the room were closed. The room 
temperature and humidity were controlled so as to become 
25 [°C] and60[%], respectively. The illumination condition 
was such that synthetic light illuminates indoors from 7 
a.m. to 10 a.m. and from 6 p.m. to 11 p.m., while natural 
light was provided through two windows. The ventilation 
rate was set to 1.5 [h** 1 ], which was minimally necessary for 
2 10 a Japanese woman and an infant to inhabit. The temperature, 

humidity, illuminance, and ventilation rate of the room 
S actually measured under such a condition are shown in Table 

: U 7, each being represented by the mean value in measurement. 

- The temperature, humidity, and illuminance are measured at 

1: 15 the center of the room. The ventilation rate was confirmed 

u 

™ before and after the measurement. 

In Example 2, the illumination condition was in the 
dark as natural light was cut. For performing measurement, 
the ventilation rate was changed to 0 . 5 to 4 [h _1 ] since the 
20 actual ventilation rate of Japanese apartment with closed 

windows is 0 . 5 to 3 [h -1 ]. 

In Example 3, immediately after spraying, a window was 
opened for 5 minutes or 2 hours in conformity to the 
description of the aerosol can. While the illumination 
25 condition was basically in the dark (illuminance 0 [lux] 

in Table 7) , natural light entered when the window was open 
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(mean illuminance 1890 [lux] in Table 7). 

The sampling processing in each of the above-mentioned 
measurements was effected as follows. 

In Examples 1 to 4 , the indoor air was sampled after 
a specified period from the spraying. Here, as for d- 
phenothrin and d-tetramethrin , at the positions 
corresponding to the respiration regions of an infant and 
an adult, air suction amount was adjusted such that particles 
having a diameter of 30 [//m] or less can be collected, 
whereby not only the compounds in a gas state but also those 
in a suspended particle state were sampled. 

In Example 1, at the centers of floor, wall, and 
ceiling, strips of materials identical thereto each having 
a size of 5 [cm] x 5 [cm] were placed before spraying, whereby 
sampling was effected at the floor, wall, and ceiling. As 
for the wall, sampling was effected at three places 
consisting of the inflow and discharge sides of the air and 
a lateral side apart therefrom. 

In the collecting tube that had sampled the air, 
flooring strip, wallpaper strip, ceiling strip, hand, and 
denim cloth, d-phenothrin and d-tetramethrin were analyzed 
by the methods shown in Table 8 . In these methods , an organic 
solvent containing BHT ( 2 , 6-di-t-butyl-p-cresol ) was used 
to extract these compounds, and cleanup was effected when 
necessary, so that analysis was effected by mass 
spectrometer (MSD) attached to a gas chromatography (GC). 
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In each process, in order to prevent the compound from 
vaporizing due to evaporation of the solvent, diethylene 
glycol or octanol is added. Both d-phenothrin (mass: 183 
[m/z]) and d-tetramethrin (mass: 164 [m/z]) were detected 
in 9 minutes under the GC condition noted in the following. 
As the apparatus, Shimadzu GC-MS QP-1100EX (Mode EI, 7 0 eV) 
was used with a DC-608 column (inner diameter: 0.53 [mm], 
length: 30 [m], thickness: 0.8 [ jum] ) programmed to raise 
temperature at a rate of 10[°C/min] from 205[°C], whereas 
the temperature of injection and ion source was 260[°C]. 
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The results of concentrations of d-phenothrin and 
d-tetramethrin in the air and the residual amounts in the 
floor, wall, and ceiling measured by the above-mentioned 
methods will be explained. 

In Example 1 , the concentration in the air was measured 
for 125 days including the time in which spraying was 
effected once (2.5 minutes ) per 2 weeks for 8 weeks . As shown 
in Fig. 7, the concentration in the air exhibited its maximum 
level immediately after spraying, where the average of 
d-phenothrin at the heights of 25 [cm] and 120 [cm] from 
the center of the floor was 752 [^g/m 3 ] and that of d- 
tetramethrin was 1040 [/^g/m 3 ], and rapidly decreased 
thereafter. Here, since the half-life of each of both 
compounds was short, i.e., about 20 minutes, there were 4 
peaks in the 8 weeks of spraying period. The concentration 
in the air after the spraying period was less than 0.10 [ ju 
g/m 3 ] in both of them. Further, there was no substantial 
difference in concentration between the heights of 25 [cm] 
and 120 [cm] from the floor. Also, no increase was seen in 
the concentration in the air due to the increase in the number 
of spraying operations. In Fig. 7, concentrations which 
were less than the detection limit (0.05 [/^g/m 3 ]) and thus 
could not be measured are plotted at the detection limit. 

The total spraying amount was 250 [g], including 1.11 
[g] of d-phenothrin and 1.38 [g] of d-tetramethrin, which 
equal 1 . 2 aerosol cans of 300 [ml ] . Since 5 million aerosol 
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cans of 300 [ml] for residual spraying are annually sold 
in Japan, annual consumption per household is calculated 
as 0.13 can. Accordingly, it is seen that the spraying 
amount used for measurement (1.3 cans per room) is 
considerably greater than the actual annual consumption. 

As shown in Table 9, the mean concentration in the air 
during the spraying period and the mean annual concentration 
in the air including 8-week spraying period were calculated 
as 2.35 [>^g/m 3 ] /0.43 [/ig/m 3 ] for d-phenothrin, and 3.25 [ ju 
g/m 3 ]/0.53 [//g/m 3 ] f or d-tetramethrin . Consequently, it is 
seen that the amount of existence of d-phenothrin with 
respect to d-tetramethrin in the air corresponds to their 
ratio of existence within the aerosol can. In Table 9, for 
those less than the detection limit value/ the mean 
concentration in the air was calculated by use of the 
detection limit value. 
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In Example 2, concentrations of d-phenothrin and 
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:. _ i 

15 



20 



25 



d-tetramethrin in the air after 3 . 7 minutes from the spraying 
with changed ventilation ratio were measured . Fig . 8A shows 
only the results of d-phenothrin. This has indicated that 
the half -lives in the air at ventilation ratios of 1,58 [ h" 1 ] , 
0.50 [h _1 ], and 4.14 [h" 1 ] are 21 minutes, 58 minutes, and 
11 minutes, respectively. 

In Example 3, immediately after the chemical was 
sprayed for 2,5 minutes, the window was opened, and the 
concentrations of d-phenothrin and d-tetramethrin in the 
air were measured. Fig. 8B shows the results of d- 
phenothrin. This has indicated that the half -lives in the 
cases where the window was opened for 5 minutes, the window 
was opened for 2 hours, and the window was kept closed are 
14 minutes, 4 minutes, and 20 minutes, respectively. Also, 
the mean concentration (2.88 [//g/m 3 ] when opened for 5 
minutes, 1.19 [//g/m 3 ] when opened for 2 hours) was less than 
1/5 that of the closed state (14.8 [//g/m 3 ]). 

The measured values of Figs. 8A and 8B are mean values 
of the concentrations in the air respectively sampled at 
heights of 25 [cm] and 120 [cm] from the floor. Also, for 
those less than the detection limit and thereby incapable 
of measurement were plotted as the detection limit value 

(0.05 [>Ug/m 3 ]). 

Table 10 lists the above-mentioned half-lives. These 
results have indicated that, in the residual spraying, the 
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ventilation of room, opening the window immediately after 
spraying in particular, is highly influential in the 
concentration in the air. For reference, Table 10 shows 
half -lives of compounds (d-resmethrin and d-tetramethrin) 
in the air when a chemical is sprayed in the space. 
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Meanwhile, in the state where the window is open as 
with Example 3, for example, the half-lives caused by 
photodegradation are estimated according to the above- 
mentioned Atmospheric Oxidation program from the structures 
of compounds to be 2 . 1 hours for d-phenothrin and 1.7 hours 
for d-tetramethrin . Since the half -lives in the air are 
shorter than the half-lives caused by photodegradation, the 
decreases in concentration of d-phenothrin and d- 
tetramethrin are considered to be mainly caused by 
ventilation rather than photodegradation. 

In Example 4, the mean concentration when the air- 
conditioner was used (11.6 [/Zg/m 3 ] for d-phenothrin, 15.5 
[//g/m 3 ] for d-tetramethrin ) was substantially the same as 
the mean concentration when the air-conditioner was not used 
(14.8 [/ig/m 3 ] for d-phenothrin, 20.2 [//g/m 3 ] for d- 
tetramethrin) . This has indicated that the influence of the 
air convection and absorption of the air into the filter 
of the air conditioner upon the concentration in the air 
is substantially negligible. 

Figs. 9A and 9B and the above-mentioned Table 9 show 
the results of measurement of residual amounts at the floor, 
wall, and ceiling for 125 days in Example 1. The values of 
the residual amounts at the wall in Fig. 9B and Table 9 are 
means of values on the air inflow and discharge sides and 
their lateral side, whereas the residual amounts of wall 
and ceiling in Fig. 9B are plotted as the mean values of 
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residual amounts at the wall and ceiling. 

As the air flows in and out, the residual amount at 
the wall on the inflow side (28 [/zg/m 2 ] for d-phenothrin, 
35 [//g/m 2 ] for d-tetramethrin) was less than the residual 
amount at the wall on the discharge side or lateral side 
(45 [//g/m 2 ] and 58 [jug/m 2 ] on the discharge side, 33 [ JUL 
g/m 2 ] and 43 [//g/m 2 ] on the lateral side). Also, while the 
residual amounts of d-phenothrin and d-tetramethrin 
enhanced as the number of spraying operations increased, 
they gradually decreased after spraying. Here, half-lives 
of d-phenothrin and d-tetramethrin at the floor were 
respectively 13 days and 20 days, whereas their half-lives 
at the wall/ceiling were respectively 31 to 41 days and 24 
to 75 days. These values are longer than those in the air, 
since the materials of floor, wall, and ceiling contain a 
large amount of antioxidants. Also, the residual amount at 
the floor, except for the site sprayed with the chemical, 
is greater than the residual amount at the wall or ceiling, 
since the suspended particles attach to the floor. 

In Example 5, a hand attachment test in which a hand 
is directly pressed against the floor where the compound 
remains, and a wiping test using a denim cloth instead of 
the hand are performed, so as to measure the skin attachment 
ratios of d-phenothrin and d-tetramethrin at the flooring. 
When Student's t-test (P ^ 0.05) was performed, the skin 
attachment ratios obtained by the wiping test (0.31[%] and 
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0.33 [%] respectively for d-phenothrin and d-tetramethrin 
after 6 hours from spraying, 0.17[%] and 0.09[%] after 24 
hours from spraying) were equal to or higher than the results 
obtained by the hand attachment test (0.22[%] and 0.32[%] 
respectively for d-phenothrin and d-tetramethrin after 6 
hours from spraying, 0.06[%] and0.07[%] after 24 hours from 
spraying). This has indicated that the safety of 
percutaneous exposure can be sufficiently estimated when 
the skin attachment ratio is determined by the wiping test 
alone. 

In the following, simulation results obtained 
according to the method of estimating an indoor behavior 
of a compound in accordance with the present invention will 
be explained- The condition here is the same as that of the 
above-mentioned Example 1. As the primary condition, the 
values described in Tables 1 to 4 of the first embodiment 
were inputted. 

When the behavior of d-phenothrin caused by residual 
spraying during the period after 2.5 minutes to 24 hours 
from the spraying was simulated, the results shown in Table 
11 were obtained. In the table, "suspension period" refers 
to a period in which the particles floating at the position 
lower than 2.4 [m] fall onto the floor or disappear in the 
air, whereas "transference amount" refers to the ratio of 
compound amount transferred within 24 hours with respect 
to the whole spraying amount. 
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According to these results, the main solvent in the 

large particles (diameter: 20 [/im]) was completely 
evaporated after 35 seconds, and disappeared from within 
the air after 2 hours. The main solvent in the small 
5 particles (diameter: 5 [//m]) was completely evaporated 

after 2.2 seconds, and existed in the air for 30 hours. The 
chemical containing the particles in the air (40[%] of the 
amount of chemical immediately after spraying ) was 
distributed to the spraying site by 1 . 0 [ % ] , floor by 7 . 0 [ % ] , 

10 wall by 0.44[%], and ceiling by 0.10[%], whereas the 

remaining 29[%] went outdoors, and 3.2[%] was degraded in 
the air. Of the chemical at the spraying site (60[%] of the 
amount of chemical immediately after spraying), 0.16[%] of 
the sprayed amount was evaporated into the air within 2 4 

15 hours, 1.7[%] was degraded by light or oxidation, and the 

remaining 58[%] was left at the spraying site. 

When the simulation results concerning the 
concentration in the air and the residual amounts at the 
spraying site and floor/wall/ceiling and the actually 

2 0 measured values which is mentioned above were compared with 

each other, relationships shown in Figs. 10A and 10B were 
obtained. In these graphs, solid lines show simulation 
results, whereas plotted points indicate actually measured 
values. These actually measured values were derived from 

25 the results of Example 1. 

From these graphs, it can be seen that the estimated 
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results according to the method of estimating an indoor 
behavior of a compound in accordance with the present 
invention are very close to the actually measured values. 
Namely, upon residual spraying, the most part of the chemical 
in the air is discharged outdoors or falls onto the floor. 
The chemical at the spraying site permeated through the floor 
or wall, and was slightly degraded or evaporated. 

In the following, the results obtained according to 
the method of evaluating safety of a compound in accordance 
with the present invention will be explained. Here, indoor 
exposure amounts of d-phenothrin and d-tetramethrin with 
respect to inhabitants including an infant when a chemical 
was residually sprayed at the boundary between the floor 
and wall for four times in eight weeks (sprayed for 2.5 
minutes for each time) were evaluated. 

The inhabitants inhale vapors or particles of d- 
phenothrin and d-tetramethrin existing in the air. Since 
humans inhale aerosol particles having a diameter not 
greater than 8 [ jura] to the lungs, sampling the particles 
having a size not greater than 30 [//m] and vapor chemicals 
is sufficient for evaluating the exposure in inhalation. 

The following Table 12 shows the estimated exposure 
amount in inhalation calculated by use of the mean indoor 
concentration in the air (see "spraying period" in Table 
9), respiration amount, and exposure time as explained in 
item (iv) in the first embodiment. Here, the mean indoor 
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respiration amount of Japanese is 0.220 [L/min/kg] for an 
adult, and 0.235 [L/min/kg] for an infant. Also, the 
exposure time was presumed to be 24 [h/day], taking account 
of the worst case where inhabitants are confined in the room 
for 24 hours after spraying. 
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COMPOUND 


SUBJECT 


INHALATION 
(WINDOW CLOSED) 


INHALATION 
(WINDOW OPEN) 


PERCUTANEOUS 


ORAL 


d-phenothrin 


INFANT 


0.795 


0.159 


2.71 


0.135 




ADULT 


0.744 


0.149 


0.553 




d-tetramethrin 


INFANT 


1.12 


0.224 


2.15 


0.107 




ADULT 


1.01 


0.202 


0.438 
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In the table, "window closed" refers to the case where 
the window is closed throughout the exposure time, whereas 
"window open" refers to the case where spraying was effected 
in the state with the closed window and then the window was 
opened for 5 minutes (Example 3). 

Also, Table 12 shows the estimated percutaneous 
exposure amount calculated by use of the mean floor residual 
amount (see "spraying period" in Table 9), skin attachment 
ratio, contact area, and body weight as explained in item 
(v) in the first embodiment. Here, it was assumed that the 
contact area for an infant and an adult was 4.0 [mVday], 
the body weight of the infant was 10,2 [kg], and the body 
weight of the adult was 50 [kg]. As the skin attachment 
ratios for d-phenothrin and d-tetramethrin at the flooring, 
values evaluated by integration and extrapolation of the 
measured values in Example 5 were used. Namely, the skin 
attachment ratio was 0.31 [%] for d-phenothrin, and 0.25[%] 
for d-tetramethrin. These values were used for calculating 
the estimated percutaneous exposure amount. 

Further, Table 12 shows the estimated oral exposure 
amount calculated by use of the estimated percutaneous 
exposure amount, hand surface area ratio, and oral 
transference ratio as explained in item (vi) in the first 
embodiment. Here, while the hand surface area ratio was 
assumed to be about 5[%], and the oral transference ratio 
was presumed to be the worst case of 100[%], only the 
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estimated oral exposure amount of the infant was calculated. 

Table 13 lists the inhalation safety coefficient, 
percutaneous safety coefficient, and oral safety 
coefficient obtained according to the estimated exposure 
amount in inhalation, estimated percutaneous exposure 
amount, and estimated oral exposure amount determined by 
the above-mentioned calculations. At the same time, the 
total safety coefficient determined from the above- 
mentioned three safety coefficients are shown in the table. 
Further shown in the table is the safety coefficient with 
respect to both of d-phenothrin and d-tetramethrin computed 
when the sum of the respective reciprocals of the safety 
coefficient (inhalation, percutaneous, or oral) of d- 
phenothrin and the safety coefficient (inhalation, 
percutaneous, or oral) of d-tetramethrin is determined and 
then is turned into the reciprocal . 
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TOTAL 
(WINDOW OPEN) 


80,400 
116,000 


12,200 
15,800 


10,600 
13,900 


TOTAL 
(WINDOW CLOSED) 


21,300 
24,400 


2,780 
3,180 


2,460 
2,820 


ORAL 


923,000 


100,000 


90,300 


PERCUTANEOUS 


369,000 
1,810,000 


465,000 
2,280,000 


206,000 
1,010,000 


INHALATION 
(WINDOW OPEN) 


116,000 
124,000 


14,400 
15,900 


12,800 
14,100 


INHALATION 
(WINDOW CLOSED) 


23,100 
24,700 


2,880 
3,190 


2,560 
2,820 


SUBJECT 
COMPOUND 


d-phenothrin 
INFANT 
ADULT 


d-tetramethrin 
INFANT 
ADULT 


BOTH 

INFANT 
ADULT 
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Here, the non-influential amount of a rat concerning 
d-phenothrin was 18 . 4 [mg/kg/day] (= non-influential in-air 
concentration 210 [mg/m 3 ] * rat respiration amount 0.3 65 
[L/min/kg] * rat exposure time 4 [h/day]; 4-week test with 
exposure of 6 days a week) in inhalation, 1000 [mg/kg/day] 
(21-day test) for percutaneous exposure, and 125 
[mg/kg/day] (6-month test) for oral exposure. On the other 
hand, the non-influential amount of a rat concerning d- 
tetramethrin was 3 . 22 [mg/kg/day] (= non-influential in-air 
concentration 49 [mg/m 3 ] * rat respiration amount 0.365 
[L/min/kg] x rat exposure time 3 [h/day]; 4-week test with 
exposure of 7 days a week) in inhalation, 1000 [mg/kg/day] 
(21-day test) for percutaneous exposure, and 10.75 
[mg/kg/day] (6-month test) for oral exposure. 

As a result, it has been confirmed that the safety 
coefficient of the infant is lower than that of the adult, 
and that the safety coefficients of d-phenothrin were higher 
than the safety coefficients of d-tetramethrin except for 
the percutaneous safety coefficient. Also, since the 
inhalation safety coefficient is the lowest of the three, 
it substantially determines the value of total safety 
coefficient when exposed by all of the three kinds of 
exposure routes . 

( 2 ) Second Embodiment 

With reference to the drawings, the second embodiment 
of the present invention will be explained in the following. 
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The method of simulating an indoor behavior of a pesticidal 
compound in accordance with this embodiment estimates the 
behavior of the compound in the case where the floor in the 
first embodiment is constituted by a rug having ears of fiber. 
For convenience of explanation, members identical to those 
shown in the drawings of the previous embodiment will be 
referred to with numerals or letters identical thereto, 
without their explanations provided. Explained here is a 
case where only the kind of floor at the spraying site is 
a carpet (e.g., the floor spraying region indicated by c 
in Fig. 2 is assumed to be a carpet) , whereas the other part 
of the floor is treated as a mere plane. Accordingly, only 
the portion different from the first embodiment will be 
explained . 

Here, though the kind of floor (whole surface) may be 
a carpet, only the suspended particles fall onto the floor 
other than the spraying site, whereby the behavior of the 
floor does not greatly differ from that in the case of 
flooring. Accordingly, it is often only necessary for 
Fugacity capacity of the floor to be corrected with the 
particle-constituting solvent component content alone. 

First, the carpet is modeled as shown in Figs. 6A, 6B, 
and 6C. Namely, the carpet is divided into a plurality of 
carpet fibers (ears of fiber) planted on a substrate, and 
space portions existing between the fibers . When a chemical 
is residually sprayed in a room in which such a carpet covers 
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the whole floor surface, the chemical infiltrates into the 
carpet fiber at the spraying site, and accumulates in the 
space portion. 

While the method of simulating an indoor behavior of 
a pesticidal compound in accordance with this embodiment 
also conforms to the flowchart of Fig. 1, at SI, as the 
primary condition, the carpet-related primary condition 
shown in Table 14 is added to the primary condition of the 
first embodiment (see Tables 1 to 4 ) . Here, the cross- 
sectional area (A c ) of the space portion of the carpet refers 
to the area in which the space portion is in contact with 
the indoor air (see Figs. 6A, 6B, and 6C). 
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The examples of values in Table 14 are based on the 
values shown in Figs • 6A, 6B, and6C. Here, one space portion 
is set to a size of 0.1 [mm] * 0.1 [mm] x 3 [mm], and one 
square of carpet fiber is set to a size of 2.5 [mm] x 2.5 
5 [mm] x 3 [mm] . Also, it is assumed that 16 [pieces ] of carpet 

fiber exist in 1 [cm 2 ]. 

In this case, the individual parameters can be 
determined as : 

ig A c = (0.1 [mm] x o.l [mm] x 16 [pieces])/l [cm] 

Q 10 = 1.6 x 10" 3 [m 2 /m 2 floor] 

ITj A cl = (0.1 [mm] x 3 [mm] x 4 [faces] x 16 [pieces] )/l 

m [cm] 

q =0.192 [m 2 /m 2 floor] 

r=% a x = 1 - a c 

J 15 = 1 - 1.6 x 10- 3 [m 2 /m 2 floor] 

A al — Aj + A cl 

= 1.19 [m 2 /m 2 floor] 
Here, A x , which is the area occupied by the spraying site 
in the first embodiment, is replaced by the area in which 
20 the carpet fibers in 1 [m 2 ] are in contact with the indoor 

air . 

As the secondary condition at S2, the carpet-related 
secondary condition shown in Table 15 is added to the 
secondary condition of the first embodiment. Here, k x and 
25 V lf which are respectively transference speed of the 

compound in the spraying site and the volume of the spraying 
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site in the first embodiment (see Table 6) , are respectively 
replaced by the transference speed of the compound in the 
carpet fiber and the volume of the whole fiber of the carpet. 
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A method of computing the definite factor conditions 
and temporal change conditions shown in the above-mentioned 
Table 15 will be explained hereinafter as exemplified by 
the case using an aqueous aerosol, i.e., case where the 
particle-constituting solvent is water. 

From the cross-sectional area (A c ) of the space portion, 
room temperature (T) , and humidity (H), the evaporation 
constant (B^) of the chemical entered into the space portion 
of the carpet is expressed as: 



g R d = A c h(H m - H)IC H 

SU .10 C H = 0.24 + 0.46// 

W h = 3.06xl0 _4 (r-7;) l/3 



..(20) 



wherein T d is the temperature of the space portion of the 
carpet (corresponding to the wet-bulb temperature), and H m 
is the humidity (corresponding to the saturated humidity 
at the wet-bulb temperature) of the space portion at T d . 

15 Here, T d and H m can be determined from the above-mentioned 

"mass-based humidity table" according to the room 
temperature (T) and humidity (H). In the case of an 
oil-based product (aerosol), while the vapor constant of 
the chemical can be computed when properties of the water 

20 and oil are compared with each other and then are multiplied 

by R d , R d = 0 may be assumed in the case where the evaporation 
of the oil is slow. The chemical entered into the space 
portion of the carpet reduces its volume upon evaporation 
and finally disappears, leaving the compound. Thus left 
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compound permeates through the carpet fiber. 

From the diffusion coefficient (D c ) in the 
floor /wall/ceiling and the area (A al ) in which the carpet 
fiber is in contact with the compound, the volume (V^ of 
the whole fiber of the carpet is determined as follows: 



The transference speed (k c ) of the compound in the 
product is expressed as k 4 /100 by use of k 4 (= GV 4 /A 4 ) which 
indicates the transference speed of the compound in the 
indoor air. 

From the above-mentioned transference speed (k c ) of 
the compound in the product, transference speed (k 4 ) of the 
compound in the indoor air, area (A t ) of the spraying site, 
cross-sectional area (A c ) of the space portion, Fugacity 
capacity (Z c ) of the chemical in the space portion determined 
by expression (24) which will be mentioned later, and 
Fugacity capacity (Z 4 ) of the air determined by expression 
(4) of the first embodiment, the transference coefficient 
(D c4 ) of the compound between the chemical in the space 
portion of the carpet and the air is expressed as follows: 



From the above-mentioned transference speed (k c ) of 
the compound in the product, diffusion coefficient (D c ), 
area (A t ) of the spraying site, area (A cl ) in which the carpet 
fiber is in contact with the compound in the space portion, 




1 



...(22) 



l/(kA c A t Z c ) + lKk 4 AA t Z 4 ) 
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Fugacity capacity (Z c ) of the chemical in the space portion, 
and Fugacity capacity (Z 1 ) of the floor determined by 
expression (2) of the first embodiment, the transference 
coefficient (D cl ) of the compound between the chemical in 
5 the space portion and the carpet fiber is expressed as: 

D cl « ...(23) 

l/(MciAZc) + l/(MdA^i) 

wherein the transference speed (K 1 ) of the compound within 

yr§ the carpet fiber is expressed as (D c /t) 05 . 

q The volume (V c ) of the chemical in the space portion 

j=y 10 is determined from the area (A t ) of the spraying site, amount 

S of emission of the particle-constituting solvent (O a ), and 

q evaporation constant (R d ) of the chemical in the space 

^ portion determined by expression (20). The volume (V c ) of 

^ the space portion decreases until the solvent in the chemical 

'" : 15 completely evaporates. 

From the vapor pressure (P 3 ), water solubility (C 3 ), 

and octanol/water distribution coefficient (K ow ) of the 

compound, and property change constant (a) of the product, 

Fugacity capacity (Z c ) of the chemical in the space portion 

20 is expressed as: 

Z c ={e- at Kow + l-e- a ')C s IP 5 ...(24) 

Here, the property change constant (a) of the product 

can be computed according to the evaporation constant (R^) 

of the chemical in the space portion determined by expression 

25 (20), content (C 3d ) of the solvent (Xylene) contained in the 
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product , and dilution (X sol ) of the product* Namely, 
Fugacity capacity (Z c ) of the chemical in the space portion 
is also expressed as: 

Z c = [ (ratio of organic solvent in product) K ow + (ratio 
5 of water in product)] x C 3 /P 8 

By simultaneously forming this expression and the 
above-mentioned expression (24), the property change 
constant (a) of the product is computed. Here, in the case 

13 . 

of the oil-based product (aerosol), Fugacity capacity (Z c ) 
p 10 of the chemical in the space portion is given by expression 

m (5) of the first embodiment from the vapor pressure (P 3 ), 

|S water solubility (C s ), and octanol/water distribution 

q coefficient (K ow ) of the compound. 

i« Consequently, the behavior of the compound in the 

15 space portion of the carpet is expressed in the form of 

differential equation concerning Fugacity (f c ) as follows: 

4f. 



^V c Z c = y t Jj n-AVtZJ., Deposition^ 

+ R d A,Z c f c V- change 

- ^d(/c - /i) Transference^) 

- D c* ( fc - /a ) Transference^) 
-KV c ZJ c Degradation 

Also, expressions (15) and (17) in the first 
embodiment are respectively replaced by: 



.(25) 
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V — change 



3 



F - change 



+2 «A4v;Zi/, 



Deposition (i) 

Transference(c) 
Transferen ce(l) 
Degradation 



...(26) 



-*>«.(/.-/«) 
-A4(/,-A) 




-gv a zj a 



Ventilation 



3 



i 7 * - change 



Transferen ce(i) 



7 



-2 D M-fiy 



Transference^ j) 



...(27) 



-K 4 V 4 Z 4 f t 

-DM -A) 



Degradation 

Transference(c) 

Transference{X) 



In expression (26), terms of V-change, Deposition (i), 
Transference (c), Transference (1), and Degradation 
respectively indicate the change in volume (increase over 
time) of the carpet fiber, attachment accompanying the 
falling of particles, amount of transference of the compound 
between the carpet fiber and the space portion, amount of 
transference of the compound between the carpet fiber and 
the air, and change in amount of photodegradation of the 
compound • 

In expression (27), terms of Ventilation, 
Transference (i), Transference (j), Degradation, 
Transference (c), and Transference (1) respectively 
indicate the amount of discharge of the compound outdoors, 
amount of transference of the compound between the air and 
the particles, amount of transference of the compound 
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between the air and the floor/ceiling/wall, change in amount 
of photodegradation of the compound, amount of transference 
of the compound between the air and the space portion at 
the spraying site, and amount of transference of the compound 
5 between the air and the carpet fiber at the spraying site. 

Accordingly, Fugacity calculation at S3 is performed 
such that eight kinds of differential equations of 
expressions (16) and (18) in the first embodiment and the 

Q 

ij§ above-mentioned expressions of (25) to (27) are solved by 

H 10 Runge-Kutta-Gill method. The subsequent steps of S4 to S17 

v\} are similar to those in the first embodiment, 

jjj As explained in the foregoing, in the method of 

i=3 estimating an indoor behavior of a pesticidal compound in 

m accordance with this embodiment, in the case where the kind 

>j 15 of the spraying site is a carpet having ears of fiber and 

a space portion between the ears, a differential equation 
for Fugacity in the above-mentioned space portion of the 
carpet is added to seven kinds of simultaneous differential 
equations concerning Fugacity of the compound in the 
20 spraying site, two kinds of suspended particles, indoor air, 

floor, wall, and ceiling, so as to calculate each Fugacity. 

Consequently, even in the case where the chemical is 
accumulated in the space portion, so that the compound 
behaves differently from the case with a flat floor, the 
2 5 behavior of the compound can be estimated accurately, 

whereby various kinds of simulations can be performed. 
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( 3 ) Third Embodiment 

With reference to the drawings , the third embodiment 
of the present invention will be explained in the following. 

The method of simulating an indoor behavior of a 
pesticidal compound (hereinafter simply referred to as 
compound) contained in a chemical (solution) such as 
insecticide in accordance with this embodiment mainly 
estimates the behavior of the compound in the case where 
the chemical is spatially sprayed in order to terminate 
indoor vermin such as mosquito or fly. Here, the chemical 
encompasses all kinds of insecticides including pyrethroid 
insecticidal compounds and organophosphorus insecticidal 
compounds . 

First, as shown in Fig.. 11, an environment is 
classified into seven kinds of media consisting of large 
particles ( aerosol particles ) 2, medium particles (aerosol 
particles) 8, and small particles (aerosol particles) 3 
having diameters different from each other, whole indoor 
air 4, floor 5, wall 6, and ceiling 7. 

The particles are divided into three kinds since their 
behavior mainly depends on their diameter. Though the 
particles are distributed in the same site immediately after 
spraying, as time passes, the large particles 2 are 
distributed into a spraying zone SZ 2 near the floor due to 
their weight, the medium particles 8 into a spraying zone 
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SZ 8 in the middle of the room, and the small particles 3 into 
a spraying zone SZ 3 near the ceiling 7 since they are lighter 
than the other two kinds of particles . 

Simultaneous differential equations concerning 
Fugacity of the compound in the above-mentioned seven kinds 
of media are formed and are solved by Runge-Kutta-Gill method, 
whereby the indoor behavior of the compound in the chemical 
when the chemical is sprayed in the indoor space, i.e., 
temporal change of compound distribution in the room, is 
estimated and analyzed. 

Here, since the floor 5 mainly comprises three kinds, 
i.e., tatami, flooring, and carpet (rug), the differential 
equation is formed so as to correspond to the kind of floor. 
Explained in this embodiment is the case where the floor 
5 is made of tatami. 

With reference to Fig. 1, a method of simulating an 
indoor behavior of the compound will be explained. The 
simulation method in accordance with this embodiment can 
be mainly divided into a step of dividing an indoor 
environment into predetermined media and forming a 
differential equation concerning a fugacity of the compound 
in each medium (SI to S2) ; a step of determining the fugacity 
of the compound in each medium from the differential equation 

(53) ; a step of determining the indoor behavior of the 
compound from the fugacity of the compound in each medium 

(54) ; a step of changing, in response to a fluctuation in 
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b 10 

5 rzr 




mass balance of the compound indoors, a minute time unit 
used when solving the differential equation (S21 to S35, 
see Fig. 5) ; and a step of evaluating safety of the compound 
with respect to a human body according to the indoor behavior 
of the compound (S5 to S17). 

First, at step (hereinafter abridged as S) 1, a primary 
condition is inputted. The primary condition comprises 
physicochemical properties of the compound (see Table 16), 
indoor environmental behavior properties of the compound 
(see Table 17), indoor environment (see Table 18), product 
properties (see Table 19) , and a spray condition (see Table 
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Here, the examples of values in Tables 16 to 20 assume 
those of a six-mat room (9.72 [m 2 ]) of a typical apartment 
in Japan in summer (see Fig. 14). 

Subsequently, using the primary condition, a 
secondary condition is determined by calculation (S2) . The 
secondary condition comprises a temporally unchangeable 
definite factor condition (see Table 21) determined by the 
primary condition alone, and a temporal change condition 
accompanying a temporal change (see Table 22). Their 
calculations will be explained later in detail. 



# # 



TABLE 21 





DbrlNllb rAClOK CUND111UN 


SYMBOL 




NUMBER OF AEROSOL PARTICLES IMMEDIATELY AFTER SPRAYING 


nOi 0-2,8,3) 




EVAPORATION CONSTANT OF AEROSOL PARTICLES 


a 




VOLUME RATIO OF COMPOUND IN PRODUCT 


Ra 




DIFFUSION COEFFICIENT OF COMPOUND IN AIR 


Dca 




SURFACE AREA OF FLOOR / WALL / CEILING 


Aj (j=5,6,7) 




DIFFUSION COEFFICIENT OF FLOOR/ WALL / CEILING 


D c 




FUGACITY CAPACITY 






OF AEROSOL PARTICLE 


Zi (i=2,8,3) 




AIR 


Z 4 


J - 


FLOOR / WALL / CEILING 


Zj (j=5,6,7) 



no 



# • 



TABLE 22 



TEMPORAL CHANGE CONDITION 


SYMBOL 


AEROSOL PARTICLE 

DIAMETER 
FALLING SPEED 
FLOATING NUMBER 


di 0=2,8,3) 
vi 0=2,8,3) 
ni 0=2,8,3) 


AEROSOL ZONE WIDTH 


Lzi 0=2,8,3) 


TRANSFERENCE SPEED OF COMPOUND 

IN AEROSOL PARTICLE 
AIR 

FLOOR / WALL / CEILING 


ki 0=2,8,3) 
k 4 

kj 0=5,6,7) 


TRANSFERENCE COEFFICIENT OF COMPOUND 
BETWEEN AEROSOL PARTICLE AND AIR 

TRANSFERENCE COEFFICIENT OF COMPOUND 

BETWEEN AIR AND FLOOR / WALL /CEILING 


D i4 0=2,8,3) 
D 4 j 0=5,6,7) 


VOLUME OF AEROSOL PARTICLE 
VOLUME OF FLOOR / WALL / CEILING 


Vi (i-2,8,3) 
Vj 0=5,6,7) 
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Using the secondary condition, seven kinds of Fugacity 
are calculated (S3). Namely, seven kinds of differential 
equations concerning the large, medium, and small particles 
2, 8, and 3, air 4, floor 5, wall 6, and ceiling 7 are 
simultaneously formed and are solved by Runge-Kutta-Gill 
method, whereby seven kinds of Fugacity are computed over 
time. Here, an estimation nick time width (minute time 
unit) set when solving the differential equations is 
automatically set so as to be varied in response to a 
fluctuation in mass balance of the compound. 

At S4, using thus computed Fugacity (fi, i = 2, 8, 3) 
of the small, medium, and large particles 2, 8, and 3, and 
Fugacity (f 4 ) of the air 4, temporal concentration of the 
compound in the air is computed; whereas temporal floor 
residual amount of the compound on the floor 5 is computed 
by use of Fugacity (f 5 ) of the floor 5. 

At S5, it is judged whether or not to perform safety 
evaluation in the case where the chemical is inhaled. If 
the safety evaluation in inhalation is to be performed, then 
an estimated exposure amount in inhalation indicating a 
degree of exposure upon inhalation of the contaminated air 
is computed by use of the above-mentioned temporal 
concentration in the indoor air (S6). Thereafter, an 
inhalation safety coefficient is computed according to the 
estimated exposure amount in inhalation (S7). At S8, the 
inhalation safety coefficient is compared with a reference 
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value defined in each country. If the inhalation safety 
coefficient exceeds the reference value, it is judged that 
"there is no problem in safety." By contrast, if the 
inhalation safety coefficient is lower than the reference 
value at S8, it is judged that "there is a problem in safety, " 
and the operation returns to SI, where alteration of the 
primary condition such as alteration of compound, 
alteration of chemical formulation, alternation of using 
condition, or the like is considered. 

If the safety evaluation in inhalation is not selected 
at S5, calculated is an estimated amount of percutaneous 
exposure indicating to what extent the skin is exposed in 
contact with the floor 5 to which the chemical is attached 
(S9). Thereafter, at S10, it is judged whether to perform 
a percutaneous safety evaluation or not. If the 
percutaneous safety evaluation is to be performed, a 
percutaneous safety coefficient is computed according to 
the estimated percutaneous exposure amount (Sll). At S12, 
as with the safety evaluation in inhalation, the 
percutaneous safety coefficient is compared with a 
reference value defined in each country, whereby the safety 
is evaluated. If there is a problem in safety at S12, the 
operation returns to SI, where alteration of the primary 
condition is considered. 

If the percutaneous safety evaluation is not selected 
at S10, calculated at S13 according to the estimated 
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percutaneous exposure amount is an estimated oral exposure 
amount indicating the degree of exposure in the case where 
the chemical attached to a hand or the like is taken orally 

(513) . Subsequently, according to the estimated oral 

5 exposure amount, an oral safety coefficient is computed 

(514) . This oral exposure may occur, in particular, when 
an infant puts a chemical-attached hand into the mouth. At 
S15, as with the safety evaluation in inhalation, the oral 

£j safety coefficient is compared with a reference value 

□ 10 defined in each country, whereby the safety is evaluated, 

ij If there is a problem in safety at S15, the operation returns 

fi to SI, where alteration of the primary condition is 

g considered - 

^ Finally, if it is judged to be safe at S8, S12, or S15, 

15 overall safety is evaluated. Here, the sum of respective 

reciprocals of the previously determined inhalation safety 
coefficient, percutaneous safety coefficient, and oral 
safety coefficient is determined, and the reciprocal of thus 
determined value is defined as an overall safety coefficient 
20 (S16). This overall safety coefficient is compared with a 

reference value, whereby an evaluation similar to the 
previous safety evaluation is effected (S17). 

In the following, the above-mentioned steps of S2 to 
S4, S6, S7, S9, Sll, S13, and S14 will be explained in detail . 
25 (i) Secondary Condition Calculation (S2) 

Calculations of the definite factor condition shown 
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b 1° 



iri 

in 15 



20 



in Table 21 and the temporal change condition shown in Table 
22 will be explained. 

The particle number (n oi/ i = 2, 8, 3) of the large, 
medium, and small particles 2, 8, and 3 immediately after 
spraying is determined from the spraying rate, particle 
diameter (d oi , i = 2, 8, 3) immediately after spraying, 
distributions of particles (three kinds), and spray time. 
The particle number (n oi ) is constant until the bottom of 
the spraying zone (SZ if i = 2, 8, 3) reaches the floor 5. 

According to surface temperature (T d ) of particles, 
properties of particle-constituting solvent (vapor 
pressure P d , molecular weight M d , and specific gravity P d ) f 
and indoor environment (air diffusion coefficient D air , vapor 
pressure P°o, and room temperature T°o), the evaporation 
constant ( a ) of the large, medium, and small particles 2, 
8, and 3 is defined as: 



wherein R is a gas constant . 

Here, the upper parts of the above-mentioned Tables 
18 and 19 exemplify the case where the spatially sprayed 
aerosol is of oil base. In this case, Too (temperature at 
the site far from aerosol particles) and T d are set to room 
temperature (T), whereas poo (vapor pressure of oil at the 
site far from aerosol particles) is set to zero. 



In the case where the spatially sprayed aerosol is of 
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S 15 



water base, it is necessary for P d , M^, p d , T d , and poo in 
expression (28) to be changed to values based on properties 
of water. Specifically, the above-mentioned factor values 
when the room temperature (T) is 298 [K] (= 25°C) and the 
relative humidity is 60[%RH] are exemplified in the lower 
parts of Tables 18 and 19. 

Here, T d can be determined, according to the room 
temperature (T) and relative humidity (H), from "mass-based 
humidity table" disclosed in "Kagaku Kikai no Riron to Keisan 
(Theory and Calculation of Chemical Machines)" (Second 
Edition) (Saburo Kamei ed., Sangyo Tosho) or the like. On 
the other hand, P d and poo can be computed from the following 
expressions : 

log 10 P d = 10.23 - 1750/(T d -38) 

log lo poo = 8.23 - 1750/(T-38) + log lo 0 
wherein 0 is an indoor relative humidity. 

From the vapor pressure (P 3 ) and melting point ( T M ) 
of the compound and surface temperature (T d ) of particles, 
Fugacity capacity (Z if i = 2, 8, 3) of the large, medium, 
and small particles 2, 8, and 3 is expressed as: 



Though the vapor pressure of the compound in a liquid state 
(P L S ) is computed by use of T M and T d here; in the case where 
T M cannot be obtained, P L 3 may be set identical to the vapor 
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pressure of the compound in a solid state (P a ) . 

From the room temperature (T) , Fugacity capacity (Z 4 ) 
of the air 4 is expressed as follows : 

•■•(*>> 

5 From the vapor pressure (P s ), water solubility (C s ), 

and octanol/water distribution coefficient (K ow ) of the 
compound, and from the oil component content ( p^, j = 5, 

3 3 

■;Q 6,7) of the material constituting the floor 5, wall 6, and 

g ceiling 7, Fugacity capacity (Z jf j = 5, 6, 7 ) of the floor 

i=y 10 5, wall 6, and ceiling 7 is expressed as: 

jyj Z } = Pj K ow C s /P s ...(31) 

q Since the volume of particles decreases over time due 

i«j to evaporation of the particle-constituting solvent 

;:~ component, during the evaporation of particle-constituting 

15 solvent component, using the diameter (d oi ) of particles 

immediately after spraying (t = 0) and the evaporation 
constant ( a ) determined from expression (28) , the diameter 
(di, i = 2, 8, 3) of the large, medium, and small particles 
2,8, and 3 at time t is represented as shown in the following 
20 expression (32): 



d.^^dl-lat ...(32) 

By contrast, the diameter (di) of particles after 
complete evaporation of the particle-constituting solvent 
2 5 component is represented as shown in the above expression 
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(32 % ) by use of the ratio of volume of the compound in the 
product (R a ) and the diameter (d oi ) of particles immediately 
after spraying. 

From the compound content (C a ) and specific gravity 
of the compound ( p ) , the ratio of volume of the compound 
in the product (R a ) is expressed as C a /p. 

Since movement of the particles is under the control 
of gravity and air resistance, according to Stokes ' law using 
kQ the specific gravity ( p d ) of the particle-constituting 

q 10 solvent and the diameter (d ± ) of particles determined by 

iy expression (32), the falling speed (v i7 i = 2, 8, 3) of the 

SSj large, medium, and small particles 2, 8, and 3 before 

^ complete evaporation of the particle-constituting solvent 

==2j component is represented by the following expression (33): 

fi 15 v t -&&-di-pd}-fHd*-2ai) ...(33) 

18?7 

S c =1+ — {6.32 + 2.01exp(-8.322xl0 6 rf.)} ...(33') 

7.6x10 d. 

^-(plPtW? ••(33") 

wherein g is gravitational acceleration, v is the 
coefficient of viscosity of the air 4, S c is a sliding 

20 correction coefficient , and/? is a speed coefficient . Here, 

while particles do not conform to the Stokes 1 law when they 
become small, the sliding correction coefficient (S c ) , which 
is a coefficient for correcting this phenomenon, is 
represented by the above expression (33 1 ) according to the 
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diameter (di) of particles determined by expression (32). 

On the other hand, the falling speed (Vi) of the 
particles after the complete evaporation of particle- 
constituting solvent component is represented by the above 
5 expression (33") using the specific gravity ( p ) of the 

compound, specific gravity ( P d ) of the particle- 
constituting solvent, speed coefficient ( j3 ) , and diameter 
(d £ ) of particles determined by expression (32 f ). 

From the spray height, aerosol zone width (L 2i ) 
10 determined by expression (35) which will be mentioned later, 

falling speed (v £ ) of particles determined by expression 
(33) or expression (33"), and the above-mentioned particle 
number (n oi ) immediately after spraying, the floating number 
(n if i = 2, 8, 3) of the large, medium, and small particles 
>j 15 2, 8, and 3 after the bottom of the spraying zone (SZi) reaches 

a 

the floor 5 is expressed as: 
»:="o$-^-t Xi )-Gt] ...(34) 

wherein t xi is the time necessary for the bottom of the 
spraying zone (SZ L ) to reach the floor 5, and t < t xi . Here, 

20 G is the ventilation rate, taking account of the fact that 

the large, medium, and small particles 2, 8, and 3 not only 
disappear as attaching to the floor, wall, and ceiling but 
also are discharged outdoors by ventilation. Also, L zi is 
the dispersion width of suspended particles upon spraying, 

25 which may be either determined by expression (35) or 
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simplified by use of the dispersion width immediately after 
dispersion as it is. 

The width of aerosol zone (L zi , i = 2, 8, 3) is the 
width (height) of the spraying zone after the bottom of the 
spraying zone (SZ^ reaches the floor 5. Since the spraying 
zone (S zi ) is absorbed by the floor 5 while changing the 
falling speed (v ± ) of particles determined by expression 
(33) or expression (33"), from the spray height and the 
above-mentioned falling speed (Vi), the aerosol zone width 
is defined as: 



wherein t yi is the time necessary for the spraying zone (SZ ± ) 
to be completely absorbed by the floor 5. 

From the transference speed (k if i = 2, 8, 3) of the 
compound in the large, medium, and small particles 2, 8, 
and 3, transference speed (k 4 ) of the compound in the air 
4, diameter (d L ) determined by expression (32) or (32 1 ), and 
Fugacity capacity (Z L ) and Fugacity capacity (Z 4 ) determined 
by expressions (29) and (30), the transference coefficient 
(D i4 , i = 2, 8, 3) between the large, medium, and small 
particles 2,8, and 3 and the air 4 is represented as follows: 



wherein Ai (= 7Tdi 2 ) is the surface area of particles. From 
the ventilation rate (G) , the transference speed ( k 4 ) of the 




1 



...(36) 



l/(Wi) + l/{(*4 + ^iM^ 4 } 
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compound in the air 4 is expressed as GV 4 /A 4 . From the 
transference speed (k 4 ) of the compound in the air, the 
transference speed (k £ ) of the compound in the particles is 
expressed as k 4 /100. Here, as shown in Fig. 14, the 
above-mentioned A 4 is the cross-sectional area of the room 
with respect to the direction of movement of the air 4 
(arrowed direction in the drawing). 

Using the diffusion coefficient (D c ) of compound in 
the floor 5, wall 6, and ceiling 7, the transference speed 
(kj) of the compound in the floor 5, wall 6, and ceiling 7 
is expressed as (D c /t) 0 " 5 . Here, D c can be computed as follows . 
Namely, D c can be determined when the diffusion coefficient 
(D ca ) in the air is multiplied by 10" 9 , whereas D ca can be 
computed from properties of the compound (e.g., structural 
formula, molecular weight, and the like) according to Wike 
and Lee method ( "Handbook of Chemical Property Estimation 
Methods," McGraw-Hill Book Company, 1982), for example. 

From the half-life of transference to air (T 4j , j = 

5, 6, 7) in each of the floor 5, wall 6, and ceiling 7, room 
size (V 4 ), volume (V-) of the floor 5, wall 6, and ceiling 
7 determined by expression (39) mentioned later, and 
Fugacity capacity (Z 4 ) and (Zj) determined by expressions 
(30) and (31), the transference coefficient (D 4j , j= 5, 6, 
7) of the compound between the air 4 and the floor 5, wall 

6, and ceiling 7 is represented as the following expression 
(37): 
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D " %,{i/(v < zX/(^z i )} " (37) 

D " " 1/(M;Z.)+1/M/ Z /) (37 ) 

The above-mentioned transference coefficient (D 4j ) 
may also be determined by the expression represented by the 
above-mentioned expression (37') from the speed of the 
5 compound in the air (k 4 ), speed (kj) of the compound in the 

floor 5, wall 6, and ceiling 7, surface area (A jf j = 5, 6, 
7) of the floor 5, wall 6, and ceiling 7, and Fugacity 
capacity (Z 4 ) and (2 5 ) determined by expressions (30) and 
(31). The surface area (A 5 ) of each of the floor 5, wall 
10 6, and ceiling 7 is determined by the size of the room (V 4 

= L 4 x W 4 x height). 

Here, the half-life of transference to air (T 4j ) in 
expression (37) is hard to actually measure and varies under 
various conditions, thus making it difficult to determine 

15 an accurate value thereof. By contrast, all the parameters 

in expression (37 ' ) can be easily determined by calculation 
without necessitating measurement. Accordingly, the 
transference coefficient (D 4j ) is more preferably determined 
by expression (37') than by expression (37). 

2 0 Using the diameter (d^ of particles determined by 

expression (32) or expression (32'), the volume (V £ , i = 2, 
8, 3) of the large, medium, and small particles 2, 8, and 
3 is expressed by : 
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V = -d 3 .»(38) 
' 6 1 y ' 

Assuming that each of the floor 5, wall 6, and ceiling 
7 before spraying the chemical is like a thin film and that, 
as the chemical infiltrates into the film, its thickness 
increases so as to enhance the volume thereof; from the 
surface area (A t ) of the floor 5, wall 6, and ceiling 7 and 
diffusion coefficient (D c ) of the compound in the floor 5, 
wall 6, and ceiling 7, the volume (Vj, j = 5, 6, 7) of floor 
5, wall 6, and ceiling 7 becomes: 
Vj-ljD^Aj -(39) 

(ii) Fugacity Calculation (S3) 

The behavior of the compound in the large, medium, and 
small particles 2, 8, and 3 is expressed in the form of 
differential equation concerning Fugacity (fi) as: 

^-KZ- = -ad.Z J. V- change 

F - change - D i4 {f t - f 4 ) Transference^) 

- KV i Z i £ Degradation • • • (40) 

Here, the terms of V-change, Transference (4), and 

Degradation respectively indicate volume change (decrease 

over time) of the large, medium, and small particles 2, 8, 

and 3, amount of transference of the compound between the 

particles and the air 4, and change in amount of 

photodegradation of the compound- Here, among the 

photodegradation constants shown in the above-mentioned 

Table 17, it is difficult to actually measure the 
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photodegradation constant (K 4 ) of the compound in the air, 
whereby it is not usually easy to determine this value. 
Therefore, when the actually measured value is unavailable, 
it may be determined as being computed by Atmospheric 
Oxidation program (Atkinson et al., 1984, Chem. Rev* Vol, 
84, pp. 437-470) using the structural formula of the compound. 
Also, the degradation constant in the particles and in the 
interior materials (floor 5, wall 6, and ceiling 7) may be 
set to zero in view of the safety of inhabitants when the 
actually measured values are unavailable. Here, the term 
of V-change is effective only while the particle- 
constituting solvent component is evaporating, whereas 
Fugacity (f L ) is effective only while the spraying zone (SZi) 
exists . 

The behavior of the compound in the air 4 is expressed 
in the form of differential equation concerning Fugacity 
(f 4 ) as: 

^V 4 Z 4 = -GV 4 Z 4 f 4 Ventilation 
at 

F - change - ^n i D i4 (f 4 - /. ) Transference^) 

i- 2.3,3 

~ J D 4j if a - fj ) Transference^) 
-K 4 V 4 ZJ 4 Degradation •••(41) 

Here, the terms of Ventilation, Transference (i), 
Transference (j), and Degradation respectively indicate 
amount of discharge of the compound outdoors, amount of 
transference of the compound between the air 4 and the 
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particles, amount of transference between the air 4 and the 
floor 5, wall 6, and ceiling 7, and change in amount of 
photodegradation of the compound. 

The behavior of the compound in the floor 5, wall 6, 
and ceiling 7 is expressed in the form of differential 
equation concerning Fugacity (f 5 ), (f 6 )r (£7) as : 



V 5 Z 5 



dt 

F - change 



-jDjtA 5 Z 5 f s 

-D A5 {f 5 -h) 
-K 5 V 5 Z 5 f 5 



V - change 

Deposition(i) 

Transference^) 
Degradation 



•(42) 



^v 6 z 6 
dt 6 6 

F - change 



= -4DjtA 6 ZJ 6 



+ y 6 y ^w. 

-*>«(/.-/,) 

-k 6 v 6 zj 6 



V - change 

Deposition(i) 

Transference^) 
Degradation 



•(43) 



10 



^v 7 z 7 
dt 7 7 



F - change + y 1 



n n 



-K 7 V 7 Z 7 / 7 



V - change 

Deposition{i) 

Transfer en ce(4) 
Degradation 



•(44) 



Here, the terms of V-change, Transference (4), and 
Degradation in each expression respectively indicate volume 
change (increase over time) of the floor 5, wall 6, and 
ceiling 7, amount of transference of the compound between 
the floor 5, wall 6, ceiling 7 and the air 4, and change 
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in amount of photodegradation of the compound. Also, 
Deposition (i) indicates attachment accompanying the 
falling of the particles* Here, y 5 , y 6 , and y 7 indicate 
ratios of the large, medium, and small particles 2, 8, 3 
attaching to the floor 5 , wall 6 , and ceiling 7 , respectively , 
and mainly depend on the ventilation rate (e.g., when the 
ventilation rate is 1.58 [1/hour], y s = 0.93, y 6 = 0.06, 
7 7 = 0.01). As y 5 , y 6 , and y 7 are introduced, the behavior 
of the compound can be estimated accurately. When y 5f y 
6 , and y 7 are unavailable, it is preferable to presume the 
worst case in terms of safety evaluation and perform 

calculation with 7 5 =1, 7 6 ~ °/ and 7 7 = 0. 

The above-mentioned seven kinds of differential 
equations (40) to (44) are simultaneously formed and are 
solved by Runge-Kutta-Gill method, so as to compute Fugacity 
(f 2 to f 8 ). 

When solving these simultaneous differential 
equations, it is necessary to set an estimation nick time 
width (dt) which is a minute time unit. Namely, the 
estimation nick time width is used such that solutions of 
the simultaneous differential equations are initially 
determined at a time (t 0 ), and then solutions of the 
simultaneous differential equations are determined at a 
time (t 0 + dt) to which the estimation nick time width is 
added . As solutions are obtained while estimation nick time 
widths are successively added, temporally changing Fugacity 
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can be determined. Theoretically, as the set time of the 
estimation nick time width is shorter, more accurate 
solutions can be obtained, though necessitating a very long 
calculation time. By contrast, when the set time is too long, 
solutions tend to diverge, thereby generating errors. 

Therefore, in the present invention, the estimation 
nick time width is set shorter when a very large change occurs 
in a chemical, whereas it is set longer when there is no 
large change. 

Specifically, mass balance is always confirmed such 
that the amount of input of the chemical and the resulting 
solution coincide with each other, and the estimation nick 
time width is set longer when the mass balance does not 
fluctuate greatly, whereas it is set shorter when the mass 
balance starts fluctuating. For example, when the 
fluctuation of mass balance is set to an accuracy of ±5%, 
the estimation nick time width is always set so as to 
constantly satisfy the relationship of: 

compound input amount/ ( existing amount + degrading 
amount + discharging amount) = 0.95 to 1.05 
namely, such that the fluctuation of mass balance lies within 
the range of ±5%. 

Here, the above-mentioned compound input amount is 
determined by the compound content (C a ), spray rate, and 
spray time. Since the temporal amounts of compound in the 
seven kinds of media are determined by the simultaneous 
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differential equations, they are summed up so as to compute 
the existing amount as follows: 



dischargin g amount = GV 4 Z 4 f 4 

The degrading amount and discharging amount are as 
represented above . 

With reference to the flowchart of Fig. 5, the method 
of setting the estimation nick time width (dt) will be 
explained . 

First, an initial value of estimation nick time width 
(dt ) is inputted (S21 ) . Then, an upper limit set value (e.g. , 
0.1[%]) which is the upper limit of difference in mass 
balance, and a lower limit set value (e.g., 10" 6 [%]) which 
is the lower limit of difference in mass balance are inputted 
(S22). Thereafter, Fugacity and mass balance at t = t 0 are 
calculated (S23 and S24), and Fugacity and mass balance at 
t = t + dt (or t 0 + dt) are calculated (S25 and S26). 

It is judged whether the fluctuation in mass balance 
is within the range of ±5% or not (S27) . If the fluctuation 
in mass balance is within the range of ±5%, it is judged 
whether the difference between the mass balance at t = t 
(or t 0 ) and the mass balance at t = t + dt (or t fl + dt) is 
at least the upper limit set value or not (S28). If it is 



7 



existing amount = 




7 



degrading amount 




•■•(45) 
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judged to be at least the upper limit set value at S28, then 
solutions become more accurate when the estimation nick time 
width (dt) is made shorter. In this case, the estimation 
nick time width (dt) is multiplied by 1/2 so as to change 
its setting (S29). When the difference is judged to be 
smaller than the upper limit set value at S28, it is judged 
at S3 0 whether the difference between the mass balance at 
t = t (or t 0 ) and the mass balance at t = t + dt (or t 0 + 
dt) is at most the lower limit set value or not. 

When the mass balance difference is not greater than 
the lower limit set value at S3 0, since solutions are not 
influenced by longer estimation nick time width (dt), the 
estimation nick time width (dt) is doubled so as to change 
its setting (S3-1). Subsequently, it is judged whether the 
estimation nick time width (dt) changed at S31 is at most 
a maximum value (e.g., 0.1 [hour]) of estimation nick time 
width ( dt ) or not ( S3 2 ) . When the estimation nick time width 
(dt) is not greater than the maximum value at S3 2, since 
solutions do not diverge, the estimation nick time width 
(dt) set at S31 is used. When the estimation nick time width 
(dt) is greater than the maximum value at S3 2, since 
solutions may diverge, the estimation nick time width (dt) 
is reset to the maximum value (S33) . When the mass balance 
difference is greater than the lower limit set value at the 
above-mentioned S3 0, namely, when it lies between the lower 
limit set value and upper limit set value, the calculations 
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are continued without changing the estimation nick time 
width (dt). 

After the step of S29, S30, S32, or S33, the operation 
returns to S25 so as to effect calculation again, and this 
process is repeated till the aimed time is attained. 

When the fluctuation in mass balance exceeds the range 
of ±5% at S27 , on the other hand, the mass balance fluctuation 
is so much that calculation is preferably effected with an 
estimation nick time width (dt) shorter than that in the 
case where the fluctuation is within the range of ±5%. 
Accordingly, the calculation is stopped once (S34), the 
lower limit set value is reset to a lower level (S35), and 
then the operation returns to the step of S23. 

Thus, when the estimation nick time width is variably 
set without being held constant, while monitoring the mass 
balance fluctuation, Fugacity can be computed accurately 
and efficiently. 

(iii) Computation of Temporal Concentration in Indoor Air 
and Floor Residual Amount (S4) 

The temporal concentration of the compound in the 
indoor air is computed when Fugacity (f L ) of the large, medium, 
and small particles 2, 8, and 3, determined by the 
above-mentioned item (ii), multiplied by their Fugacity 
capacity (Z L ) , and Fugacity ( f 4 ) of the air 4 multiplied by 
Fugacity capacity (Z 4 ) of the air 4 are summed up. Here, 
the large particles 2 or medium particles 8 may not be inhaled 
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by a human body depending on the kind of chemical. In such 
a case, calculation is effected by use of Fugacity (f 3 ) of 
the small particles 3 alone , for example. 

The floor residual amount of the compound is computed 
when Fugacity (f s ) of the floor 5 is multiplied by Fugacity 
capacity (Z s ) of the floor 5. 

(iv) Calculation of Estimated Exposure Amount in Inhalation 
and Inhalation Safety Coefficient (S6 and S7) 

The above-mentioned temporal concentration in the 
indoor air forms a curve shown in Fig. 3A, for example. This 
concentration curve is integrated, an accumulated 
concentration of the compound during a specific period (t x 
to t 2 ) is determined (see Fig. 3B), and the mean concentration 
in the indoor air is computed from thus determined value. 
While an arbitrary period is set as the specific period 
depending on the object, an appropriate period is usually 
set in view of the method of use of the product and the test 
period of toxicity data. 

Then, according to the above-mentioned mean 
concentration in the indoor air, amount of respiration, and 
exposure time, the estimated exposure amount in inhalation 
is determined. Namely, calculation of: 

estimated exposure amount in inhalation [mg/kg/day] 
= mean concentration in indoor air [mg/m 3 ] x amount of 
respiration [m 3 /kg/min] x exposure time [min/day] 
is effected. Here, as the above-mentioned amount of 
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respiration, a published value or actually measured value 
may be used. Also, when amounts of respiration are 
respectively set for adult and child, more appropriate 
estimated exposure amounts in inhalation can be obtained. 
In the case where the inhaled harmful material is not totally 
absorbed into the body but is partially discharged by 
respiration, a more appropriate estimated exposure amount 
in inhalation can be obtained when the inhalation ratio is 
taken into account . 

The inhalation safety coefficient is computed from a 
non-influential amount concerning inhalation toxicity 
examined by an animal experiment beforehand and the 
estimated exposure amount in inhalation determined above. 
Namely, it is expressed as: 

inhalation safety coefficient = 

inhalation non-influential amount [mg/kg/day ] / 
estimated exposure amount in inhalation [mg/kg/day] 
(v) Calculation of Estimated Percutaneous Exposure Amount 
and Percutaneous Safety Coefficient (S9 and Sll) 

The above-mentioned floor residual amount forms a 
curve such as that shown in Fig. 4A, for example. This 
residual amount curve is integrated, the accumulated 
residual amount of the compound during a specific period 
(t x to t 2 ) is determined (see Fig. 4B), and the mean residual 
amount is computed from thus determined value. While an 
arbitrary period is set as the specific period depending 
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on the object, an appropriate period is usually set in view 
of the method of use of the product and the test period of 
toxicity data. 

Then, according to the mean floor residual amount, 
skin attachment ratio, contact area, and body weight, the 
estimated percutaneous exposure amount is determined. 
Namely, calculation of: 

estimated percutaneous exposure amount [mg/kg/day] = 
(mean floor residual amount [mg/m 2 ] x skin attachment ratio 
[%] x contact area [mVday ] ) /body weight [kg] 
is effected. Here, as the contact area, a published value 
(e.g., 4 [mVday]) may be used. The skin attachment ratio 
is a ratio of the compound attaching to the skin when the 
latter is in contact with the floor 5 where the compound 
exists. As this value, a published value or a value 
experimentally obtained from a model may be used. 

A model experiment method for the skin attachment 
ratio is as follows. A weight (8 cm x 8 cm x 8 cm; 4.2 kg) 
is placed on a denim cloth (8 cm x io cm) with a pressure 
similar to that of an infant in contact with a floor, and 
the denim cloth is pulled on the floor at a speed (120 cm/15 
sec) similar to the moving speed of the infant. The denim 
and floor are analyzed so as to compute the compound 
contained in the denim and floor. From the ratio 
therebetween, the skin attachment ratio is obtained. It has 
been confirmed that the skin attachment ratio obtained by 
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this method is identical to or slightly higher than that 
determined from analyzed values of a hand and a floor when 
the hand is actually pressed against the floor, thereby 
proving this model experiment method to be useful for 
5 evaluating exposure of inhabitants. 

The percutaneous safety coefficient is computed from 
the non-influential amount concerning percutaneous 
toxicity examined by an animal experiment beforehand and 
the estimated percutaneous exposure amount determined above. 
Q 10 Namely, it is expressed as: 

fy percutaneous safety coefficient = 

iji percutaneous non-influential amount [mg/kg/day ] / 

□ estimated percutaneous exposure amount [mg/kg/day] 

H Nevertheless, in general, percutaneous non- 

J, 15 influential amount has not often been determined, and there 

It* 

ls ~ are not many published values . Accordingly, a more accurate 

value can be determined from the estimated percutaneous 
exposure amount, and oral non-influential amount and 
percutaneous absorption ratio for which many published 
20 values exist, according to the following expression: 

percutaneous safety coefficient = oral non- 
influential amount [mg/kg/day ]/( estimated percutaneous 
exposure amount [mg/kg/day] x percutaneous absorption ratio 

[%]) 

25 Here, when the percutaneous absorption ratio is unknown, 

employed is a national guideline (e.g., 10%) which usually 



134 



SCC 97-06 



exists . 

(vi) Calculation of Estimated Oral Exposure Amount and Oral 
Safety Coefficient (S13 and S14) 

From the estimated percutaneous exposure amount 
obtained in the above-mentioned item (v) , hand surface area 
ratio, and oral transference ratio, the estimated oral 
exposure amount from hand to mouth is determined. Namely, 
calculation of: 

estimated oral exposure amount [mg/kg/day] = 
estimated percutaneous exposure amount [mg/kg/day] * hand 
surface area ratio [%] * oral transference ratio [%] 
is effected. Here, the hand surface area ratio is expressed 
by (hand surface area/body surface area), for which a 
published value (e.g., 5 ["%']) may be used. The oral 
transference ratio is a hypothetical value, which is set 
to 100%, for example. 

In the case where oral exposure might occur via 
tableware or food contaminated with the residually sprayed 
compound, it is required that the estimated oral exposure 
amount from tableware or food to mouth be added to the 
estimated oral exposure amount from hand to mouth to yield 
the total estimated oral exposure amount. For example, the 
estimated oral exposure amount from tableware is obtained 
when, according to the tableware residual amount indicating 
the amount of the harmful material remaining in tableware, 
tableware use area which is the sum of tableware surface 
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areas, and oral transference ratio from tableware, 
calculation of: 

estimated oral exposure amount [mg/kg/day] = 
tableware residual amount [mg/m 2 ] * tableware use area 
[mVday] * oral transference ratio [%]/body weight [kg] 
is effected. Here, the tableware residual amount is 
expressed by (mean floor residual amount x tableware 
contamination ratio ) . As the tableware contamination ratio , 
an actually measured value (e.g., 9% ) or a hypothetical value 
may be used. 

The oral safety coefficient is computed from the 
non-influential amount concerning oral toxicity examined 
by an animal experiment beforehand and the estimated oral 
exposure amount determined above. Namely, it is expressed 
as : 

oral safety coefficient = 

oral non-influential amount [mg/kg/day]/ estimated 
oral exposure amount [mg/kg/day] 

As mentioned in the foregoing, in the method of 
estimating an indoor behavior of a pesticidal compound in 
this embodiment, when a chemical containing the above- 
mentioned compound is spatially sprayed, differential 
equations concerning Fugacity of the large, medium, and 
small particles 2, 8, and 3, air 4, floor 5, wall 6, and 
ceiling 7 are simultaneously formed and are solved, and the 
indoor behavior of the compound is estimated according to 
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thus obtained solution. Here, the estimation nick time 
width is variably set, while constantly confirming mass 
balance of the compound indoors after the spraying, so that 
the amount of input of the chemical indoors and the resulting 
5 solution coincide with each other. 

Accordingly, since mass balance of the compound after 
the spraying is always monitored such that the amount of 
input of the compound indoors and the resulting solution 
n coincide with each other, thereby varying the estimation 

□ 10 nick time width; the estimation nick time width is set longer 

y when the mass balance fluctuates a little, whereas it is 

fi set shorter when the mass balance starts fluctuating greatly . 

3 Namely, when solving simultaneous differential equations 

~j including a parameter accompanying temporal change, the 

15 estimation nick time width is automatically set in response 

to the fluctuation in mass balance. Consequently, when 
processed by a computer, an accurate solution can be obtained 
in a short time. 

The method of evaluating safety of a pesticidal 
20 compound in accordance with this embodiment uses the 

estimated result mentioned above to evaluate the safety of 
the compound with respect to the human body when the chemical 
is spatially sprayed. 

Accordingly, the safety of the compound with respect 
25 to the human body can be evaluated accurately in a short 

time. As a consequence, when formulating a chemical such 
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as insecticide containing the compound, simulation can be 
easily repeated while changing conditions, thereby making 
it easier to formulate a chemical having a high safety 
conforming to the aimed object. 
5 Though the kind of the floor 5 is assumed to be tatami 

in this embodiment, differential equations similar to those 
mentioned above may also be formed in the case of flooring. 
In this case, however, when determining fugacity capacity 
D (Z 5 ) of the floor 5 in equation (31), it is necessary for 

q 10 the particle-constituting solvent component content (see 

ij Table 18) to be changed from that of tatami to that of 

yi flooring. 

3 Though Fugacity is determined by use of Runge- 

3 Kutta-Gill method in this embodiment, other methods may be 

l 15 used for solving differential equations. Runge-Kutta-Gill 

method, however, is preferably used since a program for the 
above-mentioned differential equations can be easily made 
by Basic. Also in the case where differential equations are 
solved by a method other than Runge-Kutta-Gill method, 
20 similar effects can be obtained when the estimation nick 

time width is set as mentioned above. 

In the following, examples (Examples 6 to 10) of the 
method of simulating an indoor behavior of a pesticidal 
compound in accordance with the third embodiment will be 
25 explained . 

The evaluation of the method of estimating an indoor 
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behavior of a pesticidal compound in these examples is 
effected, under various experimental conditions (Example 
6 to Example 10), as the results obtained by estimation and 
the actually measured values obtained by measurement of the 
5 actual indoor behavior of the compound are compared with 

each other. The respective experimental conditions of 
Examples are as follows. 

Example 6 uses, as the spatially sprayed chemical, an 
oil-based insecticide containing pyrethroid. The indoor 

10 environment is a typical apartment room (3.6 [m] x 3.6 [m] 

x 2.4 [m]) in Japan . The floor is a flooring in which a wooden 
floor is covered with a polyurethane resin, whereas the wall 
and ceiling are covered with a polychlorovinyl wallpaper. 
The windows are completely closed during and after the 

15 spatial spraying. The spatial spraying is uniformly 

effected for 10 seconds at a height of 1 . 6 [m] from the floor. 
The other conditions are shown in Tables 16 to 20. 

The conditions of Example 7 are the same as those of 
Example 6 except that its ventilation rate is 0.5 [1/hour]. 

20 The conditions of Example 8 are the same as those of 

Example 6 except that its ventilation rate is 4 . 1 [1/hour] . 

The conditions of Example 9 are the same as those of 
Example 6 except that the windows are opened for only 5 
minutes after 5 minutes from the spraying. 

25 The conditions of Example 10 are the same as those of 

Example 6 except that the windows are opened for only 2 hours 
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after 5 minutes from the spraying. 

Figs. 12A and 12B show the simulation results and 
actually measured values of the compound concentration in 
the air under the respective conditions of the above- 
mentioned Examples. Also, Fig. 12C shows the simulation 
results of the residual amounts in floor and ceiling in 
Example 6. In these charts, solid lines, dotted lines, and 
dashed lines indicate simulation results, whereas plotted 
points indicate actually measured values. 

From the above-mentioned results, it can be seen that 
the estimated results according to the method of simulating 
an indoor behavior of the compound in accordance with the 
present invention are very close to the actually measured 
values . 

( 4 ) Fourth Embodiment 

With reference to the drawings, the fourth embodiment 
of the present invention will be explained in the following. 
The method of simulating an indoor behavior of a pesticidal 
compound in accordance with this embodiment estimates the 
behavior of the compound in the case where the floor in the 
third embodiment is constituted by a rug having ears of fiber . 
For convenience of explanation, members identical to those 
shown in the drawings of the above-mentioned third 
embodiment will be referred to with numerals or letters 
identical thereto, without their explanations provided. 
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Explained here is a case where the kind of floor is a carpet. 
Accordingly, only the portion different from the third 
embodiment will be explained. 

First, the carpet is modeled as shown in Figs. 13A, 
13B, and 13C. Namely, the carpet is divided into a plurality 
of carpet fibers (ears of fiber) planted on a substrate and 
space portions existing between the fibers . When a chemical 
is spatially sprayed in a room in which such a carpet covers 
the whole floor surface, as particles fall on the carpet, 
the chemical infiltrates into the carpet fiber and 
accumulates in the space portion. 

While the method of simulating an indoor behavior of 
a pesticidal compound in accordance with this embodiment 
also conforms to the flowchart of Fig. 1, at SI, as the 
primary condition, the carpet-related primary condition 
shown in Table 23 is added to the primary condition of the 
third embodiment (see Tables 16 to 20). Here, the 
cross-sectional area (A c ) of the space portion of the carpet 
refers to the area in which the space portion is in contact 
with the indoor air (see Figs. 13A, 13B, and 13C). 
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The examples of values in Table 23 are based on the 
values shown in Figs, 13A, 13B, and 13C. Here, one space 
portion is set to a size of 0.1 [mm] * 0.1 [mm] * 3 [mm], 
and one square of carpet fiber is set to a size of 2.5 [mm] 
5 x 2.5 [mm] x 3 [mm]. Also, it is assumed that 16 [pieces] 

of carpet fiber exist in 1 [cm 2 ]. 

In this case, the individual parameters can be 
determined as : 

=n A c = ( 0 . 1 [mm] x 0.1 [mm] x 16 [pieces ])/l [cm] 

g 10 = 1.6 x io -3 [m 2 /m 2 floor] 

jy = (0.1 [mm] x 3 [mm] x 4 [faces] x 16 [pieces ])/l 

m [ CItl ] 

3 = 0.192 [m 2 /m 2 floor] 

1 A 5 = 1 - A c 

!n 15 = 1 - 1.6 x io -3 [m 2 /m 2 floor] 

A aS = A 5 + A c5 

= 1.19 [m 2 /m 2 floor] 
Here, A 5 is the area in which the carpet fibers in 1 [m 2 ] 
are in contact with the indoor air. 

20 As the secondary condition at S2, the carpet-related 

secondary condition shown in Table 24 is added to the 
secondary condition of the third embodiment. Here, k 5 and 
V 5 , which are respectively transference speed of the 
compound in the floor and the volume of the floor in the 

25 third embodiment (see Table 22), are respectively replaced 

by the transference speed of the compound in the carpet fiber 
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and the volume of the whole fiber of the carpet. 
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TABLE 24 





CARPET - RELATED SECONDARY CONDITION 


SYMBOL 




DEFINITE FACTOR CONDITION 






EVAPORATION CONSTANT OF CHEMICAL IN SPACE PORTION 


Rd 




HUMIDITY AT TEMPERATURE (T d ) OF SPACE PORTION 






VOLUME OF WHOLE FIBER OF CARPET 


v 5 




TRANSFERENCE SPEED OF COMPOUND IN PRODUCT 


k c 




PROPERTY CHANGE CONSTANT OF PRODUCT 


a 


3 


TEMPORAL CHANGE CONDITION 






TRANSFERENCE COEFFICIENT OF COMPOUND 






BETWEEN CHEMICAL IN SPACE PORTION AND AIR 


D c4 




TRANSFERENCE COEFFICIENT OF COMPOUND 




o 


BETWEEN CHEMICAL IN SPACE PORTION AND CARPET FIBER 


D C 5 




VOLUME OF CHEMICAL IN SPACE PORTION 


v c 




FUGACITY CAPACITY OF CHEMICAL IN SPACE PORTION 
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A method of computing the definite factor conditions 
and temporal change conditions shown in the above-mentioned 
Table 24 will be explained hereinafter as exemplified by 
the case using an aqueous aerosol, i.e., case where the 
particle-constituting solvent is water. 

From the cross-sectional area (A c ) of the space portion, 
room temperature (T), and humidity (H), the evaporation 
constant (B^) of the chemical entered into the space portion 
of the carpet is expressed as: 

R d =A c h{H m -H)/C H .-(46) 
C H « 0.24 + 0.46// 
h =3.06xl0" 4 (r-7;) 1/3 

wherein T d is the temperature of the space portion of the 
carpet (corresponding to the wet-bulb temperature), and H m 
is the humidity (corresponding to the saturated humidity 
at the wet-bulb temperature) of the space portion at T d . 
Here, T d and H m can be determined from the above-mentioned 
"mass-based humidity table" according to the room 
temperature (T) and humidity (H). In the case of an 
oil-based product (aerosol), while the vapor constant of 
the chemical can be computed when properties of the water 
and oil are compared with each other and then are multiplied 
by R d/ = 0 may be assumed in the case where the evaporation 
of the oil is slow. 

From the diffusion coefficient (D c ) and the area (A o5 ) 
in which the carpet fiber is in contact with the compound, 
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the volume (V 5 ) of the whole fiber of the carpet is determined 
as : 

V 5 =2jD^A a5 -(47) 

The transference speed (k c ) of the compound in the 
5 product is expressed as k 4 /100 by use of k 4 (= GV 4 /A 4 ) which 

indicates the transference speed of the compound in the 
indoor air. 

From the above-mentioned transference speed (k c ) of 
=n the compound in the product, transference speed (k 4 ) of the 

O 10 compound in the indoor air, floor size (L 4 * W 4 ), cross- 

m sectional area (A c ) of the space portion, Fugacity capacity 

m (2 C ) of the chemical in the space portion determined by 

g expression (50) which will be mentioned later, and Fugacity 

13 capacity (Z 4 ) of the air determined by expression (30) of 

'fj 15 the third embodiment, the transference coefficient (D c4 ) of 

"~ the compound between the chemical in the space portion of 

the carpet and the air is expressed as: 

D c4 = — ( X ***( 48 ) 

i/(UWc)+i/(WM) 

From the transference speed (k c ) of the compound in 
20 the product, diffusion coefficient (D c ), floor size (L 4 * 

W 4 ) , area (A c5 ) in which the carpet fiber is in contact with 
the compound in the space portion, Fugacity capacity (Z c ) 
of the chemical in the space portion, and Fugacity capacity 
(Z 5 ) of the floor determined by expression (31) of the third 
25 embodiment, the transference coefficient (D c5 ) of the 
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compound between the chemical in the space portion and the 
carpet fiber is expressed as: 



wherein the transference speed (k 5 ) of the compound within 
the carpet fiber is expressed as (D c /t) 05 . 

The volume (V c ) of the chemical in the space portion 
is determined from the floor size (L 4 x w 4 ), spray rate, 
content (C gd ) of the solvent (Xylene) contained in the 
product, dilution (X aol ) of the product, and evaporation 
constant (Rd) of the chemical in the space portion determined 
by expression (46). 

From the vapor pressure (P 3 ), water solubility (C s ), 
and octanol/water distribution coefficient (K ow ) of the 
compound, and property change constant (a) of the product, 
Fugacity capacity (Z c ) of the chemical in the space portion 
is expressed as : 



Here, the property change constant (a) of the product 
can be computed according to the evaporation constant (R^) 
of the chemical in the space portion determined by expression 
(46), content (C ad ) of the solvent (Xylene) contained in the 
product, and dilution (X gol ) of the product. Namely, 
Fugacity capacity (Z c ) of the chemical in the space portion 
is also expressed as: 

Z c = [ (ratio of organic solvent in product) K ow + (ratio 



1 



... (49 ) 



l/(U 5 Wc)+l/(Mc 5 W 5 ) 
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of water in product)] * C s /P 8 

By simultaneously forming this expression and the 
above-mentioned expression (50), the property change 
constant (a) of the product is computed. Here, in the case 
of the oil-based product (aerosol), Fugacity capacity (Z c ) 
of the chemical in the space portion is given by expression 
(31) of the third embodiment from the vapor pressure (P s ), 
water solubility (C s ), and octanol/water distribution 
coefficient (K ow ) of the compound. 

Consequently, the behavior of the compound in the 
space portion of the carpet is expressed in the form of 
differential equation concerning Fugacity (f c ) as: 

^t V < Z < = ^ J^^A^i/i IL * Deposition® 
F - change + R d L 4 W 4 Z c f c V - change 

- D c4 (fc ~ h ) Transference^) 
- D cs (fc " fs) Transference^) 

- KV c Z c f c Degradation -(51) 

Also, expressions (42) and (41) in the third embodiment are 
respectively replaced by: 

^V S Z S = jDjJA a5 Z 5 f 5 V - change 

dt 

F - change + y 5 (l - A c )— v.l/.Z . f. Deposition® 

~ D cs (fs ~ L ) Transference^) 

- D 45 (f 5 - f 4 ) Transference^) 

- K 5 V 5 Z 5 f 5 Degradation • *(52) 



SCC 97-06 



—V 4 Z 4 « -GV 4 Z 4 f 4 Ventilation 
dt 4 4 

F - change - V /i ; Z> 4 (/ 4 - / \) Transference^) 

;-2,8,3 

- ^£> 4 ;(/4 " /; ) Transference^]) 

- K 4 V 4 Z 4 f 4 Degradation 

- D 4c (f 4 - f c ) Transfer en ce(c) • • • (53) 

In expression ( 52 ) , terms of V-change, Deposition ( i) , 
Transference (c), Transference (4), and Degradation 
respectively indicate the change in volume (increase over 
g 5 time) of the carpet fiber, attachment accompanying the 

ly falling of particles, amount of transference of the compound 

5 - ! 

m between the carpet fiber and the space portion, amount of 

q transference of the compound between the carpet fiber and 

g the air, and change in amount of photodegradation of the 

10 compound. 



-•SsS 

IS 



In expression (53), terms of Ventilation, 
Transference (i), Transference (j), Degradation, and 
Transference (c) respectively indicate the amount of 
discharge of the compound outdoors, amount of transference 

15 of the compound between the air and the particles, amount 

of transference of the compound between the air and the floor 
(carpet fiber ) /ceiling/wall , change in amount of 
photodegradation of the compound, and amount of 
transference of the compound between the air and the space 

20 portion. 

Accordingly, Fugacity calculation at S3 is performed 
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such that eight kinds of differential equations of 
expressions (40), (43), and (44) in the third embodiment 
and the above-mentioned expressions (51) to (53) are solved 
by Runge-Kutta-Gill method. The subsequent steps of S4 to 
SI 7 are similar to those in the third embodiment. 

As explained in the foregoing, in the method of 
estimating an indoor behavior of a pesticidal compound in 
accordance with this embodiment, in the case where the kind 
of floor is a carpet having ears of fiber and a space portion 
between the ears, a differential equation for Fugacity in 
the above-mentioned space portion of the carpet is added 
to seven kinds of simultaneous differential equations 
concerning Fugacity in the three kinds of large, medium, 
and small particles, indoor air, floor (carpet fiber), 
ceiling, and wall, so as to calculate each Fugacity. 

Consequently, even in the case where the chemical is 
accumulated in the space portion, so that the compound 
behaves differently from the case with a flat floor, the 
behavior of the compound can be estimated accurately, 
whereby various kinds of simulations can be performed. 

(5) Fifth Embodiment 

With reference to the drawings, the fifth embodiment 
of the present invention will be explained in the following. 

The method of simulating an indoor behavior of a 
pesticidal compound (hereinafter simply referred to as 
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compound) contained in a chemical (solution) such as 
insecticide in accordance with this embodiment estimates, 
for example, the behavior of the compound in the case where 
the chemical is heated to vaporize by use of an electric 
heating vaporizer such as electronic mosquito repellent in 
order to terminate indoor vermin. Here, the chemical 
encompasses all kinds of insecticides including pyrethroid 
insecticidal compounds and organophosphorus insecticidal 
compounds . 

First, as shown in Fig. 15, an environment is 
classified into 11 kinds of media consisting of 3 kinds of 
condensed particles 9, 10, and 11, high-concentration air 
12, medium-concentration air 13, low-concentration air 14, 
floor 5, wall 6, high-concentration ceiling 15, medium- 
concentration ceiling 16, and low-concentration ceiling 17 . 

Namely, in the case where the chemical is vaporized, 
though the whole compound is initially evaporated as 
complete vapor (high-concentration air 12), the condensed 
particles 9, 10, and 11 are successively formed when the 
emission rate of the compound (E T ) in the medium- 
concentration air 13 exceeds the saturated amount of the 
compound in the vapor. These condensed particles 9, 10, and 
11 are divided into three kinds depending on the time between 
the occurrence and extinction of particles. The indoor air 
other than the above-mentioned high-concentration air 12 
and medium-concentration air 13, which is mainly respired 
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by humans, is defined as the low-concentration air 14. The 
high-concentration ceiling 15 is a ceiling portion into 
which the condensed particles 9, 10, and 11 are absorbed; 
the medium-concentration ceiling 16 is a ceiling portion 
in contact with the medium-concentration air 13, and the 
low-concentration ceiling 17 is a ceiling portion in contact 
with the low-concentration air 14 . 

Simultaneous differential equations concerning 
Fugacity of the compound in the above-mentioned 11 kinds 
of media are formed and are solved by Runge-Kutta-Gill method, 
whereby the indoor behavior of the compound in the chemical 
when the chemical is heated and vaporized, i.e., temporal 
change in compound distribution in the room, is estimated 
and analyzed. 

Here, since the floor 5 mainly comprises three kinds, 
i.e., tatami, flooring, and carpet (rug), the differential 
equation is formed so as to correspond to the kind of floor. 
Explained in this embodiment is the case where the floor 
5 is made of tatami. 

With reference to Fig. 1, a method of simulating an 
indoor behavior of the compound will be explained. The 
simulation method in accordance with this embodiment can 
be mainly divided into a step of dividing an indoor 
environment into predetermined media and forming a 
differential equation concerning a fugacity of the compound 
in each medium (SI to S2 ) ; a step of determining the fugacity 
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of the compound in each medium from the differential equation 

(53) ; a step of determining the indoor behavior of the 
compound from the fugacity of the compound in each medium 

(54) ; a step of changing, in response to a fluctuation in 
mass balance of the compound indoors, a minute time unit 
used when solving the differential equation (S21 to S35, 
see Fig. 5) ; and a step of evaluating safety of the compound 
with respect to a human body according to the indoor behavior 
of the compound (S5 to S17). 

First, at step (hereinafter abridged as S ) 1, a primary 
condition is inputted. The primary condition comprises 
physicochemical properties of the compound (see Table 25), 
indoor environmental behavior properties of the compound 
(see Table 26), indoor environment (see Table 27), and 
product properties (see Table 28). 
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Here, A 12 , A 13 , and A 14 in Table 28 indicate cross- 
sectional areas of the respective airs 12, 13, and 14 with 
respect to the direction of movement of indoor air (arrow 
in Fig . 15 ) . 

The examples of values in Tables 25 to 28 assume those 
of a six-mat room (9.72 [m 2 ] ) of a typical apartment in Japan 
in summer, using an electric heating vaporizer employing 
an oil-based chemical (see Figs. 16Aandl6B). The electric 
heating vaporizer comprises a core portion of 1 [cm] x l 
[cm] (at a temperature of 405 [K]), a heater portion of 4 
[cm] x 4 [cm] x 0.5 [cm] (at a temperature of 373 [K]), and 
a container in which a chemical is packed (at a temperature 
of 303 [K]). The example of value of compound content (R a ) 
in Table 2 8 indicates the volume of compound relative to 
the volume of particles immediately after condensation. 

Subsequently, using the primary condition, a 
secondary condition is determined by calculation (S2) . The 
secondary condition comprises a temporally unchangeable 
definite factor condition (see Table 29) determined by the 
primary condition alone and a temporal change condition 
accompanying a temporal change (see Table 30). Their 
calculations will be explained later in detail. 
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TABLE 29 



DEFINITE FACTOR CONDITION 


SYMBOL 


FLOATING NUMBER OF CONDENSED PARTICLES 
EVAPORATION CONSTANT OF CONDENSED PARTICLES 
EMISSION AMOUNT OF CHEMICAL IN GAS STATE 
DIFFUSION COEFFICIENT OF COMPOUND IN AIR 
SURFACE AREA OF FLOOR / WALL / CEILING 
DIFFUSION COEFFICIENT OF FLOOR / WALL / CEILING 
AIR VELOCITY 


ni (i=9,10,ll) 

a 

E v 

D ca 

A k (k=5,6,15,16,17) 

D c 

Vf 


FUGACITY CAPACITY 

OF CONDENSED PARTICLE 
AIR 

FLOOR / WALL / CEILING 


Zj (i=9,10,ll) 
Zj (j=12,13,14) 
Z k (k=5,6,15,16 J 17) 


DEGRADATION CONSTANT 

OF CONDENSED PARTICLE 
AIR 

FLOOR / WALL / CEILING 


Ki (i=9,10,ll) 
Kj 0=12,13,14) 
K k (k=5,6,15,16,17) 
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TABLE 30 





TEMPORAL CHANGE CONDITION 


SYMBOL 




CONDENSED PARTICLE DIAMETER 


di (i=9,10,ll) 




TRANSFERENCE SPEED OF COMPOUND 






IN CONDENSED PARTICLE 


ki 0=9,10,11) 




AIR 


kj (j=12,13,14) 


5 ! 

<0 


FLOOR / WALL / CEILING 


k k (k=5,6,15,16,17) 


? - 

Q 


TRANSFERENCE COEFFICIENT OF COMPOUND 






BETWEEN CONDENSED PARTICLE AND AIR 


Dil3 0=9,10,11) 


? ; \, 


TRANSFERENCE COEFFICIENT OF COMPOUND 






BETWEEN AIR AND FLOOR / WALL / CEILING 


Dj k 0=12,13,14, 


s : 




k=5,6,15,16,17) 


i _ j 

"SI 

3' i 


VOLUME OF CONDENSED PARTICLE 


Vi 0=9,10,11) 




VOLUME OF FLOOR / WALL / CEILING 


V k (k=5,6,15,16,17) 
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Using the secondary condition, 11 kinds of Fugacity 
are calculated (S3). Namely, 11 kinds of differential 
equations concerning the 3 kinds of condensed particles 9, 
10, and 11, 3 kinds of air 12, 13, and 14, floor 5, wall 
5 6, and 3 kinds of ceilings 15, 16, and 17 are simultaneously 

formed and are solved by Runge-Kutta-Gill method, whereby 
11 kinds of Fugacity are computed over time. Here, an 
estimation nick time width (minute time unit) set when 
i solving the differential equations is automatically set so 

i 10 as to be varied in response to a fluctuation in mass balance 

of the compound. 

ifl At S4, using thus computed Fugacity (f 14 ) of the 

n low-concentration air 14, temporal concentration of the 

q compound in the air is computed; whereas temporal residual 

iji 

15 amounts of the compound on the floor 5 is computed by use 

^ of Fugacity (f 5 ) of the floor 5. 

At S5, it is judged whether or not to perform safety 
evaluation in the case where the chemical is inhaled. If 
the safety evaluation in inhalation is to be performed, then 

20 an estimated exposure amount in inhalation indicating a 

degree of exposure upon inhalation of the contaminated air 
is computed by use of the above-mentioned temporal 
concentration in the indoor air (S6). Thereafter, an 
inhalation safety coefficient is computed according to the 

25 estimated exposure amount in inhalation (S7). At S8, the 

inhalation safety coefficient is compared with a reference 
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value defined in each country. If the inhalation safety 
coefficient exceeds the reference value, it is judged that 
"there is no problem in safety," By contrast, if the 
inhalation safety coefficient is lower than the reference 
value at S8, it is judged that "there is a problem in safety, " 
and the operation returns to SI, where alteration of the 
primary condition such as alteration of compound, 
alteration of chemical formulation, alternation of using 
yQ condition, or the like is considered. 

O 10 If the safety evaluation in inhalation is not selected 

ry at S5, calculated is an estimated amount of percutaneous 

5 - ! 

m exposure indicating to what extent the skin is exposed in 



13 contact with the floor 5 to which the chemical is attached 

q (S9). Thereafter, at S10, it is judged whether to perform 

[fj 15 a percutaneous safety evaluation or not. If the 

m percutaneous safety evaluation is to be performed, a 

percutaneous safety coefficient is computed according to 
the estimated percutaneous exposure amount (Sll). At S12, 
as with the safety evaluation in inhalation, the 
20 percutaneous safety coefficient is compared with a 

reference value defined in each country, whereby the safety 
is evaluated. If there is a problem in safety at S12, the 
operation returns to SI, where alteration of the primary 
condition is considered. 
25 If the percutaneous safety evaluation is not selected 

at S10, calculated at S13 according to the estimated 
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percutaneous exposure amount is an estimated oral exposure 
amount indicating the degree of exposure in the case where 
the chemical attached to a hand or the like is taken orally 

(513) . Subsequently, according to the estimated oral 
exposure amount, an oral safety coefficient is computed 

(514) . This oral exposure may occur, in particular, when 
an infant puts a chemical-attached hand into the mouth. At 
S15, as with the safety evaluation in inhalation, the oral 
safety coefficient is compared with a reference value 
defined in each country, whereby the safety is evaluated. 
If there is a problem in safety at S15, the operation returns 
to SI, where alteration of the primary condition is 
considered . 

Finally, if it is judged to be safe at S8, S12, or S15, 
overall safety is evaluated. Here, the sum of respective 
reciprocals of the previously determined inhalation safety 
coefficient, percutaneous safety coefficient, and oral 
safety coefficient is determined, and the reciprocal of thus 
determined value is defined as an overall safety coefficient 
(S16). This overall safety coefficient is compared with a 
reference value, whereby an evaluation similar to the 
previous safety evaluation is effected (S17). 

In the following, the above-mentioned steps of S2 to 
S4, S6, S7, S9, Sll, S13, and S14 will be explained in detail . 

(i) Secondary Condition Calculation (S2) 

Calculations of the definite factor condition shown 
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in Table 2 9 and the temporal change condition shown in Table 
3 0 will be explained. 

The condensed particles 9, 10, and 11 are successively 
generated during time t x and are absorbed by the high- 
concentration ceiling 15. Specifically, the condensed 
particles 9 exist during the period of time ti/3 after 
evaporation, the condensed particles 10 are generated after 
the lapse of time t x /3, and the condensed particles 11 are 
generated after the lapse of time 2t x /3. Here, the floating 
number (n if i = 9, 10, 11) of the condensed particles is 
constant from their generation till absorption, and is 
expressed, from the specific gravity ( p ) of the compound, 
compound content (R a ), emission rate (E T ), condensation 
ratio (C r ), and initial diameter (d 0 ) of condensed particles 
9, 10, and 11, as follows: 



wherein G is the ventilation rate. 

According to surface temperature (T d ) of particles, 
properties of particle-constituting solvent (vapor 
pressure P d , molecular weight M d , and specific gravity p 
d ), and indoor environment (air diffusion coefficient D air , 
vapor pressure poo, and room temperature T°o), the 
evaporation constant ( a ) of the condensed particles 9, 10, 
and 11 is defined as follows: 




...(54) 
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a 



_4D air M d fP, Poo\ 



...(55) 



Q 

in io 



m 

;T1' 



wherein R is a gas constant. 

Here, the upper parts of the above-mentioned Tables 
2 7 and 2 8 exemplify the case where the heat-vaporized aerosol 
is of oil base. In this case, Tco (temperature at the site 
far from particles) and T d are set to room temperature (T), 

whereas P°o (vapor pressure of oil at the site far from 
particles) is set to zero. 

In the case where the heat-vaporized aerosol is of 
water base, it is necessary for P d , M d , p d , T d , and poo in 
expression (55) to be changed to values based on properties 
of water. Specifically, the above-mentioned factor values 
when the room temperature (T) is 298 [K] (= 25°C) and the 
relative humidity is 60[%RH] are exemplified in the lower 
parts of Tables 27 and 28. 

Here, T d can be determined, according to the room 
temperature (T) and relative humidity (H), from "mass-based 
humidity table" disclosed in "Kagaku Kikai no Riron to Keisan 
(Theory and Calculation of Chemical Machines)" (Second 
Edition) (Saburo Kamei ed., Sangyo Tosho) or the like. On 
the other hand, P d and poo can be computed from the following 
expressions : 

logio^d = 10.23 - 1750/(T d -38) 

log 10 poo = 8.23 - 1750/(T-38) + log lo 0 
wherein xfj is an indoor relative humidity. 
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From the molecular weight (M a ) of the compound, 
emission rate (E T ), and condensation rate (C r ), the amount 
of emission of chemical (E v ) to the high-concentration air 
12 is expressed as follows: 



The surface area (A k/ k = 5, 6, 15, 16, 17) of each 
of the floor 5, wall 6, high-concentration ceiling 15, 
medium-concentration ceiling 16, and low-concentration 
ceiling 17 is determined by the room size (V 4 ) and floor area 
(L 4 x W 4 ). 

Using the room size (V 4 ), ventilation rate (G), and 
cross-sectional area (A 14 ) of the room with respect to the 
direction of movement of indoor air, the indoor air velocity 
(v f ) is determined as follows: 



From the vapor pressure (P 3 ) and melting point (T M ) 
of the compound and surface temperature (T d ) of particles, 
Fugacity capacity {2> if i = 9, 10, 11) of the condensed 
particles 9, 10, and 11 is expressed as: 



E V =E T (1-C r )/M a 



...(56) 



V # 




...(57) 



z , 6xl ° 6 

' P L RT d 

PI =P s exp{6.79(T M /T d -l)} 



...(58) 



Though the vapor pressure of the compound in a liquid state 
(P L 3 ) is computed by use of T M and T d here; in the case where 
T M cannot be obtained, P L 3 may be set identical to the vapor 
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pressure of the compound in a solid state (P s ) . 

From the room temperature (T) , Fugacity capacity (Zj, 
j = 12, 13, 14) of the airs 12, 13, and 14 is expressed as 
follows : 



RT 

From the vapor pressure (P s ), water solubility (C s ), 
and octanol/water distribution coefficient (K ow ) of the 
compound, and from the oil component content ( p k , k = 5, 
6, 15, 16, 17) of the material constituting the floor 5, 
wall 6, and ceilings 15, 16, and 17, Fugacity capacity (Z k , 
k = 5, 6, 15, 16, 17) of the floor 5, wall 6, and ceilings 
15, 16, and 17 is expressed as: 
Z* = P k K ow C>lY» ...(60) 

Mainly generated in the 11 kinds of media are 
degradation reactions caused by light and oxidation. The 
degradation constant (K if i = 9, 10, 11) of the condensed 
particles 9, 10, and 11, degradation constant (Kj, j = 12, 
13, 14) of the airs 12, 13, and 14, and degradation constant 
(K k , k = 5 , 6, 15, 16, 17) of the floor 5 , wall 6 , and ceilings 
15, 16, and 17 are defined by use of the respective half-lives 
of degradation by light /oxidation ( T L ) , (T ^) , and (T k ). 
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0.693 
0.693 



0.693 



.(61) 



ri 



10 



15 



20 



WhileTi andr k are often unavailable, K L and K k may be 
nullified in this case from the viewpoint of securing safety 

of inhabitants. WhenTj of the compound in the air is 
unavailable, it may be determined by use of Atmospheric 
Oxidation program (Atkinson et al., 1984, Chem. Rev. Vol. 
84, pp. 437-470) from the structural formula of the compound . 

Since the volume of particles decreases over time due 
to evaporation of the particle-constituting solvent 
component, during the evaporation of particle-constituting 
solvent component, using the diameter (d 0 ) of particles 
immediately after spraying (t = 0) and the evaporation 
constant ( a ) determined from expression (55) , the diameter 
(d L , i = 9, 10, 11) of the condensed particles 9, 10, and 
11 at time t is represented as: 

4 = V<C - 2ar (62) 

From the diameter (d £ ) of condensed particles 9, 10, 
and 11 determined by expression (62), indoor velocity (v f ) 
determined by expression (57), and Fugacity capacity (Z i ) 
of condensed particles 9, 10, and 11 determined by expression 
(58), and Fugacity capacity (Z n ) of the medium- 
concentration air 13 determined by expression (59), the 
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transference coefficient (D il3 ) between the condensed 
particles 9, 10, and 11 and the medium-concentration air 
13 becomes: 

D il3 = — (63) 

l/(M,-Z ; ) + l/(* l3 4Z,3) 

wherein A L (= ftA*) is the surface area of particles, k 13 ( = 
v f + v i) is the transference speed of the compound in the 
medium-concentration air 13, and k L (= v f /100) is the 
transference speed of the compound in the particles. Here, 
v £ is the floating speed of the condensed particles 9, 10, 
and 11, and can be determined by use of hydrodynamics since 
it is under the control of the flow caused by the heat of 
the electric heating vaporizer, gravity, and air 
resistance. 

From the half-life of transference to air ( T jk ) in the 
ceilings 15, 16, and 17, volume (Vj, j = 12, 13, 14) of the 
airs 12, 13, and 14, volume (V k ) of the floor /wall/ceiling 
determined by expression (66) mentioned later, Fugacity 
capacity (Zj) of airs 12, 13, and 14 determined by expression 
(59), and Fugacity capacity (Z k ) of floor /wall/ceiling 
determined by expression ( 60 ) , the transference coefficient 
(D jk ) of the compound between the medium-concentration air 
13 and the ceilings 15 and 16 or between the low- 
concentration air 14 and the floor 5, wall 6, and low- 
concentration ceiling 17 becomes: 
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0.693 



(64) 



T y4 {L/(V.Z,.) + l/(^Z,)j 




1 



(64') 



Here, the total volume V 12 + V 13 + V 14 of the airs 12, 13, and 
14 equals the volume (V 4 ) of the room. 

The above-mentioned transference coefficient (D jk ) 
may also be determined by the expression represented by the 
above-mentioned expression (64 ' ) by use of the transference 
speed (kj) of the compound in the air and the transference 
speed (k k ) of the compound in the interior materials (floor 
5, wall 6, and ceilings 15, 16, and 17). 

Here, the half-life of transference to air (T jk ) in 
expression (64) is hard to actually measure and varies under 
various conditions, thus making it difficult to determine 
an accurate value thereof. By contrast, all the parameters 
in expression (64 ' ) can be easily determined by calculation 
without necessitating measurement. For example, using the 
diffusion coefficient (D c ) of compound, the transference 
speed (k k ) of the compound in the interior material is 
calculated as (D c /t) 0,5 ; whereas A jk can be determined from 
the contact area between the individual airs 12, 13, and 
14 and the individual interior materials (floor 5, wall 6, 
and ceilings 15, 16, and 17). Here, D c can be determined 
when the diffusion coefficient (D ca ) of the compound in the 
air is multiplied by 10" 9 , whereas D ca can be computed from 
properties of the compound (e.g., structural formula, 
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molecular weight, and the like) according to Wike and Lee 
method ( "Handbook of Chemical Property Estimation Methods, " 
McGraw-Hill Book Company , 1982), for example. Accordingly, 
the transference coefficient (D jk ) is more preferably 
5 determined by expression (64') than by expression (64). 

Using the diameter (d t ) of particles determined by 
expression (62) and the evaporation constant (or), change 
in volume (dVi/dt) of the condensed particles 9, 10, and 11 

==4 

=S is expressed by: 

H 10 ■3T--T arf i ( 65 > 

III ^ 

3x! 

ijl Namely, as the particle-constituting solvent component 

p evaporates over time, the volume (V L , i = 9, 10, 11) of 

I 3 

g condensed particles 9, 10, and 11 decreases. 

Assuming that each of the floor 5, wall 6, and ceilings 
15 15, 16, and 17 before heat-vaporizing the chemical is like 

a thin film and that, as the chemical infiltrates into the 
film, its thickness increases so as to enhance the volume 
thereof; from the surface area (A k ) of the 
floor/wall/ceiling and diffusion coefficient (D c ) of the 
20 compound in the floor/wall/ceiling, the volume (V k , k = 5, 

6, 15, 16, 17) of floor/wall/ceiling becomes: 
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V k = 2VD7A k (66) 

(ii) Fugacity Calculation (S3) 

The behavior of the compound in the condensed 
25 particles 9, 10, and 11 is expressed in the form of 
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differential equation concerning Fugacity (fi) as: 



— V-7>i = 
dt ' ' 



...(67) 



F - change V - change Transference(13) Degradation 

Here, the terms of V-change, Transference (13), and 
Degradation respectively indicate the volume change 
(decrease over time) of the condensed particles, amount of 
transference of the compound between the condensed 
particles 9, 10, and 11 and the medium-concentration air 
13, and amount of degradation of the compound by light and 
oxidation in the condensed particles 9, 10, and 11 (KiViZi 
indicating the degradation speed quantity) . Here, when the 
condensed particles 9, 10, and 11 are absorbed by the 
high-concentration ceiling 15, Fugacity (f L ) immediately 
disappears and is added to Fugacity (f 15 ) explained later. 

The behavior of the compound in the high-concentration 
air 12 is expressed in the form of differential equation 
concerning Fugacity (f 12 ) as: 



Here, the terms of Emission, Transference ( j ) , and 
Degradation respectively indicate the amount of emission 
of the compound as complete vapor, amount of transference 
of the compound with respect to the medium-concentration 
air 13, and amount of degradation of the compound by light 





Degradation 



...(68) 
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and oxidation (K 12 V 12 Z 12 indicating the degradation speed 
quantity) . 

The behavior of the compound in the medium- 
concentration air 13 is expressed in the form of differential 
equation concerning Fugacity (f 13 ) as: 

^ 3 Zi3 - Vffafnfn "(4a + Ai)Z l3 f l3 + 4 3 Z 14 / 14 } -(69) 
F - change Transference^]) 

- 2 "P«Afn -/,)- I A 3 *(/i3-A) 

x-9 jk-15 

Transference^]) Transference(k) 

~ ^13^13^13/13 
Degradation 

Here, the terms of Transference ( j ) , Transference (i), 
Transference (k) , and Degradation respectively indicate the 
amount of transference of the compound with respect to the 
high-concentration air 12 and low-concentration air 14, 
amount of transference of the compound with respect to the 
condensed particles 9, 10, and 11, amount of transference 
of the compound with respect to the ceilings 15 and 16, and 
amount of degradation of the compound by light and oxidation 
(K 13 V l3 Z 13 indicating the degradation speed quantity) . 

The behavior of the compound in the low-concentration 
air 14 is expressed in the form of differential equation 
concerning Fugacity (f 14 ) as: 
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= v f {A l3 Z u f u - (A l3 + A l4 )Z 14 /„ } 



(70) 



F - change 



Transference^ j) & Ven tilation 



6 



*-5 



Transfer en ce (k ) Transfer en ce(17) 



^14^14^14/14 
Degradation 



Here, the terms of Transference (j) & Ventilation, 
Transference (k), Transference (17), and Degradation 
respectively indicate the amount of transference of the 
compound with respect to the medium-concentration air 13 
and the amount of discharge thereof outdoors, amount of 
transference of the compound with respect to the floor 5 
and wall 6, amount of transference of the compound with 
respect to the low-concentration ceiling 17, and amount of 
degradation of the compound by light and oxidation (K 14 V 14 Z 14 
indicating the degradation speed quantity) . 

The behaviors of the compound in the floor 5 and wall 
6 are respectively expressed in the forms of differential 
equations concerning Fugacity (f 5 ) and (f 6 ) as follows: 



••■(71) 



F - change V - change Transference(\A) 




• 
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-^V 6 Z 6 = -4Djt\ZJ 6 -D U6 {f 6 - f u ) 



•••(72) 



F - change V - change Transference^^) 





Here, in each expression, the terms of V-change, 
Transference (14), and Degradation respectively indicate 
the volume change (increase over time) of the floor 5 or 
wall 6, amount of transference of the compound with respect 
to the floor 5 or wall 6, and amount of degradation of the 
compound by light and oxidation (K 5 V 5 Z 5 or K 6 V 6 Z 6 indicating 
the degradation speed quantity) . 

The behaviors of the compound in the high- 
concentration ceiling 15, medium-concentration ceiling 16, 
and low-concentration ceiling 17 are respectively expressed 
in the forms of differential equations concerning Fugacity 
(f 15 ), (f 16 ), and (f 17 ) as follows: 



r-KsZ* « -V^AsZls/lS ~ A3 l5 (/t5 - /l3) 



...(73) 



F - change V - change Transference(13) 
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df i6 



V l6 Z l6 = -jDJtAnZufu-Du 16 (/ 16 - / l3 ) -(74) 



dt 

F - change V - change Transference^) 
Degradation 



dt 

F - change V - change Transference^.^) 



dt V l7 Z„ = -4DjtA, n Z„f Y1 -D u 17 (/ 17 - f l4 ) -(75) 



?y - K l7 V l7 Z l7 f l7 

:IT Degradation 

fii Here, in each expression, the terms of V-change, 

ITS Transference (13), Transference (14), and Degradation 

S3 5 respectively indicate the volume change (increase over 

q time) of the ceilings 15, 16, and 17, amount of transference 

of the compound between the high-concentration ceiling 15 
or medium-concentration ceiling 16 and the medium- 
concentration air 13 , amount of transference of the compound 
10 between the low-concentration ceiling 17 and the low- 

concentration air 14, and amount of degradation of the 
compound by light and oxidation (K 15 V 15 Z 15 , K 16 V 16 Z 16 , or K 17 V 17 Z 17 
indicating the degradation speed quantity). 

When the following Fugacity (f a ) is added to the floor 
15 (k = 5), wall (k = 6), and ceiling (k = 15, 16, 17) by a 

predetermined ratio, accuracy can be improved. 
f a -ny iia Z.J iB IV k IZ k •■•(76) 

Here, V tia and f ia respectively indicate the volume and 
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Fugacity of particles immediately before the condensed 
particles attach to the interior materials. 

The above-mentioned 11 kinds of differential 
equations (67) to (75) are simultaneously formed and are 
5 solved by Runge-Kutta-Gill method, so as to compute Fugacity 

(f s/ f 6/ f 9 to f 17 ) . 

When solving these simultaneous differential 
equations, it is necessary to set an estimation nick time 
width (dt) which is a minute time unit. Namely, the 
10 estimation nick time width is used such that solutions of 

the simultaneous differential equations are initially 
J~ determined at a time (t 0 ), and then solutions of the 

5^ simultaneous differential equations are determined at a 

,~; time (t 0 + dt) to which the estimation nick time width is 

15 added. As solutions are obtained while estimation nick time 

widths are successively added, temporally changing Fugacity 
can be determined. Theoretically, as the set time of the 
estimation nick time width is shorter, more accurate 
solutions can be obtained, though necessitating a very long 
20 calculation time. By contrast, when the set time is too long, 

solutions tend to diverge, thereby generating errors. 

Therefore, in the present invention, the estimation 
nick time width is set shorter when a very large change occurs 
in a chemical, whereas it is set longer when there is no 
25 large change. 

Specifically, mass balance is always confirmed such 
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that the amount of input of the chemical and the resulting 
solution coincide with each other, and the estimation nick 
time width is set longer when the mass balance does not 
fluctuate greatly, whereas it is set shorter when the mass 
balance starts fluctuating. For example, when the 
fluctuation of mass balance is set to an accuracy of ±5%, 
the estimation nick time width is always set so as to 
constantly satisfy the relationship of: 

compound input amount/ ( existing amount + degrading 
amount + discharging amount) = 0,95 to 1.05 
namely, such that the fluctuation of mass balance lies within 
the range of ±5%. 

Here, the above-mentioned compound input amount is 
determined by the compound content (R a ) and emission rate 
( E T ) . Since the temporal amounts of compound in the 11 kinds 
of media are determined by the simultaneous differential 
equations, they are summed up so as to compute the existing 
amount as shown in the following: 



n 6 t7 

existing amount = l J nV^J i + 2 f k v k Z k + 1 fk V k^k 



<-9 Jfc-5 £-12 

It 6 17 



degrading amount = Zn.KtfZJ, + 2 K k f k V k Z t + 2 KJ k V k Z k 

i-9 Jfc-5 Jfc-12 



discharging amount - v / A l4 Z 14 / u 



•(77) 



The degrading amount and discharging amount are as 
represented above . 

With reference to the flowchart of Fig. 5, the method 
of setting the estimation nick time width (dt) will be 
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explained . 

First, an initial value of estimation nick time width 
( dt ) is inputted ( S2 1 ) ♦ Then , an upper limit set value (e.g. , 
0.1[%]) which is the upper limit of difference in mass 
balance, and a lower limit set value (e.g., 10~ 6 [%]) which 
is the lower limit of difference in mass balance are inputted 
(S22) . Thereafter, Fugacity and mass balance at t = t 0 are 
calculated (S23 and S24), and Fugacity and mass balance at 
t = t + dt (or t 0 + dt) are calculated (S25 and S26). 

It is judged whether the fluctuation in mass balance 
is within the range of ±5% or not (S27) . If the fluctuation 
in mass balance is within the range of ±5%, it is judged 
whether the difference between the mass balance at t = t 
(or t 0 ) and the mass balance at t = t + dt (or t 0 + dt) is 
at least the upper limit set value or not (S28). If it is 
judged to be at least the upper limit set value at S28, then 
solutions become more accurate when the estimation nick time 
width (dt) is made shorter. In this case, the estimation 
nick time width (dt) is multiplied by 1/2 so as to change 
its setting (S29). When the difference is judged to be 
smaller than the upper limit set value at S28, it is judged 
at S30 whether the difference between the mass balance at 
t = t (or t 0 ) and the mass balance at t = t + dt (or t 0 + 
dt) is at most the lower limit set value or not. 

When the mass balance difference is not greater than 
the lower limit set value at S30, since solutions are not 
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influenced by longer estimation nick time width (dt), the 
estimation nick time width (dt) is doubled so as to change 
its setting (S31). Subsequently, it is judged whether the 
estimation nick time width (dt) changed at S31 is at most 
a maximum value (e.g., 0.1 [hour]) of estimation nick time 
width (dt) or not (S32). When the estimation nick time width 
(dt) is not greater than the maximum value at S3 2, since 
solutions do not diverge, the estimation nick time width 
(dt) set at S31 is used. When the estimation nick time width 
(dt) is greater than the maximum value at S32, since 
solutions may diverge, the estimation nick time width (dt) 
is reset to the maximum value (S33) . When the mass balance 
difference is greater than the lower limit set value at the 
above-mentioned S3 0, namely, when it lies between the lower 
limit set value and upper limit set value, the calculations 
are continued without changing the estimation nick time 
width (dt). 

After the step of S29, S30, S32, or S33, the operation 
returns to S25 so as to effect calculation again, and this 
process is repeated till the aimed time is attained. 

When the fluctuation in mass balance exceeds the range 
of ±5% at S27 , on the other hand, the mass balance fluctuation 
is so much that calculation is preferably effected with an 
estimation nick time width (dt) shorter than that in the 
case where the fluctuation is within the range of ±5%. 
Accordingly, the calculation is stopped once (S34), the 
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lower limit set value is reset to a lower level (S35), and 
then the operation returns to the step of S23. 

Thus, when the estimation nick time width is variably 
set without being held constant, while monitoring the mass 
balance fluctuation, Fugacity can be computed accurately 
and efficiently. 

(iii) Computation of Temporal Concentration in Indoor Air 
and Floor Residual Amount (S4) 

The temporal concentration of the compound in the 
g 10 indoor air is computed when Fugacity (f 14 ) of the low- 

ly concentration air 14, determined by the above-mentioned 

jjj item(ii) , is multiplied by its Fugacity capacity (Z 14 ) . The 

« concentration in the low-concentration air 14 is used 

q because the most part in the room is filled with the 

15 low-concentration air 14, and humans are most likely to 

inhale the low-concentration air 14. The concentration in 
the air of other concentrations may also be used when 
necessary . 

The floor residual amount of the compound is computed 
20 when Fugacity (f 5 ) of the floor is multiplied by Fugacity 

capacity ( Z 5 ) . 

(iv) Calculation of Estimated Exposure Amount in Inhalation 
and Inhalation Safety Coefficient (S6 and SI) 

The above-mentioned temporal concentration in the 
2 5 indoor air forms a curve shown in Fig. 3A, for example. This 

concentration curve is integrated, an accumulated 



a e * 
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concentration of the compound during a specific period (t x 
to t 2 ) is determined (see Fig, 3B), and the mean concentration 
in the indoor air is computed from thus determined value. 
While an arbitrary period is set as the specific period 
5 depending on the object, an appropriate period is usually 

set in view of the method of use of the product and the test 
- period of toxicity data. 

Then, according to the above-mentioned mean 

Q 

m concentration in the indoor air, amount of respiration, and 

Q 10 exposure time, the estimated exposure amount in inhalation 

nj is determined. Namely, calculation of: 

If! estimated exposure amount in inhalation [mg/kg/day] 

□ = mean concentration in indoor air [mg/m 3 ] x amount of 

H respiration [m 3 /kg/min] x exposure time [min/day] 

.h 15 is effected. Here, as the above-mentioned amount of 

respiration, a published value or actually measured value 
may be used. Also, when amounts of respiration are 
respectively set for adult and child, more appropriate 
estimated exposure amounts in inhalation can be obtained. 
20 In the case where the inhaled harmful material is not totally 

absorbed into the body but is partially discharged by 
respiration, a more appropriate estimated exposure amount 
in inhalation can be obtained when the inhalation ratio is 
taken into account. 
25 The inhalation safety coefficient is computed from a 

non-influential amount concerning inhalation toxicity 
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examined by an animal experiment beforehand and the 
estimated exposure amount in inhalation determined above. 
Namely, it is expressed as: 

inhalation safety coefficient = 
5 inhalation non-influential amount [mg/kg/day]/ 

estimated exposure amount in inhalation [mg/kg/day] 
(v) Calculation of Estimated Percutaneous Exposure Amount 
and Percutaneous Safety Coefficient (S9 and Sll) 
S The above-mentioned floor residual amount forms a 

□ 10 curve such as that shown in Fig, 4A, for example. This 

y residual amount curve is integrated, the accumulated 

H residual amount of the compound during a specific period 

3 (tt to t 2 ) is determined (see Fig. 4B) , and the mean residual 

=5 amount is computed from thus determined value. While an 

=r? 

f* 15 arbitrary period is set as the specific period depending 

on the object, an appropriate period is usually set in view 
of the method of use of the product and the test period of 
toxicity data . 

Then, according to the mean floor residual amount, 
2 0 skin attachment ratio, contact area, and body weight, the 

estimated percutaneous exposure amount is determined. 
Namely, calculation of: 

estimated percutaneous exposure amount [mg/kg/day] = 
(mean residual amount [mg/m 2 ] * skin attachment ratio [%] 
25 x contact area [mVday ] ) /body weight [kg] 

is effected. Here, as the contact area, a published value 
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(e.g., 4 [mVday]) may be used. The skin attachment ratio 
is a ratio of the compound attaching to the skin when the 
latter is in contact with the floor 5 where the compound 
exists. As this value, a published value or a value 
experimentally obtained from a model may be used. 

A model experiment method for the skin attachment 
ratio is as follows. A weight (8 cm x 8 cm * 8 cm; 4.2 kg) 
is placed on a denim cloth ( 8 cm x io cm) with a pressure 
similar to that of an infant in contact with a floor, and 
the denim cloth is pulled on the floor at a speed (120 cm/15 
sec) similar to the moving speed of the infant. The denim 
and floor are analyzed so as to compute the compound 
contained in the denim and floor. From the ratio 
therebetween, the skin attachment ratio is obtained. It has 
been confirmed that the skin attachment ratio obtained by 
this method is identical to or slightly higher than that 
determined from analyzed values of a hand and a floor when 
the hand is actually pressed against the floor, thereby 
proving this model experiment method to be useful for 
evaluating exposure of inhabitants. 

The percutaneous safety coefficient is computed from 
the non-influential amount concerning percutaneous 
toxicity examined by an animal experiment beforehand and 
the estimated percutaneous exposure amount determined above. 
Namely, it is expressed as: 

percutaneous safety coefficient = 

184 



SCC 97-06 



percutaneous non-influential amount [mg/kg/day ] / 
estimated percutaneous exposure amount [mg/kg/day] 

Nevertheless, in general, percutaneous non- 
influential amount has not often been determined, and there 
are not many published values . Accordingly, a more accurate 
value can be determined from the estimated percutaneous 
exposure amount, and oral non-influential amount and 
percutaneous absorption ratio for which many published 
values exist, according to the following expression: 

percutaneous safety coefficient = oral non- 
influential amount [mg/kg/day] / (estimated percutaneous 
exposure amount [mg/kg/day] * percutaneous absorption ratio 

[%] ) 

Here, when the percutaneous absorption ratio is unknown, 
employed is a national guideline (e.g., 10%) which usually 
exists . 

(vi) Calculation of Estimated Oral Exposure Amount and Oral 
Safety Coefficient (S13 and S14) 

From the estimated percutaneous exposure amount 
obtained in the above-mentioned item (v) , hand surface area 
ratio, and oral transference ratio, the estimated oral 
exposure amount from hand to mouth is determined. Namely, 
calculation of: 

estimated oral exposure amount [mg/kg/day] = 
estimated percutaneous exposure amount [mg/kg/day] * hand 
surface area ratio [%] x oral transference ratio [%] 
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is effected. Here, the hand surface area ratio is expressed 
by (hand surface area/body surface area), for which a 
published value (e.g., 5[%]) may be used. The oral 
transference ratio is a hypothetical value, which is set 
to 100%, for example. 

In the case where oral exposure might occur via 
tableware or food contaminated with the residually sprayed 
compound, it is required that the estimated oral exposure 
amount from tableware or food to mouth be added to the 
estimated oral exposure amount from hand to mouth to yield 
the total estimated oral exposure amount. For example, the 
estimated oral exposure amount from tableware is obtained 
when, according to the tableware residual amount indicating 
the amount of the harmful material remaining in tableware, 
tableware use area which is the sum of tableware surface 
areas, and oral transference ratio from tableware, 
calculation of: 

estimated oral exposure amount [mg/kg/day] = 
tableware residual amount [mg/m 2 ] x tableware use area 
[mVday] x oral transference ratio [%]/body weight [kg] 
is effected. Here, the tableware residual amount is 
expressed by (mean floor residual amount x tableware 
contamination ratio) . As the tableware contamination ratio, 
an actually measured value (e.g., 9% ) or a hypothetical value 
may be used. 

The oral safety coefficient is computed from the 
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non-influential amount concerning oral toxicity examined 
by an animal experiment beforehand and the estimated oral 
exposure amount determined above. Namely, it is expressed 
as : 

oral safety coefficient = 

oral non-influential amount [mg/kg/day]/ estimated 
oral exposure amount [mg/kg/day] 

As mentioned in the foregoing, in the method of 
estimating an indoor behavior of a pesticidal compound in 
this embodiment, the environment is divided into the 
condensed particles 9, 10, and 11 that are classified into 
three kinds depending on occurrence and extinction, the 
high-concentration air 12, medium-concentration air 13, and 
low-concentration air 14 that are classified into three 
kinds depending on the compound concentration, the floor 
5, the wall 6, and the high-concentration ceiling 15, 
medium-concentration ceiling 16, and low-concentration 
ceiling 17 that are classified into three kinds depending 
on the compound concentration; and differential equations 
concerning Fugacity in these members are simultaneously 
formed and are solved, whereby the indoor behavior of the 
compound when heat-vaporizing a chemical containing the 
compound is estimated. Here, when solving the simultaneous 
differential equations, the estimation nick time width is 
variably set, while constantly confirming mass balance of 
the compound indoors after the spraying, so that the amount 
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of input of the chemical indoors and the resulting solution 
coincide with each other. 

Accordingly, since mass balance of the compound after 
the spraying is always monitored such that the amount of 
input of the compound indoors and the resulting solution 
coincide with each other, thereby varying the estimation 
nick time width; the estimation nick time width is set longer 
when the mass balance fluctuates a little, whereas it is 
set shorter when the mass balance starts fluctuating greatly . 
Namely, when solving simultaneous differential equations 
including a parameter accompanying temporal change, the 
estimation nick time width is automatically set in response 
to the fluctuation in mass balance. Consequently, when 
processed by a computer, an accurate solution can be obtained 
in a short time. 

The method of evaluating safety of a pesticidal 
compound in accordance with this embodiment uses the 
estimated result mentioned above to evaluate the safety of 
the compound with respect to the human body when the chemical 
is heated to vaporize. 

Accordingly, the safety of the compound with respect 
to the human body can be evaluated accurately in a short 
time. As a consequence, when formulating a chemical such 
as insecticide containing the compound, simulation can be 
easily repeated while changing conditions, thereby making 
it easier to formulate a chemical having a high safety 
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conforming to the aimed object. 

Though the kind of the floor 5 is assumed to be tatami 
in this embodiment, differential equations similar to those 
mentioned above may also be formed in the case of flooring, 
5 In this case, however, when determining fugacity capacity 

(Z 5 ) of the floor 5 in equation (60), it is necessary for 
the particle-constituting solvent component content (see 
Table 27) to be changed from that of tatami to that of 

set 

q flooring. 

3 10 Though Fugacity is determined by use of Runge- 

(j Kutta-Gill method in this embodiment, other methods may be 

S used for solving differential equations. Runge-Kutta-Gill 

method, however, is preferably used since a program for the 
£ above-mentioned differential equations can be easily made 

3? 

; I 

^ 15 by Basic . Also in the case where differential equations are 

* solved by a method other than Runge-Kutta-Gill method, 

similar effects can be obtained when the estimation nick 

time width is set as mentioned above. 

2 0 (6) Sixth Embodiment 

With reference to the drawings, the sixth embodiment 

of the present invention will be explained in the following. 

The method of simulating an indoor behavior of a pesticidal 

compound in accordance with this embodiment estimates the 
25 behavior of the compound in the case where the floor in the 

fifth embodiment is constituted by a rug having ears of fiber . 
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For convenience of explanation, members identical to those 
shown in the drawings of the previous embodiment will be 
referred to with numerals or letters identical thereto, 
without their explanations provided. Explained here is a 
5 case where the kind of floor is a carpet. Accordingly, only 

the portion different from the fifth embodiment will be 
explained . 

First, the carpet is modeled as shown in Figs. 13A, 
13B, and 13C. Namely, the carpet is divided into a plurality 

10 of carpet fibers (ears of fiber) planted on a substrate and 

space portions existing between the fibers . When a chemical 
is heated to vaporize in a room in which such a carpet covers 
the whole floor surface, as condensed particles fall on the 
carpet, the chemical infiltrates into the carpet fiber and 

15 accumulates in the space portion. 

While the method of simulating an indoor behavior of 
a pesticidal compound in accordance with this embodiment 
also conforms to the flowchart of Fig. 1, at SI, as the 
primary condition, the carpet-related primary condition 

20 shown in Table 31 is added to the primary condition of the 

fifth embodiment (see Tables 25 to 28). Here, the 
cross-sectional area (A c ) of the space portion of the carpet 
refers to the area in which the space portion is in contact 
with the indoor air (see Figs. 13A, 13B, and 13C). 
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The examples of values in Table 31 are based on the 
values shown in Figs. 13A, 13B, and 13C. Here, one space 
portion is set to a size of 0.1 [mm] * 0.1 [mm] * 3 [mm], 
and one square of carpet fiber is set to a size of 2.5 [mm] 
x 2.5 [mm] x 3 [mm]. Also, it is assumed that 16 [pieces] 
of carpet fiber exist in 1 [cm 2 ]. 

In this case, the individual parameters can be 
determined as : 

A c = (0.1 [mm] x o.l [mm] x 16 [pieces])/l [cm] 
= 1.6 x io -3 [m 2 /m 2 floor] 

A c5 = (0.1 [mm] x 3 [mm] x 4 [faces] x 16 [pieces] )/l 

[cm] 

= 0.192 [m 2 /m 2 floor] 
A 5 = 1 - A c 

= 1 - 1.6 x 10" 3 [m 2 /m 2 floor] 

^a5 = ^5 ~*~ ^c5 

= 1.19 [m 2 /m 2 floor] 
Here, A 5 is the area in which the carpet fibers in 1 [m 2 ] 
are in contact with the indoor air. 

As the secondary condition at S2, the carpet-related 
secondary condition shown in Table 32 is added to the 
secondary condition of the fifth embodiment. Here, k 5 and 
V 5 , which are respectively transference speed of the 
compound in the floor and the volume of the floor in the 
fifth embodiment (see Table 30), are respectively replaced 
by the transference speed of the compound in the carpet fiber 
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and the volume of the whole fiber of the carpet. 
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A method of computing the definite factor conditions 
and temporal change conditions shown in the above-mentioned 
Table 32 will be explained hereinafter as exemplified by 
the case using an aqueous aerosol, i.e., case where the 
particle-constituting solvent is water. 

From the cross-sectional area (A c ) of the space portion, 
room temperature (T), and humidity (H), the evaporation 
constant (R^) of the chemical entered into the space portion 
of the carpet is expressed as: 

R d -Aji(H m -H)/C„ .-.(78) 

C H = 0.24 + 0.46// 

A = 3.06xl0" 4 (r-rj l/3 

wherein T d is the temperature of the space portion of the 
carpet (corresponding to the wet-bulb temperature), and H m 
is the humidity (corresponding to the saturated humidity 
at the wet-bulb temperature) of the space portion at T d . 
Here, T d and H m can be determined from the above-mentioned 
"mass-based humidity table" according to the room 
temperature (T) and humidity (H). In the case of an 
oil-based product (aerosol), while the vapor constant of 
the chemical can be computed when properties of the water 
and oil are compared with each other and then are multiplied 
by R,,, Rj = 0 may be assumed in the case where the evaporation 
of the oil is slow. 

From the diffusion coefficient (D c ) and the area (A a5 ) 
in which the carpet fiber is in contact with the compound, 
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the volume (V s ) of the whole fiber of the carpet is determined 
as : 

V 5 = 2jDjA a5 -(79) 

The transference speed (k c ) of the compound in the 
5 product is expressed as k 14 /100 by use of k 4 (= Gv 4 /A 14 ) which 

indicates the transference speed of the compound in the 
indoor air. 

From the above-mentioned transference speed (k c ) of 
the compound in the product, transference speed (k 14 ) of the 

10 compound in the indoor air, floor area (L 4 x w 4 ) , cross- 

sectional area (A c ) of the space portion, Fugacity capacity 
(Z c ) of the chemical in the space portion determined by 
expression (82) which will be mentioned later, and Fugacity 
capacity (Z 14 ) of the low-concentration air determined by 

15 expression (59) of the fifth embodiment, the transference 

coefficient (D cl4 ) of the compound between the chemical in 
the space portion of the carpet and the low-concentration 
air is expressed as : 

D cl4 - — (80) 

l/(kAMKZJ^l/(k i4 A c L 4 W 4 Z l4 ) 

2 0 From the above-mentioned transference speed (k c ) of 

the compound in the product, diffusion coefficient (D c ), 
floor area (L 4 x W 4 ), area (A c5 ) in which the carpet fiber 
is in contact with the compound in the space portion, 
Fugacity capacity (Z c ) of the chemical in the space portion, 

25 and Fugacity capacity (Z 5 ) of the floor determined by 
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expression (60) of the fifth embodiment, the transference 
coefficient (D c5 ) of the compound between the chemical in 
the space portion and the carpet fiber is expressed as: 

D = 1 ...(8i) 

1/(M. 5 W.) + 1/(M.5W5) 

wherein k 5 (= (D c /t) 0 " 5 ) is the transference speed of the 
compound within the carpet fiber. 

The volume (V c ) of the chemical in the space portion 
is determined from the floor area (L 4 * W 4 ) , emission rate 
(E T ), condensation rate (C r ), and evaporation constant (Rj) 
of the chemical in the space portion determined by expression 
(78). The volume (V c ) of the space portion decreases until 
the solvent in the chemical completely evaporates. 

From the vapor pressure (P s ), water solubility (C 3 ), 
and octanol/water distribution coefficient (K ow ) of the 
compound, and property change constant (a) of the product, 
Fugacity capacity (Z c ) of the chemical in the space portion 
is expressed as follows: 
Z c =(e- a 'K ow +l-e- at )C s / P s -..(82) 

Here, the property change constant (a) of the product 
can be computed according to the evaporation constant (R^) 
of the chemical in the space portion determined by expression 
(78), content (C 3d ) of the solvent (Xylene) contained in the 
product, and dilution (X sol ) of the product. Namely, 
Fugacity capacity (Z c ) of the chemical in the space portion 
is also expressed as: 
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Z c = [ (ratio of organic solvent in product) K ow + (ratio 
of water in product)] x c s /P° 

By simultaneously forming this expression and the 
above-mentioned expression (82), the property change 
constant (a) of the product is computed. Here, in the case 
of the oil-based product (aerosol), Fugacity capacity (Z c ) 
of the chemical in the space portion is given by expression 
(57) of the fifth embodiment from the vapor pressure (P 3 ), 
water solubility (C 3 ), and octanol/water distribution 
coefficient (K ow ) of the compound. 

Consequently, the behavior of the compound in the 
space portion of the carpet is expressed in the form of 
differential equation concerning Fugacity (f c ) as follows: 

- Y c ^nA^ZJ, Deposition^) -(83) 

+ R d L 4 W 4 ZJ c V -change 

F - change + D cj (f c - /. ) Transference^]) 

+ D *i3 (fc - A3 ) Transfernce(k) 
- KV c Z c f c Degradation 

Also, expressions (71) and (70) in the fifth 
embodiment are respectively replaced by the following 
expressions (84) and (85): 
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V - change 



•••(84) 



11 



+ V n } V. Z . /. Deposition (/) 



+ £> c5 (f 5 ~f c ) Transference^) 
+ A 4 5 (/s " /i4 ) Transfernce(14) 
- K 5 V 5 Z 5 f 5 Degradation 




GV u Z u f u 



Ventilation 



-(85) 



u 



J«.-A»(/i4-/-) 

2A4/(/ M -/,) 



Transference^) 



Transference^ j) 




Degradation 
Transferen ce(c) 



Namely, as shown in expression (84) , the differential 
equation concerning the floor (carpet fiber) includes, in 
addition to expression (71) of the fifth embodiment, the 
terms of Deposition (i) indicating the amount of deposition 
of condensed particles onto the carpet fiber and 
Transference (c) indicating the amount of transference of 
the compound between the carpet fiber and the space portion 

In expression (85), which is differential equation 
concerning the low-concentration air, terms of Ventilation, 
Transference (i), Transference (j), Degradation, and 
Transference (c) respectively indicate the amount of 
discharge of the compound outdoors, amount of transference 
of the compound between the low-concentration air and 
condensed particles, amount of transference of the compound 
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between the low-concentration air and the medium- 
concentration air and floor, amount of degradation of the 
compound by light and oxidation, and amount of transference 
of the compound between the low-concentration air and the 
space portion. 

Accordingly, Fugacity calculation at S3 is performed 
such that 12 kinds of differential equations of expressions 
(67) to (69) and (72) to (75) in the fifth embodiment and 
the above-mentioned expressions (83) to (85) are solved by 
Runge-Kutta-Gill method. The subsequent steps of S4 to S17 
are similar to those in the fifth embodiment. 

As explained in the foregoing, in the method of 
estimating an indoor behavior of a pesticidal compound in 
accordance with this embodiment, in the case where the kind 
of floor is a carpet having ears of fiber and a space portion 
between the ears, a differential equation for Fugacity in 
the above-mentioned space portion of the carpet is added 
to 11 kinds of simultaneous differential equations 
concerning Fugacity in 3 kinds of condensed particles, air 
divided into 3 kinds depending on the compound concentration, 
floor, wall, and ceiling divided into 3 kinds depending on 
the compound concentration, so as to calculate Fugacity. 

Consequently, even in the case where the chemical is 
accumulated in the space portion, so that the compound 
behaves differently from the case with a flat floor, the 
behavior of the compound can be estimated accurately, 
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whereby various kinds of simulations can be performed. 

( 7 ) Seventh Embodiment 

With reference to the drawings, the seventh embodiment 
of the present invention will be explained in the following. 

The method of simulating an indoor behavior of a 
pesticidal compound (hereinafter simply referred to as 
compound) contained in a chemical (solution) such as 
insecticide in accordance with this embodiment estimates 
the behavior of the compound in the case where the chemical 
is sprayed over the whole floor surface, and is applied to 
_the case where a carpet (rug) made of polymer fiber is spread 
on the whole floor surface. Here, usable as the chemical 
are pyrethroid insecticidal compounds , organophosphorus 
compounds, carbamate compounds, insect growth restrainers 
(IGR), and the like. Its solvent may be water or oil. 
Explained in this embodiment is the case where the solvent 
is water. The vermin to be killed is mainly flea, mite, and 
the like. 

First, the carpet is modeled as shown in Figs. 13A, 
13B, and 13C. Namely, the carpet is divided into a plurality 
of carpet fibers (ears of fiber) planted on a substrate and 
space portions existing between the fibers. When a chemical 
is sprayed over the whole floor area covered with a carpet, 
since the solvent of the chemical is water, the chemical 
infiltrates into the carpet fiber and accumulates in the 
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space portion. The chemical not only attaches to the floor 
but also is partly suspended in the air. 

Accordingly, as shown in Fig. 17, the indoor 
environment covered with the carpet is classified into seven 
kinds of media comprising the space portion of the carpet, 
two kinds of suspended large and small particles having 
different diameters, whole indoor air, floor (carpet fiber) , 
wall, and ceiling. 

The suspended particles are divided into two kinds 
since the behavior of the particles mainly depend on their 
diameter. Namely, while the individual particles are 
distributed in the same place immediately after spraying, 
as time passes, the large particles are distributed to a 
particle zone DZ 2 close to the floor due to their weight, 
and the small particles are distributed to a particle zone 
DZ 3 close to the ceiling since they are lighter than the large 
particles . 

Simultaneous differential equations concerning 
Fugacity of the compound in the above-mentioned seven kinds 
of media are formed and are solved by Runge-Kutta-Gill method, 
whereby the indoor behavior of the compound in the chemical 
when the chemical is sprayed over the whole carpet surface, 
i.e., temporal change in the compound distribution in the 
room and amount of attachment to the interior materials of 
the floor, wall, ceiling, and the like, is estimated and 
analyzed . 
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With reference to Fig. 1, a method of simulating an 
indoor behavior of the compound will be explained. The 
simulation method in accordance with this embodiment can 
be mainly divided into a step of dividing an indoor 
environment into predetermined media and forming a 
differential equation concerning a fugacity of the compound 
in each medium (SI to S2) ; a step of determining the fugacity 
of the compound in each medium from the differential equation 

(53) ; a step of determining the indoor behavior of the 
compound from the fugacity of the compound in each medium 

(54) ; a step of changing, in response to a fluctuation in 
mass balance of the compound indoors, a minute time unit 
used when solving the differential equation (S21 to S35, 
see Fig. 5) ; and a step of evaluating safety of the compound 
with respect to a human body according to the indoor behavior 
of the compound (S5 to S17). 

First, at step (hereinafter abridged as S) 1, a primary 
condition is inputted. The primary condition comprises 
physicochemical properties of the compound (see Table 33), 
indoor environmental behavior properties of the compound 
(see Table 34), indoor environment (see Table 35), a carpet 
condition (see Table 36) , and product properties (see Table 
37) . 
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The examples of values in Table 3 5 assume those of a 
six-mat room (9.72 [m 2 ]) of a typical apartment in Japan, 

The cross-sectional area (A c ) of the space portion of 
the carpet in Table 3 6 is an area in which the space portion 
is in contact with the indoor air (see Figs. 13A, 13B, and 
13C) . Further, the examples of values in Table 3 6 are based 
on the values shown in Figs. 13A, 13B, and 13C. Here, one 
space portion is set to a size of 0.1 [mm] x 0.1 [mm] x 3 
[mm] , and one square of carpet fiber is set to a size of 
2.5 [mm] x 2.5 [mm] x 3 [mm]. Also, it is assumed that 16 
[pieces] of carpet fiber exist in 1 [cm 2 ]. 

In this case, the individual parameters can be 
determined as : 

A c .= (0.1 [mm] x o.l [mm] x 16 [pieces])/l [cm] 
= 1.6 x 10* 3 [mVm 2 floor] 

A c5 = (0.1 [mm] x 3 [mm] x 4 [faces] x 16 [pieces])/l 

[cm] 

= 0.192 [m 2 /m 2 floor] 
A 5 = 1 - A, 

= 1 - 1.6 x 10" 3 [m 2 /m 2 floor] 

^a5 = A 5 + A 15 

= 1.19 [m 2 /m 2 floor] 
Here, A 5 is the area in which the carpet fibers in 1 [m 2 ] 
are in contact with the indoor air. 

Subsequently, using the primary condition, a 
secondary condition is determined by calculation (S2 ) . The 
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secondary condition comprises a temporally unchangeable 
definite factor condition (see Table 38) determined by the 
primary condition alone, and a temporal change condition 
accompanying a temporal change (see Table 39). Their 
calculations will be explained later in detail. 
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TABLE 38 



DEFINITE FACTOR CONDITION 


SYMBOL 


EVAPORATION CONSTANT OF SUSPENDED PARTICLES 
VOLUME RATIO OF COMPOUND IN PRODUCT 
DIFFUSION rOFFFIPIFNT OF rOMPOTTND TN ATR 
EVAPORATION CONSTANT OF CHEMICAL IN SPACE PORTION 


a 

Ra 

Rd 


FUGACITY CAPACITY 

OF SUSPENDED PARTICLE 
AIR 

FLOOR / WALL / CEILING 


Zj (i=2,3) 
Z 4 

Zj 0=5,6,7) 


DEGRADATION CONSTANT 

OF CARPET SPACE PORTION 
SUSPENDED PARTICLE 
AIR 

FLOOR / WALL / CEILING 


K c 

Kj (i=2,3) 
K4 

Kj 0=5,6,7) 



V 
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TABLE 39 



TEMPORAL CHANGE CONDITION 


SYMBOL 


SUSPENDED PARTICLE 
DlAM.bl.bR 
FALLING SPEED 
SUSPENDED NUMBER 


di (i-2,3) 
vi (i=2,3) 
ni (i=2,3) 


TRANSFERENCE SPEED OF COMPOUND 

BETWEEN CHEMICAL IN SPACE PORTION AND AIR 
BETWEEN CHEMICAL IN SPACE PORTION AND CARPET FIBER 
BETWEEN SUSPENDED PARTICLE AND AIR 
BETWEEN FLOOR / WALL / CEILING AND AIR 


D C 4 
D C 5 

Di4 (i=2,3) 
D 4 j 0=5,6,7) 


VOLUME OF CHEMICAL IN SPACE PORTION 
VOLUME OF SUSPENDED PARTICLE 
VOLUME OF FLOOR / WALL / CEILING 


Dc 

Di (i=2,3) 
Dj 0=5,6,7) 


FUGACITY CAPACITY OF CHEMICAL IN SPACE PORTION 


Z c 
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Using the secondary condition, seven kinds of Fugacity 
are calculated (S3). Namely, seven kinds of differential 
equations concerning the space portion of the carpet, two 
kinds of large and small suspended particles, indoor air, 
floor (carpet fiber), wall, and ceiling are simultaneously 
formed and are solved by Runge-Kutta-Gill method, whereby 
seven kinds of Fugacity are computed over time. Here, an 
estimation nick time width (minute time unit) set when 
solving the differential equations is automatically set so 
as to be varied in response to a fluctuation in mass balance 
of the compound. 

At S4, using thus computed Fugacity of the small 
particles (taking account of Fugacity of the large particles 
when necessary) and Fugacity of the indoor air, temporal 
concentration of the compound in the air is computed; whereas 
floor residual amount of the compound is computed by use 
of Fugacity of the carpet fiber and carpet space portion. 

At S5, it is judged whether or not to perform safety 
evaluation in the case where the chemical is inhaled. If 
the safety evaluation in inhalation is to be performed, then 
an estimated exposure amount in inhalation indicating a 
degree of exposure upon inhalation of the contaminated air 
is computed by use of the above-mentioned temporal 
concentration in the indoor air (S6). Thereafter, an 
inhalation safety coefficient is computed according to the 
estimated exposure amount in inhalation (S7). At S8, the 
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inhalation safety coefficient is compared with a reference 
value defined in each country. If the inhalation safety 
coefficient exceeds the reference value, it is judged that 
"there is no problem in safety." By contrast, if the 
5 inhalation safety coefficient is lower than the reference 

value at S8, it is judged that "there is a problem in safety, " 
and the operation returns to SI, where alteration of the 
primary condition such as alteration of compound, 
alteration of chemical formulation, alternation of using 

10 condition, or the like is considered. 

If the safety evaluation in inhalation is not selected 
at S5, calculated is an estimated amount of percutaneous 
exposure indicating to what extent the skin is exposed in 
contact with the floor to which the chemical is attached 

15 (S9). Thereafter, at S10, it is judged whether to perform 

a percutaneous safety evaluation or not. If the 
percutaneous safety evaluation is to be performed, a 
percutaneous safety coefficient is computed according to 
the estimated percutaneous exposure amount (Sll). At S12, 

20 as with the safety evaluation in inhalation, the 

percutaneous safety coefficient is compared with a 
reference value defined in each country, whereby the safety 
is evaluated. If there is a problem in safety at S12, the 
operation returns to SI, where alteration of the primary 

25 condition is considered. 

If the percutaneous safety evaluation is not selected 
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at S10, calculated at S13 according to the estimated 
percutaneous exposure amount is an estimated oral exposure 
amount indicating the degree of exposure in the case where 
the chemical attached to a hand or the like is taken orally 

(513) . Subsequently, according to the estimated oral 
exposure amount, an oral safety coefficient is computed 

(514) . This oral exposure may occur, in particular, when 
an infant puts a chemical-attached hand into the mouth. At 
S15, as with the safety evaluation in inhalation, the oral 
safety coefficient is compared with a reference value 
defined in each country, whereby the safety is evaluated. 
If there is a problem in safety at S15, the operation returns 
to si, where alteration of the primary condition is 
considered. 

Finally, if it is judged to be safe at S8, S12, or S15, 
overall safety is evaluated. Here, the sum of respective 
reciprocals of the previously determined inhalation safety 
coefficient, percutaneous safety coefficient, and oral 
safety coefficient is determined, and the reciprocal of thus 
determined value is defined as an overall safety coefficient 
(S16). This overall safety coefficient is compared with a 
reference value, whereby an evaluation similar to the 
previous safety evaluation is effected (S17). 

In the following, the above-mentioned steps of S2 to 
S4 , S6 , S7 , S9 , Sll , S13 , and S14 will be explained in detail . 

(i) Secondary Condition Calculation (S2) 
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Calculations of the definite factor condition shown 
in Table 3 8 and the temporal change condition shown in Table 
3 9 will be explained. 

According to surface temperature (T d ) of particles, 
properties of particle-constituting solvent (vapor 
pressure P d/ molecular weight and specific gravity p 

d ) , and indoor environment (air diffusion coefficient D air/ 
vapor pressure poo, and room temperature Too), the 
evaporation constant ( a ) of the suspended particles is 
defined as follows: 



wherein R is a gas constant, and P d (= 2, 26 x 10 3 [Pa]), M d 
( = 18 [g/mole]), and p d (= 1 * 10 6 [g/m 3 ]) are assumed to 
be the values due to properties of water . Also , in this case , 
Too (temperature at the site far from particles) is set to 
room temperature (T), T d to indoor wet-bulb temperature (= 
292.5 [K] ) , and poo (vapor pressure of water at the site far 
from particles) to 1 .89 * 10 3 [Pa]. These values are those 
with room temperature (T) of 298 [K] and relative humidity 
of 60[%RH] . 

Here, T d can be determined, according to the room 
temperature (T) and relative humidity (H), from "mass-based 
humidity table" disclosed in "Kagaku Kikai no Riron to Keisan 
(Theory and Calculation of Chemical Machines)" (Second 
Edition) (Saburo Kamei ed., Sangyo Tosho) or the like. On 
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3=J 



the other hand, P d and P°o can be computed from the following 

expressions : 

log 10 P d = 10.23 - 1750/(T d -38) 
lpg 10 poo = 8.23 - 1750/(T-38) + log lo 0 
5 wherein is an indoor relative humidity. 

From the cross-sectional area (A c ) of the space portion, 

room temperature (T), and humidity (H), the evaporation 

constant (R^) of the chemical entered into the space portion 

of the carpet is expressed as: 

R d =A c h(H m -H)IC H ...(87) 
10 C H =0.24 + 0.46# 

h = 3.06xlO~\T-T d ) l/3 

wherein H m is a humidity at temperature (T d ) . The chemical 
entered into the space portion of the carpet reduces its 
volume upon evaporation and finally disappears, leaving the 
compound. Thus left compound infiltrates through the 
15 carpet fiber. 

From the surface temperature (T d ) of particles and the 
vapor pressure (P 3 ) of the compound, Fugacity capacity (Z if 
i = 2, 3) of the large and small particles is expressed as: 

2,-^ ...<<*) 

Pi =p'ex P {6.79(r„/:r,-i)} 

20 Though the vapor pressure of the compound in a liquid state 

(P L S ) is computed by use of the melting point (T M ) of the 
compound and T d here; in the case where T M cannot be obtained, 
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P L S may be set identical to the vapor pressure of the compound 
in a solid state (P 8 ) . 

From the room temperature (T) , Fugacity capacity (Z 4 ) 
of the air is expressed as follows: 

Since the carpet fiber and wallpaper covering the 
floor /wall/ceiling are formed by polymers; from the vapor 
*Q pressure (P s ), water solubility (C s ), and octanol/water 

0 distribution coefficient (K ow ) of the compound, and from the 

rd 10 oil component content ( Pj) of the carpet fiber and wallpaper 

1 - * 

|fi covering the floor/wall/ceiling, Fugacity capacity (Z 5 , j 

g = 5, 6, 7) of the floor/wall/ceiling is expressed as: 

q Z. = Pj K ow C s /P s ...(90) 

3 J 5 

,n The degradation constants (K c ), (K L , i = 2, 3), (K 4 ), 

15 (Kj, j = 5, 6, 7) of the carpet space portion, suspended 

particles, air, and floor /wall/ceiling are defined by use 
of their respective half-lives of degradation by 
light/oxidation (T c ), (T L ), (T 4 ), and (T^ as follows: 



K c = 0.693/r 
K. = 0.693/t,. 
K, = 0.693/t 4 
Kj = 0.693/r. 



...(91) 



2 0 Mainly generated in the seven kinds of media are 

degradation reactions caused by light and oxidation, and 
the reaction speed quantity at this time is represented by 
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the product of degradation constant, volume, and Fugacity 
capacity. Among the degradation constants shown in 
expression (91), the degradation constant (K 4 ) of the 
compound in the air is hard to actually measure, and it is 
not usually easy to determine this value. Accordingly, when 
the actually measured value is unavailable, it may be 
determined as computed by Atmospheric Oxidation program 
(Atkinson et al., 1984, Chem. Rev. Vol. 84, pp. 437-470) 
from the structural formula of the compound. Also, from the 
viewpoint of securing safety of inhabitants, the 
degradation constants in the particles and in the interior 
materials (floor/wall/ceiling) may be nullified when their 
actually measured values are unavailable. 

Since the volume of particles decreases over time due 
to evaporation, using the evaporation constant (oc) 
determined from expression (86) and the diameter (d oi ) of 
suspended particle immediately after spraying (t = 0), the 
diameter (d if i = 2, 3) of the suspended particles at time 
t is determined as: 

d. = V4-2ar ...(92) 
d; = \[R.d„ -(92 1 ) 

On the other hand, when the particle-constituting 
solvent in the suspended particle is completely evaporated, 
the diameter (d L ) becomes that of the compound alone, thereby 
being held constant. Namely, the particle diameter (di) 
after the complete evaporation of the particle-constituting 
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solvent is represented by the above-mentioned expression 
(92* ) by use of the volume ratio (R a ) of the compound in the 
product and the diameter (d oi ) of the particles immediately 
after spraying. 

From the compound content (C a ) and the specif ic gravity 
( p ) of the compound, the volume ratio (R a ) of the compound 
in the product is expressed as C a / p . 

As shown in Fig. 18, immediately after spraying the 
chemical, a particle zone (DZ ± ) having a width (W 4 ), length 
(L 4 ), and space width (H (i) ) is formed at a certain height 
from the floor. In conformity to the falling speed (Vi) of 
suspended particles, the particle zone (DZi) approaches the 
floor. When the bottom of the particle zone (DZi) reaches 
the floor, the compound in the suspended particles transfers 
to the carpet space portion and carpet fiber. 

Since movement of the particles is under the control 
of gravity and air resistance, according to Stokes ' law using 
the diameter (d ± ) of suspended particles determined by 
expression (92), the falling speed (v if i = 2, 3) of the 
suspended particles is represented by the following 
expression ( 93 ) : 



I8r/ 



£-d?=P{dl t -2ca) ...(93) 




[6.32 + 2.01exp(-8.322xl0 6 4)] ...(93') 



wherein g is gravitational acceleration, rj is the 
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coefficient of viscosity of the air, S c is a sliding 
correction coefficient , and/? is a speed coefficient . Here, 
while particles do not conform to the Stokes 1 law when they 
become small, the sliding correction coefficient (S c ) , which 
is a coefficient for correcting this phenomenon, is 
represented by the above expression (93* ) according to the 
diameter (d^ of particles determined by expression (92). 

On the other hand, the falling speed (\r L ) of the 
particles after the complete evaporation of particle- 
constituting solvent component is represented by the above 
expression (93") using the specific gravity (p) of the 
compound, specific gravity ( p d ) of the particle- 
constituting solvent, speed coefficient ( /3 ) , and diameter 
(di) of particles determined by expression (92*). 

The number of suspended particles (n if i = 2, 3) also 
decreases upon exchange of the air through a window and the 
like. Here, the small particles float in conformity to the 
indoor air flow, and the large particles, which reduce their 
diameter as their volume decreases, are assumed to conform 
to the flow of the indoor air flow. In this case, during 
the period till the particle zone (DZi) reaches the floor, 
the number of particles (n^ would decrease only upon 
ventilation of the room. Accordingly, from the room size 
( V 4 = L 4 x W 4 x H 4 ) , ventilation rate ( G) , particle space width 
(H (i) , i = 2, 3), and particle falling speed (vj determined 
by expression (93) or (93 H ), the particle number (ni) is 
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expressed as follows: 

*b..-Jh-Vi-!Lk 4 ...(94) 
dt H (i) L> 4 

wherein k 4 {= GV 4 /A 4 ) is the transference speed of the 
compound in the air. 

From the floor size (L 4 * W 4 ), ventilation rate (G), 
cross-sectional area (A c ) of the space portion, Fugacity 
capacity (Z c ) of the chemical in the space portion determined 
by expression (101) mentioned later, and Fugacity capacity 
(Z 4 ) of the air determined by expression (89), the 
transference coefficient (D c4 ) of the compound between the 
chemical in the space portion and the air is represented 
as follows : 

D c4 ...(95) 

l/{k^ 4 W 4 Z c ) + l/(k 4 AJL 4 W 4 Z 4 ) 

wherein k c ( = k 4 /100) is the transference speed of the 
compound on the water surface of the chemical in the space 
portion. 

From the diffusion coefficient (D p ) , floor size (L 4 
x W 4 ), ventilation rate (G), Fugacity capacity (Z c ) of the 
chemical in the space portion, and Fugacity capacity (Z 5 ) 
of the floor determined by expression (90) , the transference 
coefficient (D c5 ) of the compound between the chemical in 
the space portion and the carpet fiber is represented as: 
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wherein k 5 (= (D p /t) 05 ) is the transference speed of the 
compound in the carpet fiber, and is the contact area 
ratio of the carpet fiber and space portion relative to the 
floor, 

5 The diffusion coefficient (D p ) of the compound in the 

floor/wall/ceiling can be determined when the diffusion 
coefficient (D ca ) of the compound in the air is multiplied 
by 10~ 9 , whereas D ca can be computed from properties of the 
=o compound (e.g., structural formula, molecular weight, and 

£3 10 the like) according to Wike and Lee method ("Handbook of 

nj Chemical Property Estimation Methods," McGraw-Hill Book 

S - ! 

ijl Company, 1982), for example. 

S3 From the diameter (d ± ) of particles determined by 

r=3 expression (92) or (92'), transference speed (k 4 ) of the 

jJT 

*g 15 compound in the air, falling speed (v A ) determined by 

^ expression (93) or (93"), Fugacity capacity (ZJ and (Z 4 ) 

determined by expressions (88) and (89), the transference 
coefficient (D i4 ) of the compound between the suspended 
particles and the air becomes: 

20 D i4 « ...(97) 

l/(^;) + l/{(^ + ^)A^)} 

wherein Ai ( = Kd* ) is the surface area of suspended particles , 

and k £ (= k 4 /100) is the transference speed of the compound 

in the particle surface. 

From the diffusion coefficient (D p ) of the 

25 floor/wall/ceiling, room size (L 4 x w 4 x z 4 ) , ventilation rate 
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(G), and Fugacity capacity (Z 4 ) and (Zj) determined by 
expressions (89) and (90), the transference coefficient 
(D 4j ) of the compound between the air and the 
floor/wall/ceiling becomes: 

5 A -(98) 

4 ' 1/(* 4 A/Z4)+1/(*A Z /) 

wherein is the surface area of the wall and ceiling, and 

kj (= (D p /t)°* s ) is the transference speed of the compound 

,S in the floor/wall/ceiling. 

q The volume (V c ) of the chemical in the space portion 

hi 10 is determined from the floor size (L 4 * W 4 ) , spray rate, 

5 . 3 

j«j content (C 3d ) of the solvent (Xylene) contained in the 

^ product, dilution (X gol ) of the product, and evaporation 

,=2 constant (R d ) of the chemical in the space portion determined 

by expression (87). The volume (V c ) of the space portion 
*** 15 decreases until the particle-constituting solvent in the 

chemical completely evaporates . 

Since the suspended particles are mostly made of 
moisture (particle-constituting solvent), their volume 
decreases as time passes. Since the suspended particles 
20 reduce their volume till the particle-constituting solvent 

completely evaporates, using diameter (di) of suspended 
particles determined by expression (92) or (92') and the 
above-mentioned evaporation constant ( a ) , the volume 
change (dv i /dt) of suspended particles is expressed as 
25 follows: 
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§5f 



dV. 71 

— = ad. ...(99) 

dt 2 ' V } 

Assuming that the floor /wall/ceiling before spraying 
the chemical is like a thin film and that, as the chemical 
infiltrates into the film, its thickness increases so as 
5 to enhance the volume thereof; from the diffusion 

coefficient (D p ) of the compound in the floor/wall/ceiling, 
the volume (Vj) of floor/wall/ceiling is represented as 
yfj follows: 
5 V. = 2{D p tf 5 A } ...(100) 

hi 10 Fugacity capacity (Z c ) of the chemical in the space 

portion is initially expressed by K ow C s /P s as with the 
above-mentioned Fugacity of the floor/wall/ceiling. When 
the chemical begins evaporating, Fugacity capacity (Z c ) 
gradually approaches C 3 /P s (Fugacity capacity of water). 
15 Accordingly, from the vapor pressure (P s ) , water solubility 

(C 3 ), and octanol/water distribution coefficient (K ow ) of the 
compound, and property change constant (a) of the product, 
Fugacity capacity (Z c ) is expressed as follows: 
Z c = (e~ at K ow +l-e~ at )C s /P s ...(101) 
2 0 Here, the property change constant (a) of the product 

can be computed according to the evaporation constant (R^) 
of the chemical in the space portion determined by expression 
(87), content (C 3d ) of the solvent (Xylene) contained in the 
product, and dilution (X aol ) of the product • Namely, 
25 Fugacity capacity (Z c ) of the chemical in the space portion 
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is also expressed as: 

Z c = [ (ratio of organic solvent in diluted liquid) K ow 
+ (ratio of water in diluted liquid)] x c s /P s 
By simultaneously forming this expression and the 
5 above-mentioned expression (101), the property change 

constant (a) of the product is computed, 
(ii) Fugacity Calculation (S3) 

The behavior of the compound in the carpet space 
=|g portion is expressed in the form of differential equation 

q 10 concerning Fugacity (f c ) as: 

| rfT KZ < = y^ n - v ' A ^ Z ^ /H C-) Deposition^) ...(102) 

F - change 

Sj + R d L 4 W 4 Z C f c V - change 

O 5 

iS - ^ D ck (f c - f k ) Transference(4,5) 

m -KV c ZJ c Degradation 

Here, the terms of Deposition (i), V-change, Transference 
(4,5), and Degradation respectively indicate the attachment 
accompanying the falling of the suspended particles to the 

15 space portion, volume change (decrease over time) of the 

chemical in the space portion, amount of transference of 
the compound between the space portion and the air and floor 
(carpet fiber), and amount of degradation of the compound 
in the space portion. Here, Fugacity (f c ) is effective till 

20 the chemical in the space portion completely evaporates. 

Also, 7 5 indicates the ratio of suspended particles 
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attaching to the floor. 

The behavior of the compound in the large and small 
particles is expressed in the form of differential equation 
concerning Fugacity (f if i = 2, 3) as: 

&v?, = jad?,/, V- change ...(103) 

F - change 

5 - DM - / 4 ) Transference^) 

y -K.V.Z./. Degradation 

p Here, the terms of V-change, Transference (4), and 

hj Degradation respectively indicate the volume change 

: . : 

ij\ (decrease over time) of the suspended particles, amount of 

!=! transference of the compound between the particles and the 

10 air, and amount of degradation of the compound in the 

^ particles. Here, with the lapse of time, particles are 

absorbed by the carpet space portion and floor ( carpet fiber ) , 
whereby Fugacity (f L ) is lost. 

The behavior of the compound in the indoor air is . 
15 expressed in the form of differential equation concerning 

Fugacity ( f 4 ) as: 



dt 4 4 



= -GV 4 Z A f 4 Ventilation ...(104) 

F - change -T> e4 (f 4 - f e ) Transference^) 

- 2 n i°i* (/•-/) Transference^) 



- jA,,(/ 4 -/,) Transference^}) 

- K 4 V £ /, Degradation 
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Here, the terms of Ventilation, Transference (c), 
Transference (i), Transference (j), and Degradation 
respectively indicate the amount of discharge of the 
compound outdoors, amount of transference of the compound 
between the air and the space portion, amount of transference 
of the compound between the air and the suspended particles, 
amount of transference of the compound between the air and 
the floor/wall/ceiling, and amount of degradation of the 
compound in the air. 

The behavior of the compound in the floor (carpet 
fiber) is expressed in the form of differential equation 
concerning Fugacity (f 5 ) as: 

^y 5 Zs = Ys^n^A^ZJJH^ Deposition(i) ...(105) 
F - change 

~^D p ltA a5 Z 5 f 5 V- change 

" D cs (fs - fc ) Transference(c) 

- A* 5 (/ 5 " / 4 ) Transference^) 

- K s V 5 Z 5 f 5 Degradation 

Here, the terms of Deposition (i), V-change, Transference 
(c), Transference (4), and Degradation respectively 
indicate the attachment of the suspended particles 
accompanying falling to the carpet fiber, volume change 
(increase over time) of the floor, amount of transference 
of the compound between the carpet fiber and the space 
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portion, amount of transference of the compound between the 
carpet fiber and the air, and amount of degradation of the 
compound in the carpet fiber. 

The behavior of the compound in the wall and ceiling 
is expressed in the form of differential equation concerning 
Fugacity (fj, j = 6, 7) as: 

^jf-VjZ, = - ^TtAjZjfj V - change ...(106) 

F - change 

3 

+ r j ^ nyytifi l H {i) Deposition 

-D 4j ( f. - f 4 ) Transference 

- K.VjZj fj Degradation 
Here, the terms of V-change, Transference (4), and 
Degradation respectively indicate the volume change 
(increase over time) of the wall/ceiling, amount of 
transference of the compound between the wall/ceiling and 
the air, and amount of degradation of the compound in the 
wall/ceiling. Also, 7 j is the ratio of suspended particles 
attaching to the wall and ceiling. 

The above-mentioned seven kinds of differential 
equations (102) to (106) are simultaneously formed and are 
solved by Runge-Kutta-Gill method, so as to compute Fugacity 
( f 2 to f 7 , and f c ) • 

When solving these simultaneous differential 
equations, it is necessary to set an estimation nick time 
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width (dt) which is a minute time unit. Namely, the 
estimation nick time width is used such that solutions of 
the simultaneous differential equations are initially 
determined at a time (t 0 ), and then solutions of the 
simultaneous differential equations are determined at a 
time (t 0 + dt) to which the estimation nick time width is 
added . As solutions are obtained while estimation nick time 
widths are successively added, temporally changing 
Fugacity can be determined. Theoretically, as the set time 
of the estimation nick time width is shorter, more accurate 
solutions can be obtained, though necessitating a very long 
calculation time. By contrast, when the set time is too long, 
solutions tend to diverge, thereby generating errors. 

Therefore, in the present invention, the estimation 
nick time width is set shorter when a very large change occurs 
in a chemical, whereas it is set longer when there is no 
large change. 

Specifically, mass balance is always confirmed such 
that the amount of input of the chemical and the resulting 
solution coincide with each other, and the estimation nick 
time width is set longer when the mass balance does not 
fluctuate greatly, whereas it is set shorter when the mass 
balance starts fluctuating. For example, when the 
fluctuation of mass balance is set to an accuracy of ±5%, 
the estimation nick time width is always set so as to 
constantly satisfy the relationship of: 
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compound input amount/ ( existing amount + degrading 
amount + discharging amount) = 0,95 to 1.05 
namely, such that the fluctuation of mass balance lies within 
the range of ±5%. 

Here, the above-mentioned compound input amount is 
determined by the floor size (L 4 x W 4 ) , compound content (C a ) , 
and spray rate. Since the temporal amounts of compound in 
the seven kinds of media are determined by the simultaneous 
differential equations, they are summed up so as to compute 
the existing amount as shown in the following. Also, the 
degrading amount and discharging amount are as follows: 



existing amount = V c ZJ c + Y fyPx + Y fk V k Z k 

degrading amount = K c V c Z c f c + yK i n.f.V.Z, + jK k f k V k Z k --(107) 



With reference to the flowchart of Fig. 5, the method 
of setting the estimation nick time width (dt) will be 
explained , 

First, an initial value of estimation nick time width 
(dt ) is inputted (S21 ) . Then, an upper limit set value (e.g. , 
0.1[%]) which is the upper limit of difference in mass 
balance, and a lower limit set value (e.g., 10" 6 [%]) which 
is the lower limit of difference in mass balance are inputted 
(S22) . Thereafter, Fugacity and mass balance at t = t 0 are 
calculated (S23 and S24), and Fugacity and mass balance at 
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t = t + dt (or t 0 + dt) are calculated (S25 and S26). 

It is judged whether the fluctuation in mass balance 
is within the range of ±5% or not (S27) . If the fluctuation 
in mass balance is within the range of ±5%, it is judged 
whether the difference between the mass balance at t = t 
(or t 0 ) and the mass balance at t = t + dt (or t 0 + dt) is 
at least the upper limit set value or not (S28). If it is 
judged to be at least the upper limit set value at S2 8, then 
solutions become more accurate when the estimation nick time 
width (dt) is made shorter. In this case, the estimation 
nick time width (dt) is multiplied by 1/2 so as to change 
its setting (S29). When the difference is judged to be 
smaller than the upper limit set value at S28, it is judged 
at S3 0 whether the difference between the mass balance at 
t = t (or t 0 ) and the mass balance at t = t + dt (or t 0 + 
dt) is at most the lower limit set value or not. 

When the mass balance difference is not greater than 
the lower limit set value at S30, since solutions are not 
influenced by longer estimation nick time width (dt), the 
estimation nick time width (dt) is doubled so as to change 
its setting (S31). Subsequently, it is judged whether the 
estimation nick time width (dt) changed at S31 is at most 
a maximum value (e.g., 0.1 [hour]) of estimation nick time 
width (dt) or not (S32). When the estimation nick time width 
(dt) is not greater than the maximum value at S32, since 
solutions do not diverge, the estimation nick time width 
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(dt) set at S31 is used. When the estimation nick time width 
(dt) is greater than the maximum value at S32, since 
solutions may diverge, the estimation nick time width (dt) 
is reset to the maximum value (S33) . When the mass balance 
5 difference is greater than the lower limit set value at the 

above-mentioned S3 0, namely, when it lies between the lower 
limit set value and upper limit set value, the calculations 
are continued without changing the estimation nick time 
5 width (dt) . 

q 10 After the step of S29, S30, S32, or S33, the operation 

IU returns to S25 so as to effect calculation again, and this 

: S process is repeated till the aimed time is attained, 

<«* When the fluctuation in mass balance exceeds the range 

™ of ±5% at S27 , on the other hand, the mass balance fluctuation 

15 is so much that calculation is preferably effected with an 

estimation nick time width (dt) shorter than that in the 
case where the fluctuation is within the range of ±5%. 
Accordingly, the calculation is stopped once (S34), the 
lower limit set value is reset to a lower level (S35), and 
20 then the operation returns to the step of S23. 

Thus, when the estimation nick time width is variably 
set without being held constant, while monitoring the mass 
balance fluctuation, Fugacity can be computed accurately 
and efficiently . 

2 5 (iii) Computation of Temporal Concentration in Indoor Air 

and Floor Residual Amount (S4) 
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The temporal concentration of the compound in the 
indoor air is computed when Fugacity (f 3 ) of the small 
particles, determined by the above-mentioned item (ii), 
multiplied by Fugacity capacity (Z 3 ), and Fugacity (f 4 ) of 
the indoor air multiplied by Fugacity capacity (Z 4 ) are 
summed up. Here, in the case where the large particles may 
be inhaled by a human body depending on the kind of chemical, 
calculation is effected while taking account of Fugacity 
(f 2 ) of the large particles as well. 

The floor residual amount of the compound is computed 
when Fugacity (f 5 ) of the floor is multiplied by Fugacity 
capacity (Z 5 ). In the case where the solvent is water, in 
particular, it is preferably computed by summing up the 
above-mentioned computed value (f 5 , Z 5 ) and the value 
obtained when Fugacity (f c ) of the space portion is 
multiplied by Fugacity capacity (Z c ). 

(iv) Calculation of Estimated Exposure Amount in Inhalation 
and Inhalation Safety Coefficient (S6 and S7) 

The above-mentioned temporal concentration in the 
indoor air forms a curve shown in Fig. 3A, for example. This 
concentration curve is integrated, an accumulated 
concentration of the compound during a specific period (t x 
to t 2 ) is determined (see Fig. 3B), and the mean concentration 
in the indoor air is computed from thus determined value. 
While an arbitrary period is set as the specific period 
depending on the object, an appropriate period is usually 
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set in view of the method of use of the product and the test 
period of toxicity data. 

Then, according to the above-mentioned mean 
concentration in the indoor air, amount of respiration, and 
exposure time, the estimated exposure amount in inhalation 
is determined. Namely, calculation of: 

estimated exposure amount in inhalation [mg/kg/day] 
= mean concentration in indoor air [mg/m 3 ] x amount of 
respiration [mVkg/min] x exposure time [min/day] 
is effected. Here, as the above-mentioned amount of 
respiration, a published value or actually measured value 
may be used. Also, when amounts of respiration are 
respectively set for adult and child, more appropriate 
estimated exposure amounts in inhalation can be obtained. 
In the case where the inhaled harmful material is not totally 
absorbed into the body but is partially discharged by 
respiration, a more appropriate estimated exposure amount 
in inhalation can be obtained when the inhalation ratio is 
taken into account. 

The inhalation safety coefficient is computed from a 
non-influential amount concerning inhalation toxicity 
examined by an animal experiment beforehand and the 
estimated exposure amount in inhalation determined above. 
Namely, it is expressed as: 

inhalation safety coefficient = 

inhalation non-influential amount [mg/kg/day]/ 
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estimated exposure amount in inhalation [mg/kg/day] 
(v) Calculation of Estimated Percutaneous Exposure Amount 
and Percutaneous Safety Coefficient (S9 and Sll) 

The above-mentioned floor residual amount forms a 
curve such as that shown in Fig. 4A, for example. This 
residual amount curve is integrated, the accumulated 
residual amount of the compound during a specific period 
(t x to t 2 ) is determined (see Fig. 4B) , and the mean residual 
,Q amount is computed from thus determined value. While an 

□ 10 arbitrary period is set as the specific period depending 
ij on the object, an appropriate period is usually set in view 
in of the method of use of the product and the test period of 

□ toxicity data. 

•4 Then, according to the mean floor residual amount, 

15 skin attachment ratio, contact area, and body weight, the 

estimated percutaneous exposure amount is determined. 
Namely, calculation of: 

estimated percutaneous exposure amount [mg/kg/day] = 
(mean floor residual amount [mg/m 2 ] * skin attachment ratio 
20 [%] x contact area [mVday ] ) /body weight [kg] 

is effected. Here, as the contact area, a published value 
(e.g., 4 [m 2 /day]) may be used. The skin attachment ratio 
is a ratio of the compound attaching to the skin when the 
latter is in contact with the floor where the compound exists . 
25 As this value, a published value or a value experimentally 

obtained from a model may be used. 
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A model experiment method for the skin attachment 
ratio is as follows. A weight (8 cm x 8 cm * 8 cm; 4.2 kg) 
is placed on a denim cloth (8 cm * 10 cm) with a pressure 
similar to that of an infant in contact with a floor, and 
the denim cloth is pulled on the floor at a speed (120 cm/15 
sec) similar to the moving speed of the infant. The denim 
and floor are analyzed so as to compute the compound 
contained in the denim and floor. From the ratio 
therebetween, the skin attachment ratio is obtained. It has 
been confirmed that the skin attachment ratio obtained by 
this method is identical to or slightly higher than that 
determined from analyzed values of a hand and a floor when 
the hand is actually pressed against the floor, thereby 
proving this model experiment method to be useful for 
evaluating exposure of inhabitants. 

The percutaneous safety coefficient is computed from 
the non-influential amount concerning percutaneous 
toxicity examined by an animal experiment beforehand and 
the estimated percutaneous exposure amount determined above. 
Namely, it is expressed as: 

percutaneous safety coefficient = 

percutaneous non-influential amount [mg/kg/day]/ 
estimated percutaneous exposure amount [mg/kg/day] 

Nevertheless, in general, percutaneous non- 
influential amount has not often been determined, and there 
are not many published values . Accordingly, a more accurate 
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value can be determined from the estimated percutaneous 
exposure amount, and oral non-influential amount and 
percutaneous absorption ratio for which many published 
values exist, according to the following expression: 

percutaneous safety coefficient = oral non- 
influential amount [mg/kg/day ]/ (estimated percutaneous 
exposure amount [mg/kg/day] * percutaneous absorption ratio 



Here, when the percutaneous absorption ratio is unknown, 
employed is a national guideline (e.g., 10%) which usually 



(vi) Calculation of Estimated Oral Exposure Amount and Oral 
Safety Coefficient (S13 and S14) 

From the estimated percutaneous exposure amount 
obtained in the above-mentioned item (v) , hand surface area 
ratio, and oral transference ratio, the estimated oral 
exposure amount from hand to mouth is determined. Namely, 
calculation of: 

estimated oral exposure amount [mg/kg/day] = 
estimated percutaneous exposure amount [mg/kg/day] * hand 
surface area ratio [%] x oral transference ratio [%] 
is effected. Here, the hand surface area ratio is expressed 
by (hand surface area/body surface area), for which a 
published value (e.g., 5[%]) may be used. The oral 
transference ratio is a hypothetical value, which is set 
to 100%, for example. 



[%]) 
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In the case where oral exposure might occur via 
tableware or food contaminated with the residually sprayed 
compound, it is required that the estimated oral exposure 
amount from tableware or food to mouth be added to the 
estimated oral exposure amount from hand to mouth to yield 
the total estimated oral exposure amount. For example, the 
estimated oral exposure amount from tableware is obtained 
when, according to the tableware residual amount indicating 
the amount of the harmful material remaining in tableware, 
tableware use area which is the sum of tableware surface 
areas, and oral transference ratio from tableware, 
calculation of: 

estimated oral exposure amount [mg/kg/day] = 
tableware residual amount [mg/m 2 ] x tableware use area 
[mVday] x oral transference ratio [ % ] /body weight [kg] 
is effected. Here, the tableware residual amount is 
expressed by (mean floor residual amount x tableware 
contamination ratio ) . As the tableware contamination ratio , 
an actually measured value (e.g. , 9% ) or a hypothetical value 
may be used. 

The oral safety coefficient is computed from the 
non-influential amount concerning oral toxicity examined 
by an animal experiment beforehand and the estimated oral 
exposure amount determined above. Namely, it is expressed 
as : 

oral safety coefficient = 
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oral non-influential amount [mg/kg/day ] / estimated 
oral exposure amount [mg/kg/day] 

As mentioned in the foregoing, in the method of 
estimating an indoor behavior of a pesticidal compound in 
this embodiment, when a chemical containing the above- 
mentioned compound is sprayed over the whole floor area, 
differential equations concerning Fugacity of the space 
portion between carpet fibers, two kinds of large and small 
particles, indoor air, whole carpet fiber, wall, and ceiling 
are simultaneously formed and are solved, and the indoor 
behavior of the compound is estimated according to thus 
obtained solution. Here, the estimation nick time width is 
variably set, while constantly confirming mass balance of 
the compound indoors after the spraying, so that the amount 
of input of the chemical indoors and the resulting solution 
coincide with each other. 

Accordingly, since mass balance of the compound after 
the spraying is always monitored such that the amount of 
input of the compound indoors and the resulting solution 
coincide with each other, thereby varying the estimation 
nick time width; the estimation nick time width is set longer 
when the mass balance fluctuates a little, whereas it is 
set shorter when the mass balance starts fluctuating greatly . 
Namely, when solving simultaneous differential equations 
including a parameter accompanying temporal change, the 
estimation nick time width is automatically set in response 
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to the fluctuation in mass balance. Consequently, when 
processed by a computer, an accurate solution can be obtained 
in a short time. 

The method of evaluating safety of a pesticidal 
compound in accordance with this embodiment uses the 
estimated result mentioned above to evaluate the safety of 
the compound with respect to the human body when the chemical 
is sprayed over the whole floor surface. 

Accordingly, the safety of the compound with respect 
to the human body can be evaluated accurately in a short 
time. As a consequence, when formulating a chemical such 
as insecticide containing the compound, simulation can be 
easily repeated while changing conditions, thereby making 
it easier to formulate a chemical having a high safety 
conforming to the aimed object. 

Though Fugacity is determined by use of Runge- 
Kutta-Gill method in this embodiment, other methods may be 
used for solving differential equations. Runge-Kutta-Gill 
method, however, is preferably used since a program for the 
above-mentioned differential equations can be easily made 
by Basic. Also in the case where differential equations are 
solved by a method other than Runge-Kutta-Gill method, 
similar effects can be obtained when the estimation nick 
time width is set as mentioned above. 

Finally, a computer program product (storage medium) 
employed in each of the above-mentioned embodiments will 




SCC 97-06 



5 



5 10 



151 

."« 15 



20 



25 



be explained. Fig. 19 is a conf igurational view showing a 
storage medium 2 0 employed in each of the above-mentioned 
embodiments of the present invention. As shown in Fig. 19, 
the storage medium 2 0 comprises a program area 20a for 
storing a program. This program area 2 0a stores a 
pesticidal compound indoor behavior simulation program 22. 
The pesticidal compound indoor behavior simulation program 
22 comprises a main program 22a for supervising the 
processing; a differential equation forming program 22b for 
(SI) inputting a primary condition necessary for forming 
a differential equation concerning a f ugacity of the 
compound in each medium in response to an indoor environment 
divided into predetermined media and (S2) calculating a 
secondary condition so as to form the above-mentioned 
differential equation; a Fugacity calculating program 22c 
for (S3) determining the fugacity of the compound from the 
above-mentioned differential equation; an indoor behavior 
calculating program 22d for (S4) computing temporal 
concentration in indoor air and residual amount, which 
constitute an indoor behavior of the compound, from thus 
calculated fugacity; a safety evaluating program 22e for 
(S5 to S17) evaluating the safety of the compound with 
respect to a human body according to the temporal 
concentration in the indoor air and residual amount; and 
a minute time unit changing program 22f for (S21 to S35) 
changing, in response to a fluctuation in mass balance of 
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the compound, the minute time unit used when solving the 

differential equation. The above-mentioned letters and 

numerals within parentheses refer to the corresponding 

steps in Fig. 1 or 5 . 
5 As the storage medium 20, a disk-shaped storage medium 

such as flexible disk or CD-ROM is used, for example. 

Alternatively, a tape-like storage medium such as magnetic 

tape may be used as well. 
k Q The storage medium 2 0 may be used in an information 

q 10 processing apparatus shown in Figs. 20 and 21. Namely, the 

iy information processing apparatus 3 0 comprises a medium 

S drive unit 32, which can accommodate the storage medium 20. 

p As the storage medium 2 0 is thus accommodated, the 

m information stored therein can be accessed by the medium 

^ 15 drive unit 32. Consequently, the pesticidal compound 

5A3f indoor behavior simulation program 22 stored in the program 

area 20a can be executed by the information processing 

apparatus 30. 

This information processing apparatus is configured 
20 as follows. First, it comprises the above-mentioned medium 

drive unit 32, an image memory 34 for storing image data 
indicating simulation results and the like; a work memory 
(internal memory) 3 6 in which an operation system (OS) is 
resident; and a display 38 which is a display means. Also, 
25 it comprises an input device 44 which is an input means having 

a mouse 40 and a keyboard 42 for receiving input data such 
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as the primary condition; a printer 4 6 for outputting image 
data and the like; and a CPU 4 8 for controlling execution 
of the pesticidal compound indoor behavior simulation 
program 22 and the like. 

As the medium drive unit, in conformity to the storage 
medium 20, a flexible disk unit, CD-ROM drive unit, magnetic 
tape drive unit, or the like may be used. 

In the following, the processing of the pesticidal 
compound indoor behavior simulation program 22 stored in 
the program area 2 0a of the storage medium 2 0 will be 
explained. This processing is performed by executing the 
pesticidal compound indoor behavior simulation program 22 
read out by the medium drive unit 32. Upon this execution, 
the main program 2 2a of the pesticidal compound indoor 
behavior simulation program 2 2 is activated. 

The subsequent processing operations of the 
differential equation forming program 22b, Fugacity 
calculating program 22c, indoor behavior calculating 
program 22d, safety evaluating program 22e, and minute time 
unit changing program 22f correspond to those explained in 
the above-mentioned embodiments and Figs. 1 and 5. 

As explained in the foregoing, the computer program 
product (storage medium) employed in each of the above- 
mentioned embodiments of the present invention has a program 
for initially dividing the environment into predetermined 
media (constituents), simultaneously forming differential 
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equations concerning Fugacity of the compound in these 
medium, and then solving them, thereby estimating an indoor 
behavior of a compound when a chemical containing the 
compound is residually sprayed . 
5 Consequently, the environment is assumed to be a 

plurality of media, and exchanges of the compound between 
media and the like are taken into account, whereby simulation 
results close to the actual behavior of the compound can 
hS be obtained. 

O 10 The estimation nick time width set when solving the 

3SSSJ 

simultaneous differential equations is variably set, while 
constantly confirming mass balance of the compound indoors 
g after the residual spraying, so that the amount of input 

p of the chemical indoors and the resulting solution coincide 

if* 

,j3 15 with each other. 

i~f\ 

Accordingly, since mass balance of the compound after 
the residual spraying is always monitored such that the 
amount of input of the compound indoors and the resulting 
solution coincide with each other, thereby varying the 

20 estimation nick time width; the estimation nick time width 

is set longer when the mass balance fluctuates a little, 
whereas it is set shorter when the mass balance starts 
fluctuating greatly. Namely, when solving simultaneous 
differential equations including a parameter accompanying 

25 temporal change, the estimation nick time width is 

automatically set in response to the fluctuation in mass 
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balance. Consequently, when processed by a computer, an 
accurate solution can be obtained in a short time. 

Also, the computer program product (storage medium) 
employed in each of the above-mentioned embodiments of the 
5 present invention has a program for evaluating, by use of 

the above-mentioned estimated results, safety of a 
pesticidal compound with respect to the human body when the 
chemical is residually sprayed. 
k£3 Accordingly, the safety of the compound with respect 

I 10 to the human body can be evaluated accurately in a short 

time. As a consequence, when formulating a chemical such 
as insecticide containing the compound, simulation can be 
easily repeated while changing conditions, thereby making 
it easier to formulate a chemical having a high safety 
15 conforming to the aimed object. 

Industrial Applicability 

In the method of simulating an indoor behavior of a 
pesticidal compound in accordance with the present 
invention, the indoor environment is divided into 
2 0 predetermined media, and exchanges of the compound between 

media and the like are taken into account, whereby simulation 
results similar to the actual behavior of the compound can 
be obtained. Also, when solving simultaneous differential 
equations including a parameter accompanying a temporal 
25 change, a minute time unit can be automatically set in 

response to the fluctuation in mass balance, whereby an 
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accurate solution can be obtained in a short time when 
processed by a computer. 

The method of simulating an indoor behavior of a 
pesticidal compound in accordance with the present 
invention further comprises a step of evaluating, according 
to the indoor behavior of the compound, safety of the 
compound with respect to a human body, whereby the safety 
of the pesticidal compound can be evaluated accurately in 
a short time- Accordingly, when formulating a chemical such 
as insecticide containing the compound, simulation can be 
easily repeated while changing conditions, thereby making 
it easier to formulate a chemical having a high safety 
conforming to the aimed object. 



